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directly  in  the  overall  Shuttle  model.  In  this  case,  a very  simplified  model  was  chosen 
to  represent  the  critical  characteristics  (modes,  etc. ) and  forced  to  match  the  detailed 
model  shown.  The  degrees  of  freedom  using  this  simplified  model  were  reduced 
several  orders  of  magnitude  without  loss  of  accuracy  in  the  overall  systems  model. 
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0 COMPLETE  CARGO  MANIFEST  IS  ANALYZED 

0 RESULTS  ARE  PROVIDED  TO  CARGO  ORGANIZATION  (INCLUDES  LOAD  TRANSFORMATION 
MATRIX  DATA  FOR  CRITICAL  STRUCTURE) 


LOADS/DEFLECTIONS  REQUIREMENTS  AND  ANALYSIS  RESPONSIBIL'TIES 
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SPACECRAFT  LOADS  - THE  LOAD' CYCLE  APPROACH 


The  subject  of  this  presentation  is  the  approach  to  spacecraft  loads  deter- 
mination which  the  Aerospace  Corporation  recommends  to  the  Air  Force  for 
most  of  its  space  programs.  We  call  it  the  load  cycle  approach  because  it  occurs 
in  iterative  cycles,,  each  cycle  being  structured  to  accommodate  the  needs  of  the 
program  as  the  spacecraft  design  evolves. 

The  steps  enumerated  on  the  chart  call  first  for  the  establishment  of  a set 
of  design  load  factors  to  which  the  spacecraft  structure  is  initially  sized.  Sub- 
sequently, three  distinct  cycles  of  loads  analysis  are  called  for:  The  Preliminary, 

Final  Design  and  Verification  load  cycles.  These  cycles  are  distinguished  from 
each  other  largely  by  the  spacecraft  dynamic  models  employed  for  each  one.  The 
Preliminary  and  Final  Design  models  are  derived  from  finite  element  structural 
models  and  mass  data  appropriate  to  know  lege  of  the  design  at  the  time  of  the 
analysis.  The  Verification  Load  Cycle  employs  an  experimentally  derived 
dynamic  model  based  on  a modal  survey  of  the  actual  spacecraft. 

The  time  phasing  of  the  load  cycles  is  indicated  on  the  chart  by  noting  that 
the  Preliminary  Load  Cycle  should  be  completed  by  the  time  of  the  programme 
Preliminary  Design  Review  (PDR);  the  Final  Design  Load  Cycle  should  be  com- 
pleted at  the  time  of  the  Critical  Design  Review  (CDR);  and  the  Verification  Load 
Cycle  should  be  completed  as  long  before  the  first  flight  as  possible. 

The  number  and  variety  of  load  cycles  described  on  the  chart  are  typical 
of  those  performed  for  a program  of  average  complexity.  Variations  can  occur 
on  any  particular  program.  A typical  variation  consists  of  conducting  more  than 
the  specified  three  load  cycles  for  a program  which  undergoes  significant  design 
changes  after  program  initiation.  In  other  instances,  considerations  of  program 
cost  may  dictate  a reduction  in  the  number  of  load  cycles  by  eliminating  the 
Preliminary  Load  Cycle.  When  this  latter  variation  occurs  there  is  a corespon- 
ding increase  in  the  risk  of  expensive  redesigns  late  in  the  program. 
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THE  LOAD  CYCLE  PROCESS 
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This  chart  shows,  in  flow  diagram  form,  the  iterative  nature  of  the  load 
cycle  process.  The  Preliminary  Load  Cycle  results  in  loads  which  may  call 
for  redesign  of  some  structure.  These  changes,  along  with  any  others  resulting 
from  the  maturing  of  the  design  are  incorporated  into  a finite  element  structural 
model  of  what  is  hopefully  the  final  design.  This  model  is  used  in  the  FincJ 
Design  Load  Cycle.  The  loads  resulting  from  this  load  cycle  are,  character- 
istically, the  last  set  of  loads  available  be.fore  the  commitment  is  made  to 
build  the  spacecraft  hardware.  Changes  made  to  the  design  after  this  time 
are  therefore  extremely  costly.  In  addition,  the  pace  of  a typif'aJ  spacecraft 
program  is  usually  one  which  makes  it  necessary  to  use  the  results  of  the  Final 
Design  Load  Cycle  as  the  basis  for  the  static  test.  The  modal  survey  provides 
the  basis  for  the  spacecraft  dynamic  model  used  in  the  Verification  Load  Cycle. 
The  results  of  this  load  cycle  are  shown  being  compared  with  the  stetic  test 

\ 

loads  to  form  the  final  assessment  of  structura.1  qualification  and  approval  for 
flight. 
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LOAD  ANALYSIS  INTERFACES 
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LOAD  ANALYSIS  INTERFACES 


A number  of  contractors  contribute  to  the  performance  of  a load  analysis. 
The  chart  depicts  the  interaction  between  these  contractors. 

The  loads  analysis  itself  is  performed  by  the  Loads  Contractor  shown  in 
the  center  of  the  chart.  This  activity  can  be  performed  by  the  launch  vehicle 
contractor,  an  upper  stage  contractor,  or  by  any  organization  possessing  the 
basic  analysis  capability.  The  contract  for  his  work  in  this  regard  should  be 
written  in  a manner  which  distinguishes  it  from  any  other  supporting  activity 
performed  for  the  spacecraft  program. 

Basic  data  inputs  to  the  Loads  Contractor  are  made  by  the  Spacecraft, 
Launch  Vehicle,  and  Propulsion  Contractors.  These  data  are  identified  on  the 
chart  as  the  spacecraft  and  launch  vehicle  dynamic  models,  loads  transforma- 
tion matrices  and  the  forcing  functions.  Each  of  these  items  is  elaborated  upon 
in  subsequent  charts. 

The  Aerospace  Co rpo ratio"  is  indicated  on  the  chart  acting  in  support  of 
the  Air  Force,  reviewing  and  validating  the  data  inputs  and  the  loads  results. 

The  Air  Force  Space  and  Missile  Systems  Organization  (SAMSO)  has  issued  a 
Commander's  Policy  which  calls  for  the  independent  validation  of  the  last 
(Verification)  loads  analysis  performed  for  each  SAMSO  program,  Although  we 
are  not  SAMSO's  exclusive  agent  for  this  activity,  most  of  the  independent  loads 
verification  analyses  performed  to  date  have  been  conducted  by  the  Aerospace 
Corporation. 


PRELIMINARY  DESIGN  LOAD  FACTORS 
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PRELIMINARY  DESICJN  LOAD  FACTORS 


Preliminary  design  load  factors  are  the  subject  of  this  chart.  Exper- 
ience has  shown  that  care  taken  with  the  development  of  these  quantities  pays 
significant  dividends  late  in  the  program.  Along  with  stiffness  and  clearance 
requirements,  they  form  the  basis  for  the  general  structural  design  of  the 
spacecraft  and,  as  such,  have  a primary  influence  on  its  dynamic  character- 
istics and  strength.  Preliminary  Design  Load  Factors,  in  the  present  context, 
are  taken  to  mean  a specified  set  of  accelerations  which  can  be  used  to  design  the 
spacecraft  primary  structure,  the  spacecraft  adapter,  the  components  and  the 
local  structure  which  supports  the  components. 

Typically,  two  sets  of  Preliminary  Design  Load  Factors  are  prescribed. 
One  set  pertains  to  the  gross  accelerations  expected  to  be  experienced  by  space- 
craft during  the  various  transient  events  which  will  take  place  during  launch  and 
ascent.  These  gross  accelerations  form  the  basis  for  design  of  the  primary 
structure.  A second  set  of  accelerations,  which  apply  to  expected  local  accel- 
erations, is  also  prescribed  and  these  are  used  in  the  design  of  components 
and  their  supporting  structure.  The  selection  of  these  compvj-r^^t  load  factors 
is  an  area  in  which  conservatism  pays  particularly  rich  dividends.  The 
unexpected  high  accelerations  which  may  be  discovered  in  the  later  load  cycles 
are  veiy  often  local  in  nature  and  preparation  for  this,  in  terms  of  conservative 
component  design  load  factors,  can  often  prevent  costly  redesign  problems  late 
in  the  program. 

The  general  importance  of  conservatism  in  selecting  Preliminary  Design 
Load  Factors  cannot  be  overstressed.  As  noted  above,  they  have  a primary 
influence  on  the  general  structural  arrangement  selected  by  the  designer.  The 
need  to  cope,  from  the  outset,  with  significant  loads  can  force  the  designer  to 
make  structural  efficiency  a major  objective  in  his  overall  design.  Too  often, 
optimistic  initial  estimates  of  load  have  given  him  the  impression  that  he  can 
compromise  strength  considerations  in  favor  of  hi.s  many  other  mission  objectives. 
Later,  when  the  design  is  modeled  and  higher  loads  are  calculated,  he  finds  that 
additional  strength  can  only  be  obtained  at  the  expense  of  significant  weight, 
schedule  and  cost  impact. 
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PRELIMINARY  AND  FINAL 
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STATISTICAL  TREATMENTS  - "THREE  SIGMA"  LIMIT  LOADS 
"WORST  CASP'  TREATMENTS 
UNCERTAINTY  FACTORS 


PRELIMINARY  AND  FINAL 


DESIGN  LOAD  CYCLE 

After  the  spacecraft  structure  has  been  designed  and  sized  in  accordance 
with  the  Preliminary  Design  Load  Factors  and  when  the  other  functional  aspects 
of  the  design  have  been  established,  a series  of  load  cycles  is  initiated  based 
upon  analytical  models  of  the  resulting  design.  This  chart  lists  some  of  the  key 
ingredients  in  these  load  cycles. 

Analytical  Modeling.  Analytical  modeling  forms  the  basis  for  these  load  cycles 
and,  as  such,  is  an  essential  ingredient.  Analytical  models  of  both  the  space- 
craft and  the  launch  vehicle  are  required  and  these  models  must  faithfully 
reproduce  system  dynamic  characteristics  throughout  the  frequency  range  in 
which  significant  structural  loads  may  occur.  The  range  up  to  about  50  Hz 
has  been  found  to  be  important  but,  on  occasion,  higher  frequency  fidelity  is 
required.  Most  of  the  critical  transient  events  involve  significant  longitudinal 
excitation  and  the  models  must  therefore  reproduce  at  least  the  first  longitudinal 
mode  of  the  coupled  launch  vehicle/  spacecraft  system.  The  frequency  of  this 
mode,  particularly  for  upper  stage  events,  often  exceeds  50  Hz. 

Achievement  of  the  required  modeling  fidelity  requires  careful  treatment 
of  the  stiffness,  mass  and  damping  properties  of  the  system.  Extensive  finite 
element  modeling  of  the  structure  is  usually  required  and  experience  has  shown 
that  attention  to  such  details  as  separation  joints  and  bearings  is  extremely 
important  and  experimental  data  are  often  invaluable  in  these  instances.  The 
acquisition  of  configuration  peculiar  experimental  data  for  the  spacecraft 
normally  occurs  after  the  Preliminary  Load  Cycle  is  well  under  way  and  the 
use  of  test-derived  information  is  often  an  important  distinguishing  attribute  of 
the  Final  Design  Load  Cycle.  Spacecraft  mass  properties  usually  undergo 
significant  evolution  between  these  two  load  cycles  as  well.  The  spacecraft 
damping  characteristics  used  in  the  Preliminary  and  Final  Design  Load  Cycles 
are  invariably  simple  estimates  based  on  experience.  Mode  surveys  of  many 
different  spacecraft  have  shown  that  modal  viscous  damping  coefficients  on  the 
order  of  one  percent  of  critical  are  the  rule.  The  precise  magnitude  of  this 
damping  has  not  been  found  to  be  a critical  parameter  for  most  spacecraft 
transient  loading  events. 

The  dynamic  model  of  the  launch  vehicle  which  is  used  in  these  load  cycles 
(and  in  the  Verification  Load  Cycle  as  well)  is  often  a very  mature  one.  Knowledge 
of  launch  vehicle  characteristics  evolves  as  flight  experience  is  gained  and  models 
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for  the  standard  boosters  are  reasonably  well  established  today.  While  it 
has  been  found  that  many  of  the  beam -like  structural  charactertistics  of 
expendable  launch  vehicles  can  be  accurately  derived  by  relatively  simple 
analytical  means,  other  features  require  extensive  modeling  work.  For 
example,  experience  has  shown  that  analytical  prediction  of  the  dynamic 
longitudinal  behavior  of  a launch  vehicle  requires  careful  treatment  of  the 
fluid -structural  interactions  which  take  place  in  the  propellant  tanks.  The 
engine  support  structure  is  usually  another  structurally  complex  system  which 
requires  extensive  study  to  represent  adequately.  Often,  both  the  engine  and 
tankage  behavior  are  only  treated  properly  after  flight  measured  structural 
response  data  are  available  for  examination  and  the  system  has  been  remodeled. 

The  spacecraft  and  launch  vehicle  are  modeled  separately  by  their 
respective  contractors  and  analytically  coupled  to  each  other  by  the  loads 
contractor.  This  operation  requires  the  use  of  modal  synthesis  techniques 
which  take  account  of  all  significant  structural  aspects  of  the  spacecraft-to- 
launch  vehicle  interface. 

Treatment  of  the  launch  event  requires  analytical  modeling,  not  only  of 
the  spacecraft/launch  vehicle  system,  but  also  of  the  launch  stand.  Again, 
details  play  an  important  role.  The  nature  of  the  interface  between  the  stand 
and  launch  vehicle  must  be  accurately  accounted  for  in  the  model. 

Finally,  the  modeling  activity,  in  addition  to  providing  an  accurate 
representation  of  the  system  dynamic  characteristics  must  also  result  in 
load  transformation  matrices  (IjTM's)  which  relate  mass  point  inertial  loads 
to  internal  member  forces.  In  order  to  develop  an  LTM  it  is,  of  course, 
essential  that  the  acceleration  forces  applied  at  mass  points  in  the  dynamic 
model  be  relatable  to  node  point  loads  in  the  structural  model.  For  this 
reason  it  is  ideal  that  the  dynamic  model  be  derived  directly  from  the  model 
to  be  used  for  structural  analysis.  For  this  reason  and  also  because  the  LTM 
forms  the  basis  for  deriving  the  static  test  loads,  it  has  been  found  that  the 
structural  portion  of  the  modeling  effort  is  best  performed  by  structures 
personnel  in  a contractor's  organization. 
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Forcing  Functions.  Most  of  the  critical  flight  events  involve  engine  thrust 
transients  - ignitions  or  shutdowns.  It  has  been  found  that  spacecraft 
transient  response  is  highly  sensitive  to  the  detailed  character  of  thrust 
transients  and  the  only  reliable  source  of  such  detailed  information  is  high 
sample  rate  flight  data.  In  the  absence  of  such  data,  the  results  of  ground 
test  firings  can  sometimes  be  used,  but  usually  with  some  compromise  of 
confidence  in  the  details.  The  use  of  analytical  thrust-time  predictions  is 
limited  to  special  instances  and  "specification”  type  thrupc  transients  are 
rarely  satisfactory.  Propulsion  personnel  can  play  an  important  role  in 
certifying  the  validity  of  any  thrust  transient  used  as  a forcing  function. 

The  engine  ignition  event  can  be  accompanied  by  a "pressure  pulse" 
phenomenon  caused  by  reflection  of  exhaust  cases  from  various  parts  of  the 
launch  pad  facility.  The  forcing  functions  employed  to  date  to  simulate  the 
effects  of  this  phenomenon  have  been  derived  from  limited  experimental  data 
and  tailored  to  give  reasonable  agreement  with  measured  spacecraft  or  launch 
vehicle  response.  Scale  model  flow"  testing  has  been  employed  as  an  aid 

in  establishing  some  aspects  of  the  phenomenon. 

Gust  induced  forces  are  characteristically  treated  by  determining  the 
response  of  the  launch  vehicle/ spacecraft  system  to  a discrete  gust  having  an 
idealized  shape  which  is  "tuned"  over  a range  of  wavelengths. 

Some  launch  vehicle  configurations  are  subject  to  the  phenomenon  of 
periodic  vortex  shedding  in  the  presence  of  steady  ground  winds.  This 
phenomenon  poses  a potential  load  problem  for  the  launch  vehicle  and  can  be 
of  significance  for  the  spacecraft  in  that  it  establishes  initial  conditions  for 
the  launch  event.  Wind  tunnel  data  provides  the  only  source  of  information 
available  for  treating  this  problem. 

The  problem  of  establishing  forcing  functions  to  represent  buffeting  is 
a frustrating  one.  A sizable  amount  of  experimental  work  was  done  on  the 
subject  during  the  early  days  of  the  space  program.  "Hammerhead"  shapes 
were  investigated  in  wind  tunnels  and  a considerable  volume  of  data  were 
collected  for  specific  shapes.  Unfortunately  the  data  is  highly  specialized  and 
its  application  to  new  shapes  is  always  questionable.  Fortunately,  acceleration 
response  data  for  the  cylinder-cone  shapes  characteriscic  of  most  current 
configurations  show  only  modest  response  due  to  buffeting.  Data  from  Titan  III 
flights,  however,  indicates  response  during  transonic  flight  which  has  been 
attributed  to  the  effects  of  interaction  between  the  large  solid  rocket  motors 
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and  the  launch  vehicle  center  body.  Since  the  only  data  relevant  to  the  pheno- 
menon are  flight  response  measurements,  a procedure  has  been  developed  for 
scaling  these  flight  data  for  application  to  new  configurations.  It  involves 
identifying  "similar"  system  modes  for  the  flight  tested  and  new  spacecraft 
systems.  The  procedure  is  not  completely  satisfying  and  therefore  requires 
the  use  of  a safety  factor. 

Transient  Response  Analyses.  After  the  dynamic  models  and  forcing  functions 
have  been  assembled  the  remaining  task  is  determining  the  transient  response 
information  which  is  needed  to  verify  the  spacecraft  design.  The  chart  lists 
several  considerations  which  bear  on  this  design  load  issue.  It  was  noted 
above  that  flight  data  should  be  utilized  to  obtain  thrust  transients.  These 
transients  vary  from  flight  to  flight  and  a sizeable  family  of  them  is  necessary 
to  form  a representative  ensemble.  Analyses  are  conducted  using  each  thrust 
transient  and  accelerations  for  each  mass  point,  along  with  forces  for  each 
critical  spacecraft  member  are  calculated.  Examination  of  these  responses 
shows  significant  variability  from  case  to  case  and  statistical  treatment  is 
required  to  establish  bounds  for  the  design  loads.  Design  limit  load  is  usually 
taken  as  the  "mean-plus -three -sigma"  value  of  the  calculated  response.  Used 
in  this  sense,  mean -plus -three-sigma  is  meant  to  imply  a probability  of  non- 
exceedance comparable  to  the  three  sigma  level  for  a normally  distributed 
statistical  variable. 

In  instances  where  only  limited  amounts  of  forcing  function  data  are 
available  the  small  size  of  the  sample  should  be  accounted  for.  In  instances 
where  very  small  amounts  of  forcing  function  data  are  available,  statistical 
treatment  is  not  possible.  In  this  event,  v/orst  case  responses  are  multiplied 
by  a factor  to  obtain  limit  load. 

Confidence  in  the  spacecraft  analytical  model  should  increase  as  the 
program  matures.  In  the  early  stages  of  the  program,  when  the  Preliminary 
Load  Cycle  model  is  developed,  knowledge  of  the  structural  and  mass  property 
details  of  the  design  is  necessarily  limited.  It  is  recommended  practice  to 
account  for  this  uncertainty  by  applying  a multiplying  factor  of  ).  5 to  the 
calculated  responses  from  this  load  cycle.  When  the  analytical  model  for  the 
Final  Analytical  Load  Cycle  is  developed,  knowledge  of  the  design  should  be 
significantly  improved  and  this  "uncertainty  factor"  is  usually  reduced  to  1.  25. 


or 

a. 

o 

a: 

a.  z 
D_  o 
< r- 


= 8 


S <x 

Qi  O 

LU  U. 

tu  LU 
o CJ> 


— o 

or  ~ 

LU  uo 

a. 

X LU 

Q o 


§ o z 

< z:  o 
bi  < — 

S □□  tu 

(/)*“  = 
LU  to  LO 
0^00  . 

ZD  _J  t— 

Of 


O < 

' — I 

Q ZD 


LU 

> 

>- 

1 

O 

o 

MM 

1 

_l 

LlJ 

H- 

H- 

<r 

< 

X 

</) 

o 

o 

O 

h— 

o 

mJ 

CO 

ZD 

s 

Qi 

LL. 

O 

• 

• 

< 

• 

1— 

h— 

oo 

00 

>OJ 

UJ 

OO  tj 

^ >-  Q 

LO 

U-  P 

^ s-  s 

< 9r  — t/^ 


— z:  00 


R o 

o o 


< < 
S Q 

ZD  Z 

g 8 

LO  00 


0 LU 

1 I V) 
=>< 
00  S 


LU 

> 

uu 

a- 

o 

> 

LU 

h- 

oo 

-J 

& 

< 
1 — 

00 

1— 

z 

h- 

oo 

n 

2 

LU 

< 

LU 

z 
1 1 1 

< 

z 

< 

>- 

o 

o 

oo 

ZD 

LU 

LU 
— J 

o 

i 

1 

Qi 

F— 

0^ 

OL 

o 

x: 

00 

o 

h~ 

Z 

O 

o 

O 

O 

• 

z 

• 

• 

'1 


61 


THE  SPACECRAFT  MODE  SURVEY 

The  commitment  to  construct  the  spacecraft  is  made  at  the  CDR  follow- 
ing completion  of  the  Final  Design  Load  Cycle.  This  chart  addresses  the 
spacecraft  mode  survey  which  occurs  as  soon  as  a suitable  test  article  can  be 
built  and  made  available  for  testing. 

Ebcperience  has  shov/n  that  a spacecraft  often  experiences  loads  which 
depend  upon  the  detailed  character  of  its  higher  order  modes,  modes  in  the  20 
to  50  Hz  range.  Further,  it  has  been  found  that  modifying  analytical  models 
to  give  good  agreement  with  modal  test  results  in  this  range  is  extremely  , 
difficult.  As  a consequence,  the  objective  of  the  modal  survey  is  the  measure- 
ment of  all  modes, in  the  frequency  range  of  interest,with  sufficient  accuracy  to 
permit  their  direct  use  in  the  Verification  Load  Cycle. 

The  test  article  should  be  of  flight  quality  in  so  far  as  its  dynamic  charac- 
teristics are  concerned.  Some  simulation  of  sensitive  components  is  usually 
necessary  but  the  mass  and  stiffness  properties  of  the  simulators  must  accurately 
replicate  those  of  the  real  equipment.  Black  boxes  and  solar  cells  are  typical 
candidates  for  simulation.  Occasionally,  it  may  be  desirable  to  remove  certain 
subsystems  from  the  test  article  to  simplify  the  testing  process.  Solar  arrays 
are  a typical  example.  When  this  is  done,  the  subsystem  modes  are  determined 
in  a separate  test  and  coupled  analytically  to  the  measured  spacecraft  normal 
modes.  The  spacecraft  test  article  must,  of  course,  be  suitably  ballasted  and 
instrumented  so  that  modes  will  be  obtained  which  will  make  this  coupling 
practicable  after  the  test. 

For  spacecraft  which  are  launched  on  expendable  launch  vehicles,  the 
mode  survey  is  usually  conducted  with  the  spacecraft  mounted  on  its  adapter 
with  the  adapter  grounded  at  its  launch  vehicle  interface.  Some  Shuttle - 
launched  spacecraft  will  require  mode  testing  while  mounted  in  their  "cradles”. 

When  the  test  is  completed  but  before  the  set-up  is  torn  down,  the  mea- 
sured results  should  be  subjected  to  as  much  scrutiny  as  possible  to  provide 
assurance  that  an  adequate  set  of  data  have  been  acquired.  Completeness  of  the 
set  of  measured  modes,  in  the  frequency  range  of  interest,  is  judged  by  review- 
ing all  of  the  sinusoidal  sweep  data  obtained  during  the  test.  All  indications  of 
modal  presence  in  the  sweep  data  must  be  accounted  for  in  the  set  of  modes  which 
were  measured.  Questions  concerning  missing  modes  must  all  be  resolved.  The 
quality  of  the  measured  modes  themselves  can  be  judged  by  a number  of  criteria. 
The  precision  of  the  "tuning"  achieved  during  mode  acquisition  is  a matter  of 
judgement  exercised  by  the  test  engineer  during  the  test.  After  the  set  of  modes 
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has  been  recorded,  however,  an  orthogonality  check  is  made  to  test  the 
independence  of  the  measured  modes  with  respect  to  each  other.  An 
analytically  derived  mass  matrix  is  necessary  for  perfojrmance  of  this 
orthogonality  check.  The  matrix  should  be  derived  in  a manner  which  is 
compatible  with  the  instrumentation  employed  during  the  test.  The 
orthogonality  criterion  applied  to  the  measured  modes  requires  that  all  off- 
diagonal  terms  in  the  generalized  mass  matrix  [m]  be  less  than  0.  10,  where 

[Kl]  = C<^]  ^ [m]  [♦] 

and  [♦]  = matrix  of  measured  modes  normalized  to  unity  for  each 
modal  generalized  mass 

[m]  = test  article  mass  matrix 
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OPTIMIZED  TESTING/ DATA  REDUCTION  TIME 
ENHANCED  CLOSE  MODE  DISCRIMINATION 


' MODE  SURVEY  TESTING  METHODS 

The  experimental  determination  of  the  normal  modes  of  a complex 
spacecraft  is  a technically  demanding  problem.  This  chart  lists  some  key 
features  of  current  mode  testing  methods  and  comments  on  the  current  state 
of  the  testing  technology. 

The  conventional  method  v/as,  until  relatively  recently,  the  only 
1 technique  employed  for  major  mode  surveys.  It  consists  of  applying  sinusoidal 

excitation  to  the  test  article  using  several  shakers  which  are  physically  l<>cated 
and  driven  in  a manner  which  excites  only  one  mode  at  a time.  This  technique, 
while  difficult  to  execute,  has  the  advantage  that  the  modes  are  individually 
excited  and  can  be  subjected  to  close  quantitative  and  qualitative  examination, 
one  at  a time.  It  has  the  further  advantage  that  data  reduction  is  substantially 
! completed  as  each  mode  is  surveyed,  so  that  a complete  set  of  checks  can  be 

made  before  the  test  set-up  is  torn  dovm.  The  technique  has  been  employed 
with  considerable  success  for  many  spacecraft  programs. 

The  conventional  technique  has  two  significant  shortcomings*. 

(a)  its  use  places  major  demands  on  the  skill  and  perceptiveness 

of  the  test  engineer,  particularly  in  the  matter  of  discriminating 
modes  which  are  closely  spaced  in  frequency,  and 

(b)  a considerable  amount  of  time  is  required  for  its  execution. 

The  time  requirement  is  a particularly  import-  r»t  disadvantage  since  conduct  of 
the  test  requires  the  exclusive  use  of  a test  article  which  must  be  released  at 
the  earliest  possible  time  for  other  testing  purposes. 

A secc.id  technique,  called  here  the  single  point  random  (SPR)  method, 

' has  received  considerable  attention  of  late.  It  makes  use  of  an  approach 

which  is  entirely  different  from  the  conventional  one.  Instead  of  exciting  one 
I mode  at  a time,  the  method  calls  for  exciting  several  modes  at  once  using  a 

single  shaker  which  delivers  random  force  excitation.  The  measured  excitation 
and  response  are  then  analytically  processed  to  obtain  transfer  functions  for 

\ 

each  measurement  station  on  the  jpacecraft.  Modal  parameters  are  then 
[ deduced  by  further  analytical  processing  of  these  transfer  functions.  The 

t 

test  is  relatively  simple  to  pei*furm.  In  principle,  shaker  locations  can 
> be  selected  before  the  test  and  application  of  the  excitation  is  a routine  process 

I which  can  be  performed  in  a short  period  of  time.  The  data  reduction  process 

I can  be  performed  after  the  testing  itself  (i.  e.  the  excitation)  has  been  completed. 

y 
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The  test  set-up  can  therefore  be  torn  down  at  an  early  date  and  the  test  article 
can  be  released  for  other  purposes. 

At  the  present  time  the  SPR  method  can  be  said  tc  have  two  significant 
shortcomings: 

(a)  Completion  of  th^  data  reduction  process  requires  a significant 
amount  of  time  (more  than  a month  for  a complex  spacecraft)  and 

(b)  The  effectiveness  of  the  method  has  not  been  established  under  a 
sufficient  variety  of  circumstances  to  demonstrate  that  difficult 
modal  discrimination  problems  can  be  solved  with  its  use. 

When  SPR  techniques  are  employed,  for  a modal  survey,  the  Aerospace 
Corporation  currently  recommends  that  verification  of  the  results  be  obtained, 
after  the  fact,  using  a conventional  sinusoidal  search  and  dwell  technique. 

It  is  thus  seen  that  improved  modal  testing  technology  is  definitely 
needed.  Research  in  this  field  should  be  aimed  at  improving  existing  techniques 
and  at  the  development  of  new  techniques  as  well.  The  objectives,  as  noted  on 
the  chart,  must  include  optimizing  both  the  testing  and  the  data  reduction  time 
and  must  address  the  problem  of  close  mode  discrimination.  There  is  work 
currently  underway  in  this  area.  Modal  testing  problems  have  been  the  subject 
of  a number  of  projects  which  have  been  pursued  by  the  Air  Force  Flight 
Dynamics  Laboratory,  the  Jet  Propulsion  Laboratory,  the  Langley  Research 
Center  and  by  several  contractors.  Also,  we  at  Aerospace  have  had  the  sub- 
ject under  consideration  for  some  time.  Investigations  by  the  various  organi- 
zations range  from  seeking  improvements  to  existing  methodology  to  the 
development  of  entirely  new  techniques. 

A common  feature  of  most  new  methods  is  the  avoidance  of  exciting 
individual  modes,  one  at  a time.  Instead,  modes  are  deduced  by  analytical  means 
from  forced  or  free  decay  response  measurements.  This  approach  treats  the 
testing  time  problem  admirably  but  leaves  the  mode  discrimination  problem  open 
to  question.  Mode  discrimination  is  a serious  problem  in  the  conventional, 
sinusoidal  sweep  and  dwell  test  as  well,  but  when  modes  are  separately  excited 
and  observed  the  existence  of  the  problem  is  apparent  and  measures  can  be  taken 
to  find  a solution.  The  fundamental  issue,  then  is  how  assurance  can  be  acquired 
that  valid  modes  have  been  obtained  with  the  new  methods  in  the  absence  of 
direct  modal  excitation  and  observation. 
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Aa  noted  later  in  this  presentation,  modal  testing  is  a major  contri- 
butor to  the  cost  of  the  load  cycle  process.  Improvements  which  would 
reduce  this  cosi:  are  sorely  needed  and  it  is  recommended  that  current 
research  to  this  end  be  both  continued  and  enhanced. 


THE  VERIFICATION  LOAD  CYCLE 
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THE  VERIFICATION  LOAD  CYCLE 


The  ultimate  determination  of  spacecraft  loads  occurs  in  the  Verifi- 
cation Load  Cycle.  This  chart  describes  the  salient  features  of  this  analysis. 

The  Verification  Load  Cycle  occurs  as  soon  as  possible  after  the  mode 
survey  has  been  completed  and  the  data  has  been  reduced  to  suitable  form 
for  use  in  a loads  analysis.  This  requires  that  the  modes  of  separately  tested 
subsystems  be  coupled  to  the  measured  spacecraft  modes;  that  any  other  test 
article  deviation  from  the  flight  condition  be  corrected  and  that  the  measured 
modes  be  ortho  go  nalized.  It  has  been  found  that  this  last  step  is  desireable 
to  insure  mathematical  compatibility  of  the  resulting  spacecraft  representation 
with  the  model  of  the  launch  vehicle. 

The  launch  vehicle  model  and  the  forcing  functions  employed  in  the  Veri- 
fication Load  Cycle  should  incorporate  any  data  necessary  to  account  for  know- 
ledge gained  after  the  Final  Design  Load  Cycle.  This  is  particularly  true  of 
any  thrust  transient  data  which  may  have  been  added  to  the  statistical  famM 
as  the  result  of  additional  flights. 

The  load  transformation  matrix  is  often  updated  between  the  Final 
Design  and  Verification  Load  Cycles.  In  particular,  the  timing  of  the  Veri- 
fication Load  Cycle  often  permits  the  use  of  structural  information  acquired 
during  the  static  test. 

An  uncertainty  factor  of  1.  0 can  be  employed  in  the  Verification  Ijoad 
Cycle  if  the  spacecraft  mode  survey  is  considered  adequate. 

The  results  of  the  Verification  Load  Cycle,  being  the  ultimate  load 
determination,  provide  the  criteria  by  which  the  adequacy  of  the  spacecraft 
static  test  can  be  judged.  Comparison  of  the  member  forces  obtained  from 
the  load  cycle  with  those  induced  during  the  static  test  provides  the  basis  for 
this  judgement  and  for  final  commitment  to  flight. 
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• IMPROVED  MODAL  TESTING  TECHNIQUES 


V 


OBSERVATIONS  ON  THE  LOAD  CYCLE  PROCESS 


This  chart  treats  several  final  observations  on  the  Load  Cycle  Process. 
The  intent  here  is  to  note  that  the  process  has  the  potential  of  fulfilling  its 
goals  provided  adequate  attention  is  given  to  its  planning  and  execution.  On  the 
other  hand,  it  is  also  noted  that  the  process,  in  its  present  state  of  development, 
leaves  significant  room  for  improvement. 

The  need  for  improvement  is  most  dramatically  underscored  by  the 
current  high  cost  of  the  process.  The  quoted  cost,  in  the  neighborhood  of 
one  million  dollars,  is  representative  of  a program  which  calls  fcr  three 
load  cycles.  An  approximate  distribution  of  these  costs  is  shown  on  the  chart. 

As  mentioned  earlier,  the  Preliminary  Load  Cycle  is  occasionally 
eliminated  in  the  interests  of  cost  reduction.  From  the  technical  point  of  view, 
this  entails  some  risk  that  excessive  loads  will  be  discovered  late  in  the  program. 
This  risk  can  be  partially  offset,  however,  by  employing  more  conservative 
preliminary  design  load  factors  than  would  otherwise  be  used.  It  is,  of  course 
conceivable  that  the  entire  load  cycle  process  could  be  by-passed  if  sufficiently 
conservative  design  criteria  were  employed.  This  would  require  the  use  of 
conservative  crite  .-ia  for  both  loads  and  stiffness.  A simple  mode  survey  would 
still  be  required  to  demonstrate  achievement  of  the  stiffness  (frequency)  goals. 

The  cost  breakdown  indicates  that  modeling,  together  with  loads  calcula- 
tion (which  also  involves  a significant  amount  of  modeling)  accounts  for  seventy 
percent  of  the  cost.  Clearly,  improvements  in  modeling  technology  could  have 
a significant  impact  on  load  cycle  costs.  The  need  In  the  modeling  area  is  for 
improvement  in  both  accuracy  and  cost.  There  are  strong  indications  that  the 
number  of  modes  required  to  achieve  convergence  in  certain  loads  analyses  can 
be  dramatically  reduced  by  employing  residual  stiffness  techniques.  This  is  an 
example  of  a step  toward  the  goal  of  reducing  costs  while  improving  accuracy  at 
the  same  time. 

Finally,  the  mode  survey  is  seen  to  account  for  nearly  a third  of  tt.e  total 
load  cycle  cost.  Again,  technical  improvement  and  cost  reduction  are  both  im- 
portant goals.  There  are  a number  of  promising  activities  now  in  progress  in 
the  modal  survey  field.  These  have  already  been  alluded  to  and  they  deserve 
continued  pursuit  and  encouragement. 
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The  SPAR  System 


- SPAR  provides  a high-order  language  for  solving  a broad  spectrum  of 
engineering  problems.  The  system  consists  of  an  array  of  independent 
processors,  communicating  automatically  through  a highly  efficient  generatl 
purpose  engineering  data  base  management  system. 

- Advanced  sparse  matrix  solution  methods  provide  low  execution  costs, 
minimal  central  memory  requirements,  and  large  size  capacity;  permitting 
extremely  fine  meshes  to  be  used.  Static,  buckling,  and  vibrational  problems 
in  the  10,000  to  20,000  DOF  range  are  solved  routinely.  Maximum  capacity, 
without  substructuring,  is  typically  in  the  50,000  DOF  range. 

- Highly  effective  in  both  interactive  and  batch  operation.  Major  studies 
generally  are  best  performed  using  a combination  of  both  batch  (RJE)  and 
interactive  runs.  Graphics  terminals,  30  cps  thermal-printer  terminals, 
and  line  printers  all  are  used  in  various  phases  of  typical  applications. 

- All  input  is  free- field.  Extensive  facilities  are  provided  in  areas  such 
as  mesh  generation,  data  checking,  automation  of  all  aspects  of  problem 
definition,  data  base  interrogation,  and  custom  report  generation. 

- Restarting  and  mixed  batch/ Interactive  operation  is  fully  automated.  To 
restart,  the  user  needs  only  to  reassign  the  file,  or  files,  in  which  the 
data  base  resides,  and  resume  execution  as  though  the  prior  run  had  not 
been  terminated. 

- The  data  base  management  system  provides  an  automatic  means  of  communicating 
with  external  programs,  through  ordinary  sequential  files,  for  purposes  of 
obtaining  source  data  or  furnishing  results.  All  Information  produced  by 
the  system  is  accessible  in  this  way. 

- Effects  of  pre-stress  may  be  included  in  all  forms  of  static,  dynamic,  and 
buckling  analysis. 

- The  EIG  eigensolver  directly  solves  sparse  high-order  eigenproblems  without 
using  any  form  of  DOF  condensation.  10,000-t-  DOF  vibrational  and  buckling 
eigenproblems  are  solved  routinely.  Costs  are  very  low,  particularly  when 
modes  from  prior  analyses  are  used  as  initial  approximations. 

Extensive  facilities  for  substructure  and  other  forms  of  Rayleigh-Rltz 
analysis  are  provided  for  use  where  appropriate.  Fully  coupled  400+  DOF 
vibrational  eigenproblems  are  solved  routinely. 

- Extensive  dynamic  analysis  capability  is  provided:  e.g.  transient,  random, 
shock  spectrum,  steady  state.  Any  system  state  quantities  may  be  tracked 
and/or  recorded  in  the  data  base  for  use  in  subsequent  studies.  For  example, 
the  results  of  a complete  vehicle  response  analysis  may  be  recorded  in  the 
data  base,  and  later  used  to  define  payload  base  excitation. 

- Extensive  thermal  element  repertoire,  including  conduction,  convection, 
radiation,  and  mass  transport  elements.  Steady  state  and  transient  analysis 
of  linear  and  nonlinear  problems  is  performed.  Common  utilities,  e.g.  mesh 
generation,  plotting,  data  entry  routines,  are  shared  by  structural  and 
thermal  functions. 
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Example  SPAR  Runstrean 


Load  case  1: 


I Load  case  2: 

^ , Z . 10000. 

/ 


Commands: 


Explanation: 


eXQT  ikB 
STARl  5 

TITLE ’BEAM  EXAMPLE 
JOINT  LOCATIONS 
$ X Y Z 
1 .0  .0  .0 
2 .0  .0  10. 

3 .0  .0  20. 

4 .0  .0  30.  . 

5 .0  .0  40. 

MATERIAL  CONSTANTS 

1 10. +6,  .3,  .101,  1.-4 

BEAM  ORIENTATION 
1 111  1. 

E21  SECTION  PROPERTIES 
TUBE  ■!  2.0,  2.25 

CONSTRAINT  CASE  1 
ZERO  1,2, 3, 4, 5, 6:  1 

eXQT  ELD 
E21 

1,2:  2,3:  3,4;  4,5 
eXQT  TOPO 
€XQT  E 
eXQT  EKS 
eXQT  K 
@XQT  INV 
eXQT  AUS 

ALPHA:  CASE  TITLES 
1 'TRANSVERSE  LOAD 
2 'AXIAL  LOAD 
SYSVEC:  APPLIED  FORCES 
CASE  1:  1=2;  J=5:  1000, 

CASE  2;  1=3:  J=5:  10000. 
eXQT  SSOL 
6XQT  GSF 
eXQT  PSF 
§XQT  VPRT 

PRINT  STATIC  DISPLACEMENTS 
PRINT  STATIC  REACTIONS 
STOP 


$ The  model  has  5 Joints. 

$ Create  data  set  containing  Joint  locations. 


$ Create  table  of  material  constants 
$ E,  nu,  rho,  alpha  for  material  1. 

$ Create  table  of  orientation  reference  data. 

$ Create  table  of  cross  section  properties. 

♦ Tube  #1  radii=  2.0,  2.25. 

$ Create  data  set  defining  constraint  case  1. 

$ Zero  all  6 motion  components  of  Joint  1. 

. Define  all  elements. 

$ Define  all  type  E21  (general  2-node)  elements. 
$ Beams  connect  Joints  1 and  2,  2 and  3,  etc. 

. Analyze  element  interconnection  topology. 

. Analyze  element  geometry. 

. Create  element  stiffness  matrices,  etc. 

. Create  assembled  system  K. 

. Create  factored  system  K. 

. Enter  Arithmetic  Utility  System. 

. Create  data  set  named  "CASE  TITLES" 


$ Create  a data  set  named  "APPLIED  FORCES" 

$ Page  1 (case  1)  of  "APPLIED  FORCES" 

$ Page  2 (case  2)  of  "APPLIED  FORCES" 

. Create  data  sets  containing  Joint  motions. 
. Create  data  sets  containing  stresses. 

. Produce  Printed  display  of  stresses. 


i- 

Exaaple  - Data  Base  Table  of  Contents  ^ 

J 

\ 

J 

Through  the  following  conunandSi  data  base  tables  of  contents  may  be  obtained: 

6XQT  DCU  . DCU  is  the  Data  Complex  Utility  Program  - 

TOC  1 

The  above  commands  cause  a display  of  the  following  type  to  be  produced: 


TABLE  OF  CONTENTS,  LIBRARY  1 


BEAM  EXAMPLE 
SEQ  RR  DATE 

TIME 

E 

R WORDS 

NJ  N1«NJ 

T 

Y 

DATA  SET 
N1  N2 

NAME 

N3 

N4 

Prooesaor  by 
whioh  created 

1 

17 

1 10778 

191116 

0 

18 

1 

18 

0 

JDF1 

BTAB 

1 

8 

TAB 

2 

18 

1 10778 

191116 

0 

5 

5 

5 

0 

JREF 

BTAB 

2 

6 

it 

3 

19 

110778 

191116 

0 

12 

1 

12 

1 

ALTR 

BTAB 

2 

4 

it 

4 

20 

110778 

191116 

0 

18 

1 

18 

4 

NDAL 

0 

0 

it 

5 

21 

110778 

191116 

0 

15 

5 

15 

1 

JLOC 

BTAB 

2 

5 

it 

6 

22 

110778 

191116 

0 

10 

1 

10 

i 

KATC 

BTAB 

2 

2 

it 

7 

23 

110778 

191116 

0 

5 

1 

5 

1 

MREF 

BTAB 

2 

7 

it 

8 

24 

1 10778 

191116 

0 

31 

1 

31 

1 

BA 

BTAB 

2 

9 

it 

9 

26 

110778 

191116 

0 

5 

5 

5 

0 

CON 

1 

0 

it 

10 

27 

110778 

191116 

0 

45 

5 

45 

1 

QJJT 

BTAB 

2 

19 

it 

11 

29 

110778 

191120 

0 

72 

49 

882 

0 

DEF 

E21 

1 

2 

ELD 

12 

61 

110778 

191120 

0 

2 

1 

2 

0 

GD 

E21 

1 

2 

if 

13 

62 

1 10778 

191120 

0 

15 

1 

15 

4 

GTIT 

E21 

1 

2 

it 

14 

63 

1 10778 

191120 

0 

20 

1 

20 

0 

DIR 

E21 

1 

2 

it 

15 

64 

1 10778 

191120 

0 

1 

1 

1 

4 

ELTS 

NAME 

0 

0 

if 

16 

65 

1 10778 

191120 

0 

1 

1 

1 

0 

ELTS 

NNOD 

0 

0 

ft 

17 

66 

1 10778 

191120 

0 

1 

1 

1 

0 

ELTS 

ISCT 

0 

0 

It 

18 

67 

110778 

191120 

0 

15 

1 

15 

0 

NS 

0 

0 

If 

19 

68 

1 10778 

191122 

0 

896 

5 

896 

0 

KMAP 

9 

3 

TOPO 

20 

100 

110778 

191122 

0 

1792 

5 

1792 

0 

AMAP 

@88§ 

9 

3 

tt 

21 

164 

110778 

191126 

0 

560 

4 

140 

4 

E21 

EFIL 

1 

2 

E,  EKS 

22 

184 

1 10778 

191124 

0 

30 

5 

30 

-1 

OEM 

DIAG 

0 

0 

tt 

23 

186 

110778 

191127 

0 

2240 

5 

2240 

1 

K 

SPAR 

36 

0 

K 

24 

266 

110778 

191129 

0 

3584 

5 

3584 

1 

INV 

K 

1 

0 

INV 

25 

394 

110778 

191131 

0 

30 

1 

15 

4 

CASE 

TITL 

1 

1 

AUS/  ALPHA 

26 

396 

110778 

191131 

0 

60 

5 

30 

-1 

A PPL 

FORC 

1 

1 

AUS/SYSVEC 

27 

400 

110778 

191133 

0 

60 

5 

30 

-1 

STAT 

DISP 

1 

1 

SSOL 

28 

404 

1 10778 

191133 

0 

60 

5 

30 

-1 

STAT 

REAC 

1 

1 

tt 

29 

408 

110778 

191135 

0 

208 

4 

208 

-1 

STRS 

E21 

1 

1 

GSF 

30 

416 

110778 

191135 

0 

208 

4 

208 

-1 

STRS 

E21 

1 

2 

tt 

DCU,  the  Data  Complex  Utility  Program,  enables  user's  to  perform  many  data 
management  functions,  e.g.  copying  data  sets  from  one  library  to  another, 
disabling  data  sets,  re-enabling  previously  disabled  data  sets,  etc. 
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TYPICAL  DYNAMIC  ANALYSIS 


Data  exchange  between  all  models  is  accomplished  automatically  via  data  base. 


Analysis  of  frame  model! 


Vibrational  modes 


Joint  force/moment  elgenmatrices 


Dynamic  response,  e.g.  histories 
of  generalized  coordinates  and 
time  derivatives  as  required.  ' — ~ 


Joint  motions,  reactions,  gross 
resultants  (e.g.  overturning 
moments) , or  any  other  state 
quantities. 


Analysis  of  multiple  detailed  shell 
and/or  models  of  components; 


•Joint  flexibility  coefficients 


Stress  transformation  matrices  associ- 
ated with  unit  force/moment  resultants 
applied  as  membrane  loads  at  edges. 


Dynamic  stress  transformation  matrices 


Detailed  dynamic  stress  response  of  the 
entire  model,  or  of  any  number  of  sel- 
ected key  quantities. 


Data  base  scanners, 
graphic  display  generators, j 
custom  report  generators.  V 


A 
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SPAR  Structural  Element  Repertoire 
Name  Description 

E21  General  beam  elements  such  as  channels, 
wide> flanges,  angles,  zees,  tubes,  etc. 

E22  Beams  for  idiich  the  intrinsic 
stifftiess  matrix  is  given. 

E23  Bar  - Axial  stiffness  only. 

E2M  Plane  beam . 

E25  Zero-length  element  used  to 

elastically  connect  geometrically 
coincident  Joints. 

Two-dimensional  (area)  elements: 

E31  Triangular  membrane. 

E32  Triangular  plate . 

E33  Triangular  combined  membrane  and 
bending  eleaent. 

E41  Quadrilateral  membrane. 

E42  Quadrilateral  plate. 

E43  Quadrilateral  combined  membrane  and 
bending  element. 

E44  Quadrilateral  shear  panel. 

Three-dimensional  solids: 

S41  Tetrahedron  (pyramid) . 

S61  Pentahedron  (wedge) . 

S6l  Hexahedron  (brick). 

Compressible  fluid  elements: 

F41  Tetrahedron  (pyramid) . 

F6l  Pentahedron  (wedge) . 

F8l  Hexahedron  (brick). 


Notes: 

- Aeolotropic  constitutive  relations  permitted,  all  area  elements. 

- Laminated  cross  sections  permitted  for  E33,  £43. 

- Membrane/ bending  coupling  permitted  for  E33,  E43. 

- E41,  E42,  E43,  E44  may  be  warped. 

- Aeolotropic  constitutive  relations  permitted  for  3-D  solids. 

- Non-structural  mass  permitted  for  line  and  area  elements. 
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SPAR  Processor  Functions 


Name  Function 

TAB  Translates  user  inputs  into  data  sets  containing 

basic  tables  of  Information  such  as: 

••  Joint  locations. 

- Material  constants. 

- Element  section  properties. 

> Joint  reference  frame  orientations. 

- Constraint  conditions. 

ELD  Produces  data  sets  containing  basic  element 

definitions,  i.e.  connected  Joints,  integers 
pointing  to  applicable  lines  in  tables  of  section 
properties,  material  constants,  etc. 

E Generates  a system  of  daica  sets  called  the  'E-state,' 

consisting  of  individual  element  information  packets 
containing  data  such  as  clement  geometry  (dimensions, 
orientation),  and  literal  section  properties. 

E also  foms  the  system  diagonal  mass  matrix. 

EKS  Computes  element  stiffness  and  stress  influence 

matrices,  and  inserts  them  into  the  'E-state' . 

TOPO  Analyzes  element  interconnection  topology,  and 

produces  data  sets  used  to  guide  other  SPAR  processors 
in  forming  and  factoring  assembled  system  matrices. 

K Forms  system  elastic  stiffness  matrix. 

M Forms  system  consistent  mass  matrix. 

KG  Forms  system  geometric  (pre-stress)  stiffness  matrix. 

FSM  Forms  system  matrices  (dilitational  strain  energy, 

gravitational  energy,  kiretlc  energy)  associated  with 
fluid  elements. 

INV  Factors  system  matrices  in  SPAR's  standard  sparse- 

matrix  format,  e.g.  K,  K-t-Kg,  K-cN, 

AUS  The  Arithmetic  Utility  System,  containing  an  array  of  subprocessors 

in  the  following  categories: 

- Source  data  table  construction  and  editing.  For  example,  the 
following  commands  cause  a 2 x 3 matrix  to  be  created  and  stored 
in  the  data  base  in  a data  set  named  XYZ: 


TABLE ( NI 

=2,  NJ*3): 

J=1:  2.3 

5.7 

Js2:  3.4 

4.2 

J=3:  1.2 

8.0 

6 
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SPAR  Processor  Funccions  (continued) 


Nane 

AUS(cont) 


RMK 

EQNF 

SSOL 

GSF 

PSF 

ES 

EIG 

CEIG 

DR 


Function 

• Matrix  arlthmetis  operations.  For  example,  the  following 
command  causes  a new  data  set,  named  KSUM,  to  be  created  and 
stored  in  the  data  base.  KSUM  is  the  sum  of  two  matrices  residing 
within  data  sets  named  K and  KG.  These  matrices  may  be  in  any  of 
a number  of  forms,  e.g.  total  system  matrices  in  sparse  format: 

KSUM=  SUM(  K,  KG) 

- Special  functions,  including  subprocessors  used 
in  performing  substructure  analysis. 

Translates  arbitrary  M and  K data  into  SPAR  format 
system  matrices. 

Computes  fixed- Joint  forces  associated  with  thermal, 
dislocational , and  pressure  loading.  Computes 
element  generalized  initial  strain  arrays. 

Computes  Joint  motions  and  reactions  due  to  static 
loading. 

Produces  data  sets  containing  element  stresses  and 
internal  loads. 

Produces  tabular  stress  reports  from  data  sets 
generated  by  GSF. 

The  Element  State:  Processor,  which  will  supersede 
GSF  and  PSF.  ES  performs  an  array  of  functions, 
including  stress  scans  and  automated  production 
of  dynamic  stress  transformation  matrices. 

Solves  high-order  eigenproblems  involving  system 
matrices  in  SPAR's  sparse  matrix  format.  Used  to 
solve  both  vibrational  and  buckling  eigenproblems. 

Computes  ocxnplex  modes  and  frequencies  of  damped, 
spinning  structures.  System  matrices  are  in  SPAR's 
standard  sparse  matrix  format,  permitting  analysis 
of  systems  of  very  high  order. 

Computes  linear  transient  modal  response. 


Name 

SIN 

STRP 

SSBT 

SH 

DCU 

VPRT 

PR 

PS 

PLTA 

PLTB 


SPAR  Processor  Functions  (continued) 

Function 

Synthesizes  system  M and  K from  substructure  data 
In  the  form  produced  by  AUS  subprocessors  SSPREPi 
SSH,  and  SSK. 

General  purpose  elgensolver,  full  mass  and  stiffness 
matrices.  Used  primarily  In  analyzing  systems 
synthesized  by  SYN. 

Substructure  back- transformation  processor. 

Computes  Joint  motions  in  individual  substructures 
from  system  state  data  in  the  form  generated  by 
SYN  and  STRP. 

•<» 

The  System  Modification  Processor.  SM  alters  the 
basic  definition  of  thM  structure  to  cause  modes 
and  frequencies  to  apprcach  target  values  defined 
by  the  user.  Typical  applications  include  tuning 
finite  element  models  to  agree  with  dynamic  test 
results,  and  design  of  vibration  attenuators. 

The  Data  Complex  Utility  Program.  DCU  performs 
utility  operations  such  as  printing  data  base 
tables  of  contents,  copying  data  sets  from  file 
to  file,  printing  selected  items  from  data  sets, 
and  transferring  data  t6  or  from  programs  outside 
the  SPAR  system. 

Prints  reports  of  data  in  SPAR's  SYSVEC  (system 
vector)  format,  e.g.  static  displacements, 
reactions,  vibrational  or  buckling  eigenvectors. 

Generates  reports  of  the  results  of  dynamic  response 
analyses.  A variety  of  display  formats  are  provided 
for  the  results  of  transient  and  ranciom  response 
studies. 

Prints  designated  parts  of  SPAR-format  system 
matrices . 

Transforms  user  inpu^'i;  into  data  sets  detailing  the 
composition  of  plots  to  be  produced  by  PLTB. 

Produces  plots  of  deformed  or  undeformed  structure, 
stresses,  etc. 


Ill 
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Name 


K21 

K31 

K41 

K61 

K81 


C21 

C31 

C41 


C32 

C42 

C62 


MT21 


MM2 

HT62 


R21 

R31 

R41 


Name 

TGEO 

SSTA 


nRA 
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<• 


Thennal  Element  Repertoire 


Description 

Conducting  Elements: 

2 node  line  element 

3 node  area  element 

4 node  area  element 

6 node  volume  element 
8 node  volume  element 

Convection  to  a Known  Temperature: 

2 node  line  element 

3 node  area  element 

4 node  area  element 

Fluid-Surface  Convective  Exchange: 

3 node  line  element 

4 node  line  element 
6 node  area  element 

Mass-Transport : 

2 node  line  element 

Integrated  Mass-Transport,  Convective  Exchange: 
4 node  line  element 
6 node  area  element 

Radiating  Elements: 

2 node  line  element 

3 node  area  element 

4 node  area  element 


Thermal  Processor  Functions 


Function 

Creation  of  primary  thermal  element  network  data  structures. 
Computes  steady  state  solutions,  both  linear  and  nonlinear. 
Computes  transient  solutions,  both  linear  and  nonlinear. 
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DSCS  in  STACKED  LAUNCH  CONFIGURATIONS 
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•XltAV  SEPARATIONS 
NO  RAILS  RCQUIREO 


MODAL  TEST  OBJECTIVE  AND  CRITERIA 
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WITH  VERIFICATION  LO^.D  CYCLE 


OVERALL  APPROACH 
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VALIDATE  SPR  TECHNIQUE 
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STRUCTURAL  TRANSFER  FUNCTION  DATA 
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APPROXIMATELY  6% 
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TWO  EXCITATION  SYSTEMS  WITH 
COMMON  DATA  ACQUISITION 


APPROXIMATELY  40 
MODES  BETWEEN  10 
AND  50  HERTZ 


MULT(-SHAKER  SINE  TECHNIQUE 


SINGLE  SHAKER  RESPONSE 


C9  lii 
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FREQUENCY  HZ 


MODE  C NULLED 
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ORTHOGONALITY  OF  REPEATED  MODES 


4m,  .1 


U 


ca  uj 


130 


REPEATED  MODES  SHOW 
IMPROVED  ORTHOGONALITY 


SUMMARY  OF  PRESENT  RESULTS 
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OBJECTIVES 
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APPIOACH 
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WHITE’S  FORMULATION  (contd) 
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• UNDERLINE  ■ TEST 


METHOD 
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MEMBER  FACTOR  CONVERGENCE 
OF  PERTURBED  ANALYTICAL  MODEL 
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a 890  1.059  1.020  0.989 


FREQUENCY  CONVERGENCE 
F PERTURBED  ANALYTICAL  MODEL 
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CONSIDERED  AS  THE  TEST  DATA  TO  WHICH  THE  PERTURBED  MODEL  IS 
TO  CONVERGE. 


X* 


FREQUENCY  COMPARISON  OF  TEST  RESULTS 
AND  ORIGINAL  MODEL 


MODE 

NO. 

FREQUENCY, 

Hz 

ORIGINAL  MODEL 
FREQUENCY, 

Hz 

DIFFERENCE, 

Hz 

DIFFERENCE, 

% 

1 

12.95 

12.37 

0.58 

4.5 

2 

17.66 

15.85 

1.81 

10.2 

3 

20.80 

19.37 

1.43 

6.9 

4 

22.97 

26.49 

-3.52 

-15.3 

5 

28.31 

27.90 

0.43 

1.5 

6 

32.76 

34.19 

-1.43 

4.4 

7 

42.80 

42.06 

0.74 

1.7 

8 

50.67 

47.69 

2.98 

5.9 

9 

50.40 

52,39 

-1.99 

3,9 

10 

52.51 

11 

60.47 

12 

65.38 

60.96 

4.42 

6.8 

RSS  OF  DIFFERENCE 

7.29  Hz 

12 
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DIFFICULT  IF  NOT  IMPOSSIBLE  FOR  COMPLEX 
STRUCTURE. 


IF  IttSS  AND  STIFFNESS  MATRICES  ARE  AS  FOLLOWS: 
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EVALUATION  OF  JACOBIAN 
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OPERATIONAL  PROCEDURE 
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FIGURE  1.  3 DEGREES-OF-FREEDOM  SAMPLE  PROBLEM 
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COMPARISONS  OF  FREQUENCY  AND  MODE  SHAPE  FOR  HASSELMAN'S  .METHOD 
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WHITE'S  METHOD 
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Application  of  Perturbation 
Methods  to  improve 
Anaiytical  Model  Correlation 

with  Test  Data 

J.  A.  Garba  and  B.  K.  Wada 

Structures  and  Materials  Section, 
Jet  Propulsion  Lab. 


APPRECIABLE  EFFORT  Is  being  expended  on  the 
•reflearch  and  development  of  computer  programs 
for  the  analytical  simulation  of  structures.  A 
parallel  effort  consists  of  the  improvement  of  test 
methods  to  measure  the  dynamic  characteristics  of 
complex  structures.  Both  of  these  have  provided 
the  engineer  with  valuable  tools.  The  need  has 
existed  in  the  past  and  will  exist  in  rhe  future  for 
capabilities  to  systematically  update  a mathematical 
model  to  more  accurately  represent  the  test  data. 
Currently,  most  organizations  use  a ’’trial  and 
error”  approach  to  improve  mathematical  models. 
The  demand  in  the  future  to  correlate  mathematical 
models  with  test  daca  will  be  greater  as  the  models 
become  more  ccmple?:.  More  reliance  will  be 
placed  on  test  correlated  analytical  models  to 
verify  modifications,  and  some  future  requirements 
are  for  large  structural  systems  which  cannot  be 
fully  ground  tested.  For  such  systems  confidence 
must  be  established  by  analyses  with  test  verifica- 
tion of  analytical  models  of  subsystems  or 
components. 

The  paper  is  limited  to  analysis/test  correlation 
rather  than  ’’system  identification,  ” which  aims  to 
create  a mathematical  model  that  will  reproduce  the 


test  results.  Analysis/test  correlation  is  considered 
important  because: 

(1)  The  model  is  usually  used  to  obtain  struc- 
tural design  loads.  Thus,  it  contains  information 
vital  to  the  structure,  such  as  modal  force 
coefficients. 

(2)  The  model  can  be  made  physically  meaning- 
ful to  simulate  the  test  conditions.  Thus  modifica- 
tions to  the  model  to  simulate  flight  coriditions  are 
possible.  The  modifications  may  be  result  of 
anticipated  design  changes  or  to  eliminate  some 
nnuosirable  ground  test  condition,  such  as  a 
gravitational  force. 

(3)  The  use  of  analysis/test  correlated  subsys- 
tems in  obtaining  system  models  using  modal  syn- 
thesis techniques  result  in  better  simulations. 

The  overall  plan  for  the  development  of  a method 
for  the  systematic  correlation  of  anah^ical/test 
models  is  to  evaluate  the  applicability  of  the  pub- 
lished methodologies  to  ’’real  structures”,  to 
apply  new  approaches  as  they  arc  evolved  during 
the  effort,  to  select  the  most  promising  ap.  roach 
based  upon  the  experiences,  aiid  ihun  generate  a 
computer  program  which  is  compatilile  to  user 
oriented  programs  such  as  NASTRAN.  The  above 


ABSTRACT 


There  are  current  and  future  requirements  to 
develop  a systematic  method  to  update  a mathe- 
matical model  of  a structure  to  more  closely  match 
the  test  data.  The  effort  is  cost  effective  since  the 
number  of  reanalyses  of  a lai^je  structure  will  be 
reduced.  Additionally,  the  mathematical  model 
will  more  closely  represent  its  test  data. 
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This  activity  describes  the  positive  and  negative 
experiences  in  using  a method  putilishod  oy  C.  W. 
White  and  an  extension  of  the  method  The  results 
are  based  iq>on  our  understanding  of  the  method  as 
published. 


/•/ 


steps  are  considered  necessary  to  use  technology 
developed  to  date,  apply  the  methodology  to  a 
realistic  problem  (often  the  methodology  has  been 
shown  to  be  successful  only  on  simple  problems 
or  a mathematically  constructed  set  of  "test  data"), 
and  to  obtain  user  experience  prior  to  commitment 
of  funds  to  develop  a computer  program.  Many 
problems  with  various  proposed  methods  are  only 
fully  realized  after  an  attempt  at  its  utilization. 

This  paper  is  an  application  of  an  approach 
published  by  C.  White  to  the  Viking  Propulsion 
Subsystem  (VPS)  for  which  anal3rtical  models  and 
modal  test  data  are  available.  The  VPS  allows 
variation  of  the  mass  and  stiffness,  and  is  suffi- 
ciently complex.  In  addition,  another  method  that 
evolved  during  the  work  will  be  presented.  The 
attempt  is  to  describe  the  experiences  that  were 
both  successful  and  unsuccessful.  The  authors 
have  applied  C.  White's  method  to  the  selected 
problem  to  the  best  of  our  ability.  The  theoretical 
development  of  C.  White  is  repeated  only  to  docu- 
ment our  interpretation  of  his  effort.  In  addition 
to  using  the  actual  VPS  test  data  for  correlation,  a 
set  of  mathematically  generated  test  data  was  also 
employed  to  evaluate  the  methodology.  The  result- 
ing perturbed  analytical  model  is  representative  of 
actual  flight  hardware  and  hence  more  lealistic 
than  the  simple  models  usually  used  in  such  studies. 
Other  methods  such  as  those  proposed  by  members 
of  Wiggins  and  Kaman  Corp.  are  in  the  review 
process  and  will  be  applied  to  the  VPS  to  understand 
their  merits. 

The  "best"  methodology  may  not  exist,  but  may 
be  dependent  on  the  accuracy  desired,  or  the 
dynamic  characteristics  of  significance  for  the 
model  usage.  The  paper  assumes  that  a correlation 
of  most  dynamic  related  characteristics  is  important. 
These  characteristics  include  eigenvalues,  eigen- 
vectors, kinetic  energy  distribution,  strain  energy 
distribution,  cross-orthogonality,  effective  weight, 
and  modal  force  coefficients. 

LINEAR  PERTURBATION  METHOD 


a transformation  from  discrete  to  normal  mode 
coordinates  can  be  introduced  as 


(2) 


adiere 

{{}=  normal  mode  or  generalized  coordinate 
vector 


= eigenvector  matrix  obtained  as  the 
-*  solution  of  the  following  equation 


[»o]  <»> 


where 


rWg  I = diagonal  matrix  of  the  system 
* circular  frequencies  squared 
Substituting  Eq.  (3)  into  Eq.  (1)  and  premultiplying 
by  [p  o]  ^ leads  to 


{■«■>»  m (4). 

by  virtue  of  the  normalization  of  * 

A linear  perturbation  in  the  nominal  stiffness 
and  mass  matrices  can  now  be  introduced. 


Q 

[“] ' [«o]  H ‘ [“»]  ^ H«] 

q=l  ^ 


The  perturbation  method  follows  the  development 
by  C.  W,  White,  Refs.  (1),  ♦ (2).  and  (3).  Only 
the  pertinent  equations  of  the  derivation  are  pre- 
sented here. 

Starling  witli  the  homogeneous  equations  of 
motion  for  the  structure  with  negligible  damping 
as 

[M  J {H}  + [K  J {h}  = {0}  (1) 

where 

{h}  = global  system  absolute  discrete 
displacement  vector 

'Moj  = the  initial  model  mass  matrix 
j - the  initial  model  stiffness  matrix 


where  [k]  and  [m]  are  the  perturbed  system  stiff- 
ness and  mass  matrices,  respectively.  The  factors 
6p  and  6q  denote  the  linear  variations  of  the  affected 

element  stiffness  and  mass  matrices.  A total  of  P 
stiffness  elements  and  Q mass  elements  are 
assumed  to  have  been  perturbed. 

Using  Eqs.  (5)  and  (6),  a perturbed  eigenvalue 
problem  similar  to  Eq.  (2)  can  be  formulated  as 


([o-w'Nw)[.]r'-j- 


‘J  f *.A  ' ^ 

K V/ 


’•‘Numbers  In  parentheses  designate  References 
at  end  of  paper. 
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// 


' *■  / 


vdiere 


[^2  m eigenvector  matrix  of  the  perturbed 
problem 


diagonal  matrix  of  the  perturbed 
system  circular  frequencies 
squared 

The  transformation  from  discrete  to  normal  mode 
coordinates'implied  by  Eq.  (7)  is 


(8) 


where 

[h  1 = a set  of  experimentally  measured 

modes 

Examining  Eq.  (9),  the  authors  in  Ref.  (1)  note  that 
the  right-hand  side  represents  the  potential  and 
kinetic  energy  of  the  perturbed  stiffness  and  mass 
elements.  Eq.  (9)  can  be  rewritten  as 


where 

{y}  ^ the  normal  modes  coordinate  vector  of 
the  perturbed  system 
Equation  (7)  can  be  rewritten  as 


q-1  ^ 

In  terms  of  model  correlation,  the  terms  in 
Eq.  (9)  can  be  interpreted  as  follows: 

is  related  to  the  frequencies  of  the 
analytical  model  to  be  correlated 

is  related  to  the  measured  frequencies 

is  the  set  of  eigenvectors  of  the 
analytical  model 

is  the  set  of  eigenvectors  relating 
the  analytical  to  the  experimental 
modes 

The  matrix  [ ip  ] has  the  property 

[it]  - [I]  (10) 


for  an  unperturbed  problem,  meaning  that  there 
exists  perfect  correlation  between  test  and  analysis. 
For  go^,  but  not  perfect  correlation,  the  matrix 
( fp  )is  strongly  diagonal,  meaning  that  the  size  of 
the  diagonal  terms  are  much  larger  than  the  off- 
diagonal  terms. 

In  terms  of  measured  modes,  the  matrix  ( 0 1 
can  be  obtained  from 


(11) 


where  every  term  in  matrix  [r]  contains  the  P + Q 
unknown 

In  theory  Eq,  (12)  leads  to  a set  of  equations 
that  can  be  solved  for  the  unknown  factors  6 . In 
practice  there  are  several  difficulties,  namely: 

(1)  The  solution  of  Eq.  (11)  for[^]  requires 
the  inversion  of  a typically  very  large  nonsymmetric 

matrix  [^o]. 

(2)  The  measured  modes,  [h^xp]  usually  contain 
fewer  measurements  than  the  original  vector  {h}  of 
the  analysis  requires.  This  leads  to  an  incompati- 
bility in  the  solution  for[^]  in  Eq.  (11). 

(3)  Eq.  (12)  can  be  expanded  into  1/2  (n^  + n) 
equations  in  P + Q unknowns,  where  n is  the  number 
of  normal  modes  being  considered.  If  the  number 
of  unknowns  equals  the  number  of  equations,  the 
solution  can  be  obtained  by  matrix  inversion. 
Otherwise,  a least  squares  fit  or  a linear  program- 
ming solution  must  be  used  depending  on  if  the  sys- 
tem is  underdetennined  or  overdetermined. 

Although  Eq.  (12)  can  lead  to  an  exact  solution, 
practical  considerations  require  Ihe  search  for  a 
simplified  approach,  possibly  leading  to  an  iterative 
solution  to  the  problem  of  determining  the  required 
factors  6 . The  authors,  in  Ref.  (1),  develop  such 
a simplification  as  follows 

If  in  Eq.  (9)  only  the  i-th  column  is  considered, 
the  equation  simplifies  to 

where  columns  of  [e]^  are 

or 
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The  solution  of  Eq.  (13)  is  somewhat  simplified  from 
Eq.  (12)  but  it  is  still  restricted  by  the  solution  to 
Eq<  (11).  Further  simplifications  can  be  made  by 
assuming  that  the  perturbed  modes  are  close  to  the 
original  modes,  then,  from  £q.  (10), 

[*] " 

If,  furthermore,  in  the  formulation  of  the  solution 
only  the  diagonal  terms  are  considered,  the  set 
of  equations  simplify  to 

p-"S|=[E„]W  (15) 

i^lch  leads  to 

{«)=[eJ  jn*-»51  (16) 

for  the  case  where  the  number  of  the  unknowns 
equals  the  inodes  under  consideration,  or 

for  a least  squares  fit  solution  where  the  number  of 
equations  exceed  the  unknown  6*8. 

Equations  (15)  and  (16)  form  the  basis  of  the 
model  correlation  used  herein.  It  should  be  noted 
that  the  simplifications  used  in  the  above  derivation 
have  eliminated  the  use  of  eigenvectors  and  rely 
solely  on  frequency  correlation.  While  frequencies 
can  be  measured  accurately,  the  method  proposed 
by  Eqs.  (15)  and  (16)  is  essentially  nothing  more  than 
an  eneigy  balance.  The  drawbacks  of  this  will  be 
discussed  later  in  the  application  of  this  method. 
Since  energy  is  a scalar  and  the  eigenvectors  are 
vectors,  it  is  obvious  that  by  using  only  frequency, 
rather  than  frequency  and  mode  shape,  the  modal 
test  data  is  not  used  to  its  full  potential. 

Before  proceeding  to  apply  the  above  method, 
let  us  consider  an  alternate  formulation  of  the  per- 
turbation method  as  formulated  by  J,  C,  Chen, 

Kefs.  (4)  and  (5).  Consider  a svstem  of  equations 
for  the  structural  system  analogous  to  Eq.  (1): 

[M]  {x}  +[K]  {x}  = {0}  (18) 

and  Introduce  the  linear  perbirbations: 

M ■■=  [Mq]  [mJ 

[K]  + ej-Kj]  (19) 


t:j 


= N J 

|q}  (20) 

where  }q|  is  the  normal  mode  generalized  coordi- 
nate. In  the  context  of  correlation,  CMq^,  CKq] 
andC^o3  refer  to  the  analytical  model,  whereas 
Chfl,  IK],  and  [0]  refer  to  the  test  values.  Hence 
it  is  desirable  to  find  eMi,  and  eK2,  by  perturbing 
a set  of  masses  and  stiffnesses  and  matching  fre- 
quencies and  mode  shape.  Substituting  Eqs.  (19) 
ana  (20)  into  Eq.  (18),  we  obtain  two  sets  of 
equations. 

[0]  [M][0]  jq|  + W^[K][0]  jq}  = {o| 
and 

[00  ^ ^ [^0  ^^1]  [^0  ^^1]  1^'  f 

f*  "1 X 

Expanding  Eq.  (21)  and  neglecting  the  second-order  ' 
terms  and  noting  that 


•'.V  ?Ot.  «» 

hrwH 

1^1  =[0j 


[0]  [M]  [0]  = [I] 

W'NN-w 

[♦f  [K]  [,].[■ 


(22) 


where  as  before  the  SI  and  u>q  related  to  the  test  and 
analytical  frequencies,  respectively,  leads  to 

[s]  l<i  I *[' {,).{»} 

and 


(23) 
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for  the  two  equations  to  be  equlvalentt  the  coeffi- 
cients are  equated: 

[•.rhiw-([-.rww| 

-(W’HN 

•WWW) 


and 


HW-r-’-J-r-SJ  L 

■^!JW 

•[••rr-ij)  ) 

Eq.  (28)  can  now  be  solved  for  either  {fig}  or  {6  } 
by  evaluating  Eq.  (26)  with  ** 


[‘*i]=[«]-W  <»»> 


Blnce  the  eigenvectors  are  linearly  independent  of 
each  other,  any  vector  can  be  expressed  as  a 
linear  combination  of  the  complete  set,  thus 


where  the  matrix  (a]  contains  the  linear  factors. 
Using  Eq.  (25),  we  can  establish  the  following 
relationships: 

[*,]"[“«]['».]■[“]  ) 

> (26) 

Substituting  Eqs.  (26)  and  (24)  and  noting  that  in 
the  notations  used  earlier 


W’(-.1W-E‘W’W>] 

P»1 


WK} 


(27) 


[*o]  ['“i][*o]-E‘a[*o]  ["'IjM 

q»l  '* 

-WK} 


Since  in  the  set  of  Eqs.  (28)  the  effect  of  pertur- 
bations in  mass  and  stiffness  matrices  are  imcoupled, 
one  has  a choice  of  varying  one  or  the  other  at  a 
time  in  an  iteration  cycle,  or  one  can  separate  the 
modes  effected  by  mass  and  stiffness  and  use  both 
of  the  equations  in  the  same  iteration  cycle. 

Each  of  the  two  Eqs.  (28)  lead  to  1/2  (n^  + n) 
equations  in  either  P or  Q unknowns,  similar  to  the 
solution  of  Eq.  (12). 

On  the  surface  it  appears  that  the  alternate 
formulation  of  Eqs.  (28)  has  the  advantage  over 
Eq.  (12)  by  making  use  of  the  measure.^  mode 
shape  information  (cPi) . The  practical  problems 
of  applying  Eq.  (28)  in  model  correlation  will  be 
discussed  later. 

CORRELATION  MODEL 

The  Viking  Propulsion  Subsystem,  which  was  an 
important  structural  subsystem  of  the  Viking 
Orbiter,*  was  selected  for  correlating  an  analytical 
model  with  test  results.  It  acco"  its  for  70%  of  the 
Viking  Orbit  er  weight,  and  the  . ght  configuration 
contains  large  quantities  of  fluids  with  ullage. 

The  Viking  Orbiter  was  designed  by  loads 
analysis.  This  process  relies  cn  a rsprese.itative 
mathematical  model  to  obtain  flight  loads  to  achieve 
a reliable  design.  A good  correlation  between  the 
test  results  and  analytical  predictions  is  of  para- 
mount Importance.  Data  obtained  during  the 
powered  phase  of  the  two  Viking  flights  have  shown 
very  good  correlation  between  predicted  loads 
and  flight  measurements,  Ref.  (6). 

Other  considerations  that  make  this  subsystem 
appropriate  for  the  application  of  the  perturbation 
technique  are: 

(1)  The  structure  was  thoroughly  tested  and 
analyzed;  the  results  are  documented  in  Kefs.  ^) 

♦The  Viking  Orbiter  <s  part  of  the  Viking  Space- 
craft, which  was  flown  successfully  in  1975. 
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Pig.  1 - Schematic  of  the  Viking  propulsion 
subsystem  finite  element  model 

PERTURBED  ANALYTICAL  MODEL 

Before  the  perturbation  method  was  applied  to 
the  analytical  model  to  obtain  an  imprcA  ed  test/' 
analysis  correlation  with  test  data,  the  existing 
analytical  model  was  perturbed.  This  perturbed 
model  was  then  treated  as  the  analytical  model, 
and  the  existing  analytical  model  was  considered 
as  the  test  data.  The  method  developed  earlier 
was  applied  to  the  assumed  analytical  model  to  try 
to  reconstruct  the  original  answers. 

The  purpose  of  this  was  two-fold: 


and  (8).  The  eigenvalues,  eigenvectors,  and 
modal  damping  were  experimentally  determined 
by  a modal  test. 

(2)  An  acceptable  correlation  between  test  and 
analysis  was  obtained  previously.  This  correlation 
was  achieved  in  a project  environment  with  severe 
schedule  constraints;  it  was  obtained  by  inspection 
rather  than  a systematic  fashion.  The  remaining 
differences  appear  small  enough  that  perturbation 
theory  is  applicable  for  further  correlation  attempts. 

(3)  The  structure  is  of  moderate  complexity. 

It  contains  thin~wall  pressure  vessels  supported  by 
beams.  The  local  tank  tab  areas  were  modeled 
using  the  Guyan  stiffness  matrix  reduction  to  gen-* 
erate  superelements  (Ref.  (9)).  It  represents  a 
sufficiently  complex  three-dimensional  structure; 

a good  test  case  for  the  various  proposed  methods 
for  col  I elation.  Usually  the  examples  consist  of 
either  continuous  beams  representing  launch 
vehicles,  or  three-dimensional  structures  using 
axial  members  only, 

(4)  The  design  of  the  major  load  carrying 
memters  was  effected  mainly  by  the  lowest  six 
normal  modes.  Some  higher  modes  were  of  impor- 
tance only  in  the  design  of  local  structure.  The 
correlation  thus  could  concentrate  mainly  on  the 
lowest  six  modes. 

(5)  The  variations  of  large  masses  were 
Important  unknown  parameters.  The  effective 
masses  of  the  fluids  were  important  considerations. 
Extensive  testing  (Ref.  (10))  was  performed  for 
their  determination. 

Several  configurations  of  the  propulsion  subsys- 
tem were  used  in  the  modal  test,  these  v ere  (1)  the 
flight  configuration  containing  appreciable  ullage 
volume  in  both  the  fuel  and  oxidizer  tank,  ^^)  both 
tanks  completely  filled  with  referee  fluid  wiA 
no  ulkige,  and  (3)  both  tanks  empty.  These  con- 
figurations were  chosen  to  verify  the  model 
inaependently  of  the  effect  of  the  fluid  on  the 
normal  mooes.  The  model  discussed  herein  repre- 
sents tike  second  case,  both  tanks  full. 

The  atiai>i:ical  model  contains  677  static  (stiff- 
ness) degrees  of  freedom  and  84  dynamic  (mass) 
d^rees  of  freedom.  There  are  a total  of  5 types 
of  tank  tab  supereleinent  stiffness  models, 

117  plate  elements  and  184  beam  elements. 

Figure  1 shows  the  schematic  of  the  aiuilytical 
model  including  the  coordinate  syi?teni.  Figure  2 
is  a photograph  of  the  test  hardware  used  in  the 
jrnt»di'.l  test. 

The  frequency  correlation  between  the  teet 
values  and  the  l)C8t  post  test  model,  hereafter 
called  the  original  model,  is  listed  in  Table  1. 

Note  that  there  exists  a one-to-one  correspondence 
In  all  the  modes  except  modes  10  and  11.  The 
analysis  shows  that  these  are  local  thrust  plate 
and  engine  modes  that  were  not  excited  during 
the  modal  test.  Since  they  are  of  no  consequence 
In  tlie  design  ut  primary  structural  members,  these 
two  modes  will  be  dropped  from  any  further 
discussion.  The  root  sum  squared  (RSS)  value  of 
the  frequency  difference  in  Table  1 is  7. 29  Hz. 


(1)  Check  the  method.ology  using  tije  same 
analytical  model  which  is  to  be  used  for  the  test/ 
analysis  correbtion. 

(2)  Establish  memoer/mode  grouping  for 
model  convergence. 

The  perturbed  model  was  generated  by  varying 
a set  of  member  stiffnesses  and  retaining  the 
masses  unchanged.  The  members  were  selected 
based  on  the  strain  energy  contribution  to  the 
lowest  12  normal  modes.  In  the  process  of 
choosing  the  members  to  be  perturbed,  the  strain 
energy  contribution  for  the  Iteam  elemeutb  wau 
separated  by  stiffness  component.  Thus  Uie  axial 
contribution  was  separated  from  Ikendlng  mui  torsion. 
For  the  tank  tabs  modeled  by  superelements,  a 
single  perturbation  constant  was  assigned  for  each 
of  the  groups  varied. 
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Fig,  2 - Propulsion  subsyst.ni  modal  test  setup 
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A total  of  seven  groups  of  members  were 
selected  to  be  perturbed.  The  member  groups 
and  the  amount  of  original  perturbation  is  listed 
In  Table  2.  This  perturbation  is  similar  to  that 
used  by  C.  W.  White  in  Ref.  (1),  except  tliat  the 
perturbation  in  the  referenced  article  was  applied 
to  axial  members  only  on  a relatively  simple  dynamic 
model.  The  method  outlined  in  Refs.  (1)  to  (3)  was 
now  used  to  systematically  improve  the  model  by 
applying  Eqs.  (16)  and  (17). 

After  carefully  matching  the  strain  energy  and 
the  affected  normal  modes,  the  model  was  indeed 
forced  to  converge  in  four  iterations.  The  fre- 
quency convergence  is  shown  in  Table  3.  The 
corresponding  factors  for  the  perturbed  members 
are  listed  in  Table  4.  The  factors  are  defined  as 
the  cumulative  multipliers  to  be  applied  to  the 
properties  of  the  unperturbed  ciiginal  model. 

Thus  a factor  of  1. 000  means  that  the  group  has 
the  same  properties  as  the  original  model. 

Several  important  observations  made  in  this 
exercise  of  applying  the  method  to  a perturbed 
analytical  model  are: 

(1)  Conveigence  is  not  guaranteed.  As  a matter 
of  fact,  the  method  diverged  twice  to  factors  that 
had  no  physical  significance.  Convergence  was 


Table  1 • Frequency  Comparlacn  of  Teat  Reeulu 
and  Orlflnal  Model 


Oriftnal  Mode! 


Mode 

Freqijency, 

Frequency. 

Difference, 

Difference, 

No. 

Hx 

Hx 

Hx 

% 

Mode  DeecripUon 

1 

12. 93 

12.37 

0.  38 

4.3 

Oxldlxer  and  fuel 
In  Y,  In  pbaee 

3 

17.66 

IS.  83 

1.81 

10.2 

Oxldlxer  tank  in  Z 

3 

20. 80 

19.37 

1.43 

6.8 

Oxldlxer  and  fuel 
In  Z 

4 

2Z,  97 

26.49 

-3.32 

-15.3 

Oxldlxer  and  fuel 
In  Y out-of-phaee 

S 

28.33 

27.90 

0.43 

1.3 

Fuel  tank  in  Z 

6 

32.76 

34.19 

-1.43 

4.4 

Prcisurant  tank 
In  6 , 

y 

7 

42.80 

42.06 

0.74 

1.7 

Local  preaaurant 
control  aaacmbly 
In  X 

8 

30.67 

47.69 

1.98 

3 9 

Preaaurant  tank 
In  Y 

8 

30.40 

82.39 

-1.88 

3.9 

Local  preaaurant 
control  asee  nbly 
In  Z 

10 

32.81 

txKal  thruat  plate 
and  Ttoxxle 

11 

60.47 

loLui  thruat  plate 
and  noixle 

11 

83.38 

60. 9G 

4 41 

6.9 

Local  preaaurant 

control  eeeembly 


nSSeldUffemnc^  7 If  Hi 


f 


TfthU  S • Ptiturbed  An*tyUeal  Modt)  M«nb*r  Orsuflm 


Tfttoli*  3 - Frequency  Conveifence  o(  Perturbed  Analytical  Model 


Oroup  Dcecription 

A Side  bipods.  4 members  supportlnf  the  taidts 

In  a vertical  plane 

B Top  bipods,  4 members  supporting  the  tanks 

In  a boiisontal  plane 

C 3*hole  tank  tabs,  4 superelemcnts  that  (crm 

the  link  between  the  side  bipods  and  the 
rigid  tank 

D Top  Siamese  tank  tabs,  2 superelementa 

oonnecUng  the  two  unks  at  the  top 

E Bottom  Slameac  tank  taba,  2 superelementa 

connecting  the  two  planes  at  the  bottom 

P Preaaurant  lank  round  members,  3 membera 

supporting  the  presiurant  Unk 

O Pressurant  control  asaembly  (PCA)  support 

•tniU,  4 members  supporting  the  PCA 


Table  4 - Member  Factor  Convergence  of 
Perturbed  Analytical  Modal 


♦40% 

-30% 

♦10% 

-w% 

♦40% 

♦30% 

♦40% 


Mode 

JiSi, 


Original 

Model 

Frequency, 

Hi 


Perturbed 

Model  

Frequency,  lat, 
Ht  H* 


Iteration 


2nd. 

H*. 


3rd, 

Hi. 


ill. 


Sth, 

Jl3- 


Member 


Original 


_DT 


•poo^ 


Oroup 

Perturbation 

let 

tod 

3rd 

4th 

A 

1.4 

0.863 

0.994 

1.015 

1.002 

0.997 

B 

0,7 

0.894 

1.019 

1.020 

1.000 

0.H9 

C 

1.3 

0.972 

0.972 

0.946 

0.993 

1.005 

D 

0.1 

0.902 

1.023 

0. 997 

1.000 

1,000 

E 

1.4 

1.774 

1.298 

0.923 

0.996 

0.999 

r 

1.3 

0.821 

0.969 

0.991 

1.030 

1.002 

0 

1.4 

0. 190 

1.059 

1.026 

0.999 

0.979 

attained  by  not  allowing  the  member  factors  to  vary 
more  than  an  arbitrary  factor  in  each  iteration. 

Thus,  for  example,  if  in  the  second  iteration 
member  group  C was  allowed  to  vary  by  6S%,  the 
system  soon  diverged  to  a point  where  some  other 
member  group  was  removed  totally,  resulting  in 
an  unstable  structure.  This  situation  was  corrected 
by  regrouping  the  members  and  modes,  even  though 
the  original  group  seemed  perfectly  logical.  This 
arises  because  Eq.  (15)  is  indeed  nothing  more 
than  an  energy  balance.  Thus  unless  the  member 
groups  are  matched  properly  with  the  proper  modes, 
physically  unrealistic  results  are  obtained.  Mathe- 
matically a solution  of  Eq.  (15)  can  be  found  that 
adds  a multiple  of  one  set  of  members  and  removes 
another  set  altogether.  The  energy  is  balanced,  but 
the  result  is  physically  not  reasonable.  Unfortunately, 
the  matching  of  modes  and  member  groups  is  not 
that  clear-cut,  and  logical  matching  does  not 
assure  physically  reasonable  results  without 
some  prior  knowledge. 

(2)  Simple  models  such  as  the  one  used  in 
Ref.  1 apparently  can  be  made  to  converge  not 
only  uniformly,  but  also  monotonically.  Even  with 
a most  careful  matching  of  modes  and  member 
groups  the  perturbed  model  did  not  converge 
uniformly. 

(3)  The  alternate  formulation  of  Eq.  (18)  failed 
to  produce  physically  reasonable  results  in  the 
first  iteration,  even  though  the  matching  of  modes 
and  member  groups  was  the  same  as  that  used  in 
Eq.  (17).  This  method  was  not  further  pursued. 

The  speculation  is  that  the  perturbations  were  too 
large  for  the  alternate  formulation  leading  to 
unreasonable  perturbations  in  the  n)ode  shape 

that  in  turn  affected  the  solution  of  Eq.  (28). 


1 

12.37 

13.49 

11.87 

12.34 

13.36 

12.37 

12.37 

1 

15. 95 

16.34 

15. 07 

15.63 

15.85 

15. 86 

15.65 

3 

19.37 

20.29 

19.51 

19.94 

19  37 

19  36 

19.37 

4 

29. 49 

26.25 

26.49 

27.15 

25. 24 

26.49 

26. 49 

9 

27.90 

29.30 

27.50 

29  13 

27.95 

27.91 

27.90 

9 

94. 19 

34.29 

34.09 

34.  :6 

34.17 

34.20 

34.19 

7 

42.06 

42.64 

41.66 

42.09 

42.09 

42.08 

42.06 

9 

47.99 

49.43 

45. 10 

47.49 

47.62 

47.93 

47.69 

9 

52.99 

32.61 

51.56 

32.49 

32.48 

52.41 

32.39 

10 

52.51 

54.21 

52.41 

52.65 

52. 65 

52.50 

32.51 

11 

60.47 

61.94 

88.65 

60.34 

60. 30 

60.  52 

60.47 

12 

90. 99 

64.23 

59.49 

61.21 

61.26 

61.05 

60.96 

♦Coneldered  the  test  data  to  whldi  the  perturbed  model  la  to  convene. 

« Comparison  of  Teat  Frequencies  to  the 
Original  and  the  Improved  Models 

Teat 

Original  Model 

Model  A 

Model  B 

Mode 

Frequency, 

Frequency, 

Frequency. 

Frequency, 

No. 

Hz 

Hz 

Hz 

Hz  * 

1 

12.95 

12.37 

13  50 

13.42 

(0  59)* 

(-0,  55) 

(-0.47) 

2 

17.60 

15. 65 

17,59 

17.28 

(1.91) 

(0.07) 

(0.38) 

3 

20.  90 

19.97 

20. 67 

20.98 

(1.43) 

(0. 13) 

(0.18) 

4 

22.97 

26. 49 

22.94 

22. 97 

(-3.52) 

(0.03) 

(0.00) 

5 

29.33 

27.90 

28. 35 

28.34 

(0.43) 

(-0.02) 

(-0.01) 

9 

32.76 

34. 19 

32.74 

32.76 

(-1.43) 

(0. 02) 

(0.00) 

t 

42.90 

42.06 

42.85 

42.85 

(0.  74) 

(-0.05) 

(-0.03) 

9 

50.67 

47.69 

50.32 

50.48 

(2.98) 

(0.08) 

(0. 19) 

9 

50.40 

52.39 

52.23 

52.40 

(-1.99) 

(-1.56) 

(-2.00) 

12 

65.38 

60. 96 

62.38 

62.43 

(4.42) 

(3. 10) 

(2.95) 

RSS  of  Difference: 

All  10  modta 

7.29 

3.02 

3.62 

Flrat  6 modea 

4.50 

0 57 

0.63 

*Th«  values  In  parentheaea  refer  to  the  difference: 

teat  frequency  less  model 

frequency. 

ANALYSIS/TEST  CORRELATION 

The  previous  section  has  established  the  feasi- 
bility of  the  perturbation  method  and  has  also 
revealed  some  of  the  shortcomings  of  the  method 
when  applied  to  a real-life  structure.  The  method 
will  now  be  used  to  try  to  improve  the  corr-^lation 
between  the  test  results  and  the  analytical  model 
of  the  Viking  Propulsion  Subsystem. 

The  approach  will  be  that  of  applying  Eqs.  (16) 
and  (17),  or  (28)  in  an  iterative  scheme.  In 
proceeding  in  this  manner  it  is  fully  realized  that 
the  convergence  criterion  is  based  on  frequency 
alone  and  that  the  method  consists  of  an  energy 
balance.  Once  a reasonable  frequency  match  is 
established  the  mode  shape  is  checked  to  see  if 
an  improvement  in  the  analytical  model  has  been 
achieved. 

The  first  step  in  applying  the  perturbation 
technique  of  Ref.  1 is  to  identify  the  major  modei 
energy  contributors,  both  potential  and  kinetic. 
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Then  the  major  contributing  elements  are  used  with 
the  affected  modes  to  solve  Eq.  (IG)  or  (17). 

It  was  found  that  the  lowest  5 modes  all  had 
appreciable  strain  energy  contribution  by  the  side 
bipods,  top  bipods,  the  3-hole  tank  tab,  and  the 
top  and  bottom  Siamese  tank  tabs,  corresponding 
to  groups  A through  E in  the  perturbed  analytical 
model,  Table  2,  The  kinetic  energy  in  the  first 
5 modes  was  dominated  by  the  fluid  masses. 

Numerous  combinations  of  modes  and  member 
and  mass  groupings  were  used  in  the  solution  of 
Eqs.  (16)  or  (17).  The  groupings  allowed  for 
variations  in  (1)  stiffness  only,  (2)  mass  only,  and 
(3)  both  stiffness  and  mass.  None  of  the  trials 
produced  physically  reasonable  results.  The 
solutions  t>T)ically  would  indicate  variations  in 
member  stiffnesses  of  in  excess  of  il00%  and 
fluid  masses  of  ±30%.  Use  of  Eq.  (28)  did  not 
improve  the  results.  None  of  these  variations 
were  considered  physically  reasonable.  The 
systematic  search  for  the  perturbation  factors 
using  Eqs.  (16)  and  (17)  was  hence  abandoned. 

Instead  of  a systematic  search,  the  following 
procedure  was  used:  one  group  of  stiffness  or 
mass  elements  was  varied  at  a time  trying  to 
match  the  frequencies  of  from  3ne  up  ♦o  four 
normal  modes. 

The  group  selected  had  the  maximum  energy 
content  in  these  modes.  The  factor  for  the  member 
group  under  consideration  was  calculated  using  a 
least  squares  fit  approach,  which  is  equivalent 
to  solving  Eq.  (17)  with  one  unknown.  Consequently 
other  member  groups  were  rsed  to  improve  the 
correlation. 

Using  the  procedure  described  above,  two 
improved  models  were  obtained. 

Model  A:  Started  with  variations  in  the 

member  stiffnesses.  Once  this 
was  exhausted,  mass  variations 
were  used;  twelve  iterations  were 
required. 

Model  B:  Started  with  variations  in  the  fluid 
masses.  Once  this  was  exhausted, 
member  stiffness  variations  were 
used;  seven  iterations  were  required. 
For  either  case  the  criteria  for  convergence  was 
the  minimization  of  the  root  sum  square  (RSS)  or 
the  frequency  difference  between  the  test  model  and 
tlie  analytical  model. 

Table  5 shows  the  frequency  comparison  of  the 
various  models  and  the  test  data.  Table  6 gives 
the  perturbation  factors,  defined  earlier,  that  are 
required  to  obtain  the  frequency  improvement.  It 
now  remains  to  be  seen  if  the  model  has  indeed  been 
Improved,  Since  frequency  by  itself  is  not  a 
measure  of  model  correlation,  the  mode  shape, 
cross  orthogonality,  kinetic  energy  distribution, 
effective  weight  (Refs.  (11)  and  (12)),  and  force 
coefficients  of  the  first  6 normal  modes  will  b^ 
examined.  Table  7 compares  the  modal  deflec- 
tions for  the  major  masses  of  the  test  to  the 
analytical  deflcctionji  of  the  original  analysiSi 
Model  A and  Model  B. 
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Table  6 - Perturbation  Factors 


riement 

Model  A 

Model  n 

SUflncts 

Side  bipods 

1.307 

1.227 

Top  blpuds 

1 328 

1.000 

3'hole  tabs 

1.170 

1.089 

Top  Siamese  tabs 

0 754 

0.865 

Prcssur.mt  i.inK  supiwrt  round 

1.401 

1.685 

Pressuranl  UnK  support  square 

0 813 

0.813 

PC  A support  struts 

0.996 

1.000 

Shear  tie 

0.218 

0.19C 

Mass 

Oxidizer  lateral  mass 

1 000 

0.933 

Oxidizer  lontiiitudlnal  mass 

1.000 

0.900 

Fuel  later.d  ni.iss 

1.000 

0.933 

Fuel  longitudin.'il  mass 

1 000 

0.936 

Pressurant  Control  Assembly  Inertia  1^ 

3.199 

J.199 

Pressurant  Control  Assembly  Inertias  ly, 

1.050 

1 050 

Table  8 gives  the  cross-orthogonality  matrix 
between  the  first  six  test  modes  and  the  correspond- 
ing analytical  modes  for  the  three  analytical 
models.  The  cross  orthogonality  matrix  (CO)  is 
defined  as 


(30) 

% 


where 

prp  the  test  modes  or  eigenvectors 

^A  analytical  modes  or  eigenvectors 

L.p  Choleski  decomposition  of  the  test  mass 
matrix,  M*p 

Choleski  decomposition  of  the  analytical 
mass  matrix,  M^ 

thus 


Mt  “ ; 


(31) 


The  local  kinetic  energy  distribution  (Ref.  (13))  of 
the  major  masses  for  the  first  six  modes  is  listed 
in  Table  9. 

Table  10  lists  the  effective  weight  (Refs.  (11) 
and  (12))  as  a percentage  of  the  total  weight  for 
the  test  and  the  analytical  models. 

For  purposes  of  member  load  calculati  ns,  the 
force  coefficients  are  the  most  important  correla- 
tion parameter.  The  force  coefficients  for  the  main 
load  carrying  members  for  the  first  five  normal 
modes  are  shown  in  Table  11. 


DISCUSSION  OF  RESULTS 


Examining  the  results  of  the  analysis/test  corre- 
lation, the  following  observations  can  be  inadt*: 

(1)  The  perturbation  technique  of  Ret.  (1)  did 
not  converge  to  produce  a single  unique  mutlel. 
Numerous  other  combinations  of  variations  in 
stiffness  and  mass  elements  could  produce  similar 
frequency  correlation  to  that  shown  for  Models  A 
and  B, 


T;iblt  7 


Mode  Shape  Coinparlaon 


M.ift  Point  and 

Original 

Model 

Moclct 

Original 

Model 

Mulel 

Direction 

lifi 

Analyeis 

A 

Teet 

An.tlysle 

■ 

B 

2isna 

Mode  4 

X 

-0  063 

0.014 

0.014 

0.013 

1.000 

0 96? 

0.974 

0.966 

Rncine 

Y 

1.000 

0.948 

0.923 

0. 926 

0.075 

-0.040 

-0  026 

-0.028 

Z 

o.ou 

0.052 

0.056 

0.036 

0.221 

0. 150 

0.104 

0. 118 

X 

‘0.004 

O.OOH 

0.008 

0.008 

0.418 

0.423 

0.459 

0 420 

OxidUer  tank 

Y 

0.  G34 

0,318 

0.391 

0.  379 

-0.032 

-0. 026 

-0  021 

-0.023 

Z 

0.009 

0.037 

0.038 

0. 038 

0.344 

0.390 

0.269 

0 337 

X 

0.0*0 

0.005 

0.003 

0.003 

0.407 

0.365 

0.409 

0.383 

Fuel  uuk 

Y 

0.4^8 

0.437 

0.386 

0.392 

-0.010 

-0.013 

-0.  005 

-0. 002 

Z 

0.031 

0.038 

0.031 

0. 030 

-0.090 

-0. 136 

-0.114 

-0.J48 

^ 

X 

0.007 

-0.008 

-0.011 

-0.011 

-0  322 

-0. 178 

-0. 006 

-0. 139 

Y 

‘0.131 

-0. 177 

-0. 193 

-0  IV7 

0.029 

I)  007 

0.004 

0.006 

Z 

0 Oli 

0.036 

0. 037 

0. 037 

0,243 

0. 126 

0.074 

0.093 

RS8  u(  difference 

0.134 

0. 179 

i.  179 

0.264 

0.462 

0.355 

Mode  4 

X 

•l.OOU 

-0.662 

-0. 290 

-0  432 

0.033 

-0. 043 

-O.OOj 

-0.007 

Cnflne 

Y 

-0  060 

-0.080 

-0. 008 

-0.014 

0.636 

0 425 

0 446 

0.442 

Z 

0. 809 

0,708 

0.471 

0.363 

0.039 

-0.012 

-0  005 

0 007 

X 

-0.682 

-0.472 

-0.298 

-0.378 

-0.014 

-0. 030 

U 025 

0.021 

ORldUcr  Unk 

Y 

-0.044 

-0.051 

-0.073 

-0.  079 

-0.435 

-0.336 

-0.2G7 

-0.  279 

Z 

0.787 

0.733 

0.643 

0.672 

-0.027 

0.033 

-0.046 

-0.043 

X 

-C. 

•0. 260 

-0.  IGl 

-0.213 

0.017 

-0.101 

-0.023 

-0.023 

Fuel  tank 

Y 

-0.036 

-0.239 

0.060 

0.061 

1.000 

0.661 

0.636 

0.644 

Z 

0. 626 

0. 673 

0.294 

0.440 

-0.030 

-C.173 

-0.031 

'0.030 

X 

-0.443 

-0.410 

-0. 333 

-6.497 

0.127 

-0.129 

0 097 

0. 086 

PrvMurttAt 

Y 

0.180 

0.048 

0. 030 

0.042 

-0.114 

-0. 109 

-0. 105 

•0  098 

1AM 

Z 

0.830 

0.742 

0.497 

0.390 

-0.078 

-0. 084 

-0,  037 

-0.055 

RS8  uf  difference 

0,  343 

1.152 

0.680 

0.4G7 

0.456 

0.446 

Mude  } 

^lodc  B 

X 

-0.007 

0.004 

-0.071 

-0.110 

-0.060 

-0.007 

-0.013 

-0.012 

engine 

Y 

0 

0.059 

0.005 

0.005 

-0.367 

-0  055 

-0. 042 

-0.042 

Z 

0.314 

0.212 

0.265 

0.242 

0.068 

-0. 007 

0.005 

0.007 

X 

0. 193 

0, 139 

0.080 

0. 123 

0.041 

-0.005 

0.003 

0 

Qnidiaer  lank 

Y 

-0. 040 

-0. 062 

-0.025 

-0.023 

0.205 

0.043 

0.031 

0.032 

Z 

-0. lOi 

-0. 100 

-0.071 

-0.079 

-0.062 

-0. 007 

-0. 009 

-0.011 

X 

0.228 

0.144 

0.129 

0.153 

0.077 

o.ooe 

0.  003 

0 007 

Fuel  tank 

Y 

0.073 

0. 109 

0.031 

0. 029 

-0. 528 

-0.  078 

-0. 059 

-0.  U60 

Z 

0 613 

0.435 

0.545 

0.494 

0.309 

0.015 

0.029 

0 031 

n 

X 

1.000 

0.792 

0.770 

0.757 

1.000 

0.167 

0. 172 

0. 1G8 

Prceiurunl 

Tnjdc 

Y 

0. 024 

-0.016 

0.011 

0.008 

0.212 

0.027 

0.016 

0.016 

Z 

0.266 

0 173 

0 253 

0.218 

-0  140 

0.007 

0.012 

0.012 

RSS  of  difference 

0.334 

0.505 

0.331 

1.074 

1.096 

1 098 

0?’'  ^ 
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(2)  The  quality  of  the  improved  models  is 
questionable.  The  only  clearcut  improvement 
achieved  is  in  the  frequency.  It  is  not  clear  if  it 
would  not  have  been  just  as  effective  to  retain  the 
masses  and  stiffnesses  of  the  original  model  and  to 
adjust  the  frequencies  (generalized  stiffnesses)  to 
match  the  test.  The  data  show  that  the  mode  shape 
correlation  of  the  original  model  with  test  data  is 
generally  better  than  that  of  the  improved  models. 
This  is  also  time  of  the  local  kinetic  energy,  the 
effective  freights , and  the  force  coefficients.  It 

is  interesting  to  note  that  cross-orthogonaiity, 
which  is  often  used  as  an  indication  of  correlation, 
is  an  inconclusive  Indicator.  The  cross- 
orthogonality between  the  test  modes  and  all  three 
analytic  modes  is  quite  good  and  apparently 
insensitive  to  variations  in  mode  shape. 

(3)  It  is  difficult  to  reconcile  the  perturbation 
factors  listed  in  Table  6 with  systematic  model 
errors.  Thus  the  required  Increase  of  the  bipod 
areas  by  40%  does  not  seem  physically  reasonable, 
rhese  are  axial  load  carrying  members  that 
should  not  contain  that  high  an  inconsistency.  The 
shear  tie  decrease  in  stiffness  by  80%  docs  have 

a physical  explanation.  During  ^e  modal  test  it 
was  found  tliat  the  shear  link  did  act  as  a nonlinear 


Ttbla  • - Kinetic  Cnarsy  Olatiibutlon. 
% of  Total  Kinetic  Energy 


Table  10  - EffccUvc  Weight  Comparleon. 
% of  Total  Weight  and  Inertia 


Maaa  Point  and 

Original 

Model 

Model 

Original  1 

Model 

Model 

Direction 

Teat 

Analyata 

A 

B 

Teet 

Analyala 

A 

B 

Model 

klode  2 

X 

28.6 

84  0 

46.0 

87.8 

OBtidUer  tank 

Y 

76.  t 

66.3 

74.6 

78.6 

Z 

50.1 

35.1 

16.7 

27.8 

X 

18.8 

16.  f 

81.8 

18.0 

Fuel  tank 

y 

16. 0 

24.8 

16.5 

19.1 

z 

0.7 

8.2 

1.8 

2.8 

ModeJ 

Mode  4 

X 

27.0 

19.3 

14.5 

19.8 

Oxldlxer  tank 

Y 

81.6 

33.3 

19.6 

20.6 

2 

28.6 

49.6 

69.5 

59.5 

X 

6.2 

8.6 

2.1 

8.0 

0 

l.t 

0 

0 

Fuel  Uidt 

Y 

64.2 

50.1 

61.4 

60.6 

z 

32.8 

22.0 

7.6 

18.6 

0.1 

3.7 

0 

0 

Preaaurimt 

tank 

2.7 

0 

4.8 

4.3 

2.8 

0 

8.9 

3.9 

Mode  5 

Mode  6 

X 

11.7 

11.3 

3.4 

6.1 

OxJdUer  tank 

Y 

0.4 

2.0 

0 

0 

2.3 

8.3 

1.9 

1.6 

Z 

3.S 

5.7 

2.4 

3.1 

X 

8.9 

6.4 

4.4 

6.7 

Fuel  tank 

Y 

0.6 

3.4 

0 

0 

6.3 

5.8 

3.5 

3.4 

z 

56.9 

56.6 

74.9 

66.4 

Preeeurant 

tank 

•y 

3.6 

0 

1.7 

1.9  61.6 

65.4 

89.1 

89.2 

j 

h 

1.2 

2.1 

1.9 

1.8 

Direction 

Teat 

Original 

Analyala 

Model 

A 

Model 

B 

Teat 

Original  1 

Analysis 

Model 

A 

Model 

0 

Mode  1 

Mode  2 

31.6 

44.7  I 

50,7 

46.5 

83,8 

64.9 

03.0 

82.2 

1.2 

1.0 

0 

0 

0 

0.8 

0.8 

0.5 

27.2 

15.5 

7 3 

10.5 

95.0 

96.0 

96.5 

96.6 

0 

0.3 

0 

0 

V 

80.3 

92.1 

95.1 

94.7 

17.4 

11.9 

16.1 

10.3 

Mode  8 

Mpde^ 

91.6 

81.7 

16.9 

20.0 

0 

1.4 

0 

0 

1.2 

2.4 

0 

0 

4.8 

2.4 

1.5 

6.4 

62.1 

74.9 

72.6 

72.9 

1.0 

0.3 

1 2 

0.8 

0 

0.3 

0 

0 

1.9 

0.6 

2.  7 

2.4 

V 

14.6 

6.5 

0 

9.2 

0 

0 

1.4 

1.1 

71. 1 

67.  1 

67.  1 

68  8 

Mode  5 

Mode  6 

23.7 

21.7 

10.4 

18.6 

2.4 

0 

U 

0 

1.2 

0 

0 

0 

9.3 

6.8 

16.0 

13.2 

V 

0 

0 

2.4 

0 

0.7 

4.9 

0 

0.8 

7.7 

8 4 

3 8 

4.1 

spring  due  to  monoball  chatter.  Modal  test  data 
was  obtained  in  the  linear  region.  Such  a non- 
linear phenomenon  does  indeed  reduce  the  effective 
spring  constant,  and  it  was  found  that  the  shear 
link  had  to  be  reduced  in  stiffness  by  80%  to 
match  the  test  results  of  the  fourth  mode,  the 
only  mode  which  was  affected  bj  his  member, 

A mass  decrease  of  up  to  10%  as  shown  for 
Model  B is  physically  reasonable  due  to  possible 
fluid  slosh  effects,  even  though  the  tanks  were 
completely  full  and  pressurized.  The  increase 
in  the  moment  of  inertia  of  the  Pressurant  Control 
Assembly  (PCA)  about  the  X axis  by  300%  is 
indeed  acceptable  because  the  model  was  found  to  be 
in  error  in  the  PCA  moment  of  inertia  and  the 
calculated  increase  is  well  within  the  physically 
meaningful  range. 

(4)  The  data  presented  shows  that  it  would  be 
desirable  to  correlate  using  mode  shape  in  addition 
to  frequency.  Since  the  mode  shape  is  a vector 
rather  than  a scalar,  this  approach  might  be  more 
fruitful.  The  use  of  mode  shape  data  as  suggested 
in  Ref.  (1)  is  impractical.  Refs.  (14),  (15),  and  (16) 
use  frequency  and  mode  shape  for  model  improve- 
ment. This  method  requires  the  evaluation  of  the 
derivatives  of  the  eigenvalues  and  eigenvectors. 

In  addition,  it  considers  statistical  error  distribu- 
tions for  these  test  data  and  analyses.  The  method 
is  successfully  applied  to  a beam  model  of  the 
Saturn  V space  vehicle.  The  feasibility  of 
applying  the  method  proposed  by  Refs.  (14),  (15), 
and  (16)  to  a three-dimensional  structure  such  as 


Table  11  - Force  Coefficients,  Newtons 


Member 

No. 

Teet  ve. 
Analysis 

1 

-2527.5 

2 

343.0 

Mode  No, 
3 

-1183.7 

4 

-237. 1 

4 

Test 

-194.4 

Original  Analysis 

-2.'i48  4 

502.2 

-ior.5  6 

-491  5 

-iHfi.  4 

Model  A 

-2977.8 

768.4 

-535.5 

-321  2 

-437  1 

Model  li 

-2930  3 

689  4 

-Kyti  8 

-307  2 

-324  9 

3 

Test 

-297  1 

-387.4 

379  4 

402.6 

327  8 

Original  Analysis 

-328.3 

-488.0 

323  4 

552  0 

433  I 

Model  A 

-303.7 

-619.8 

,7.0 

34  5 2 

543.  1 

Model  B 

-332.0 

-629.  1 

1m7  1 

351.:i 

4ho.  5 

41 

Teat 

<120  0 

-1123.2 

-4H7  b 

9ir,  3 

-255. 3 

Original  Analyais 

1900.3 

-1086.7 

-657  4 

815  8 

-137  5 

Model  A 

2089. 6 

-1338.8 

-528  4 

718  3 

-217.2 

Model  B 

2094  6 

-1289.  1 

-6.J8.8 

8.U.9 

-270.3 

40 

Teat 

915.9 

643.7 

365  2 

251  3 

-235  8 

Original  Analysis 

1108  5 

639  2 

425,3 

179  7 

-438. 1 

Model  A 

1328.9 

724  5 

474.  1 

24ti  7 

-342  1 

Model  B 

1280.0 

758  3 

500.0 

260  b 

-.353.2 

12 

Test 

2213.9 

367.4 

-13K'  9 

191.7 

237  1 

Original  Analysis 

2442.5 

426.  1 

-1294.0 

279  8 

Model  A 

2845.2 

723.6 

-742.0 

r,(<  1 

4.’b  4 

Model  B 

2799.  7 

617.6 

-1174.3 

73.9 

-357  4 

11 

Test 

673.5 

-442.2 

355  4 

-43,’  7 

440  H 

Origin.!!  Analysis 

678,8 

-481.7 

258  4 

-It, 4 4 

5'0..3 

Model  A 

797  0 

-57y.3 

5JH.3 

-312  7 

t 5 

Model  B 

744.1 

-604  5 

-J64.  1 

-323  b 

9 j2  5 

the  Viking  Propulsion  Subsystem  has  not  been 
established. 

CONCLUSION 

The  application  of  the  method  proposed  !)v 
C.  Wliite  to  the  VPS  providetl  the  authors  with  a 
good  insight  of  the  method’s  merits  ami  hnutatiuas. 

A comparison  of  this  mcUu/l  to  otlu‘r  published 
methods  has  not  Ixjcn  made  because  other  mctliods 
have  not  Ixjcn  tested  as  yet. 
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The  lai:ger  objective  of  this  work  is  to  obtain 
insight  into  the  merits  and  limitations  of  various 
methods,  and  to  select  one  of  these  for  an  algorithm 
to  be  used  for  a computer  program  for  model 
correlation.  The  program  is  to  be  comparable 
with  NASTRAN. 

Activities  described  herein  are  invaluable 
in  a research  and  development  activity  leading 
to  a user-oriented  program.  Especially  since 
the  proposed  methods  in  the  literature  have 
been  applied  to  a very  limited  number  of  real 
problems. 

The  use  of  analytically  generated  "test"  results 
to  Illustrate  different  methodologies  can  be  mis- 
leading. Algorithms  that  are  successful  on 
analytically  generated  "test"  results  may  not  be 
successful  on  real  test  data. 


*Thls  paper  presents  'the  results  of  one  phase 
of  research  carried  out  at  the  Jet  Propulsion 
Laboratory,  California  Institute  of  Technology, 
under  Contract  No.  NAS  7-100  sponsored  by  the 
National  Ae^’cnautlcs  and  Space  Administration. 
The  effort  was  supported  by  Dr.  A.  Amos, 
Materials  and  Structures  Division,  Office  ot 
Aeronautics  and  Space  Technology,  National 
Aeronautics  and  Space  Administration. 
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MacNeal’s  Method 
of 

Component  Representation  for  Modal  Synthesis 
As  Employed  at  Rockwell  International/Space  Division 
for  Shuttle  Payload  Anal:,ses 

by  M.  A*  Martens 


The  shuttle  liftoff  or  landing  vehicle,  exclusive  of  the  payload, 
is  treated  as  one  component  which,  for  the  synthesis  of  system  modes, 
is  represented  in  the  form  described  by  MacNeal,  Each  payload  compo- 
nent is  represented  in  any  convenient  form  which  contains  the  orbiter- 
to-payload  boundary  degrees  of  freedom  in  physical  coordinates  (usually 
Craig-Bampton  or  physical  stiffness  and  mass), 

Rubin’s  extension  to  MacNeal’s  method  is  not  used  because,  for 
th€  specific  problem  at  hand,  improvement  in  results  sufr'icient  to 
justify  the  increased  cost  have  not  been  observed  in  our  investigations. 
However,  Rubin’s  paper  (reference  1)  provides  a lucid  explanation  of 
the  basic  method  and  should  be  referred  to  for  a more  thorough  deve- 
lopment of  this  method.  The  equation  numbers  in  this  paper  refer  to 
the  corresponding  equations  in  the  reference  paper,  though  the  inertial 
and  dissipative  terms  have  been  deleted  here. 

Basically,  the  component  is  described  by  a subset  of  its  free- 
free  modes.  However,  to  partially  account  for  the  degradation  in  the 
mathematical  representation  caused  by  modal  truncation,  the  flexibi- 
lity terms  corresponding  to  the  discarded  modes  are  computed  and  in- 
cluded in  the  component  representation.  These  residual  flexibility 
terms  (Gp)  are  found  as  the  difference  between  the  total  flexibility 
(G)  and  the  flexibility  matrix  (GN)  obtained  from  the  generalized 
stiffness  matrix. 


Gp  = G - GN 


(24) 


Since  the  component  is  in  an  unconstrained  state,  special  considera- 
tion must  be  given  to  obtain  the  "freed”  total  flexibility  matrix. 

It  has  been  shown  that  it  may  be  computed  from  the  rigid  body  modes 
(0R)>  physical  mass  matrix  (M),  and  physical  flexibility  matrix  (Gc), 
Gc  is  obtained  by  inverting  the  physical  stiffness  matrix  with  any 
arbitrary  set  of  statically  determinant  constraints  applied.  Assuming 
the  component  eigenvectors  are  normalized  to  unit  generalized  mass, 
the  total  flexibility  may  be  found  as: 

G = A^GcA  (21) 


(in 


where 


A = I . M0r0^ 
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The  Gj^  matrix  may  be  found  by  inverting  the  non-zero  diagonal  parti- 
tion (Kn)  of  the  generalized  stiffness  matrix  and  transforming  it 
back  to  physical  coordinates. 

GN  = (24) 

The  N partitions  of  the  generalized  stiffness  and  modal  matrices 
refers  to  those  partitions  corresponding  to  all  retained  flexible 
body  modes. 

With  these  equations,  Gp  may  be  computed.  The  component  may 
now  be  represented  with  a mass  and  stiffness  matrix  expressed  in  a 
coordinate  system  consisting  of  the  selected  free-free  component 
modes  (Q)  and  the  physical  boundary  degrees  of  freedom  such 

that  the  equation  for  eigenanalysis  is  as  follows: 


a 

' Ub 

Ma[  0 

0 I 0 


where 

Mq  = I (assuming  modes  normalized  to  unit 
generalized  mass) 


Kq 

0b 


= Generalized  stiffness 

= Partition  of  modal  matrix  containing  the  rows 
corresponding  to  boundary  freedoms  for  all 
retained  modes. 


G_  = Boundary  partition  of  residual  flexibility 
matrix  (square,  symmetric). 

This  component  representation  is  combined  with  the  equations 
for  all  other  components,  adding  equations  in  the  same  unknown  (Ut,) 
in  direct  stiffness  fashion.  The  boundary  loads  sum  to  zero.  Eigen- 
values ( y\*)  and  eigenvectors  (Q*,  Ub*)  extracted  in  the  usual 

manner  from  the  system  equations.  At  this  point,  the  transient  solu- 
tion could  be  performed  by  passing  the  non-boundary  applied  loads 
through  two  coordinate  transformations  and  the  boundary  loads  through 
one.  The  residual  flexibility  terms  would  appear  in  the  transforma- 
tions. A more  convenient  approach  is  to  transform  the  system,  modes 
back  to  the  physical  coordinate  system  before  the  transient  solution. 
In  order  to  obtain  the  non-boundary  terms  (subscript  i = interior  or 
non-boundary)  of  the  system  modes  in  the  physical  coordinate  system, 


I 


i 
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the  residual  flexibility  gives  rise  to  additional  terms  in  the  trans- 
formation equation: 

= 0iQ*  + Gp.j^  J2(bQ*  j (32  & 33) 

Once  the  U^*  and  U^*  portions  of  the  system  modes  are  available,  the 
transient  solution  may  proceed  exactly  as  if  these  modes  had  been 
obtained  without  the  use  of  a synthesis  method. 

In  order  to  employ  this  method  of  component  representation, 
inspection  of  equations  30,  32,  and  33  reveals  that  the  only  component 
data  required  are  the  retained  free-free  modes,  the  corresponding 
generalised  stiffness  and  mass  matrices,  and  the  columns  of  the  resi- 
dual flexibility  matrix  that  correspond  to  the  boundary  degrees  of 
freedom.  Considerable  savings  can  be  achieved  in  the  multiplication 
operations  of  equations  21  and  24  by  computing  only  the  required 
columns  of  the  residual  flexibility  matrix. 


Reference  1:  Rubin,  S.,  "Improved  Component-Mode  Represonxation  for 

Structural  Dynamic  Analysis,"  AIA.\  Journal,  Volume  13, 
Number  8,  August  1975,  pp.  995-1006. 


M.  Martens,  6/78 
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A method  has  been  developed  by  S.  R.  Ibrahim  of  Old  Dominion  University,  Norfolk 
Virginia,  for  analyzing  free-decay  data  for  complex  structures.  The  method  analyzes 
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COMPARISON  OF  MODAL  FREQUENCIES 


IDIOQ. 

flPSD 

SPR-JPI 

SPR-LRC 

1 

10.5-10.7 
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2 

10.9-11.3 

11.29 
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17.1-17.2 
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17.18 
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5 
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22.6-23.5 
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11 
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1^1 
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2i\ 
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10.95  1st  BENDING.  Y 

11.19  1st  bending.  X 

17.25  1st  torsion 

21.80  RTG 

23.12  2nd  torsion  + LECP 

23.30  2nd  torsion 

29.10  2nd  torsion  + MAGNETOMETER 

25.13  2nd  bending.  X 

26.05  2nd  bending.  Y 

27.81  SCAN  PLATFORM.  AXIAL 

30.92  ANTENNA.  2 

31.56  3rd  BENDING.  Y + ANTENNA 

32.50  3rd  bending.  X 

32.98  3rd  torsion 

33.96  LECP.  RTG  TORSION 

39.79  SCAN  PUTFORM.  X 

35.99  RTG.  X 

36.81  SCIENCE  APPENDAGE.  X 

38.75  SCIENCE  APPENDAGE.  X.  ANTENNA 

39.19  SCIENCE  APPENDAGE.  X.  2 

ROTATION 

90.29  SCIENCE  APPENDAGE.  X 

90.69  SCIENCE  APPENDAGE.  X 

92.21  ANTENNA  ROTATION.  X 

99.59  TANK  ROT..  SCIENCE  APPENDAGE 

99.80  HYDRA2ENE  TANK  ROTATION 

95.65  1st  axial 

98.09  TANK.  APPENDAGE 

50.69  RTG.  SCAN  PLATFORM.  2 
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A COMPARISON  OF  TEST  TECHNIQUES  USED  DURING 
MODAL  TESTING  OF  ET  LOX  TANK 

G,  0.  Johnston,  Houston  M.  Hammac,  A.  D.  Coleman,  NASA/MSFC 


This  presentation  will  briefly  show  the  results  of  test  data  obtained  from 
the  modal  test  of  the  ET  Lox  Tank.  A comparison  will  be  made  of  data 
obtained  using  Multi-Point  Sine  tuning  of  modes  with  data  obtained  from 
Single-Point  Random  tests  for  the  same  condition.  For  the  propose  of 
this  presentation,  only  the  liftoff  condition  will  be  used  to  make  this 
comparison.  A more  detailed  comparison  of  the  remaining  conditions  and 
description  of  the  test  setup  will  be  made  in  the  published  article. 

This  test  vvas  a major  link  in  the  verification  chain  of  hydro-elastic 
analysis  and  test  technique,  and  the  results  are  being  used  to  mathematically 
predict  the  modal  characteristics  of  the  entire  Exteinal  Tank.  The  general 
objective  of  this  test  program  was  to  determine  the  symmetric  and  anti- 
symmetric hydro-elastic  modal  properties  of  a fliglit  configuration  test 
article.  Specifically,  the  program  was  performed  lo  experimentally 
determine  vibrational  frequencies,  mode  shapes,  damping,  and  modal  energy 
distribution  of  targeted  modes  selected  from  the  pre-test  analysis. 

The  test  article  consisted  of  a Support  Ring,  Intertonk,  and  Lox  Tank, 

Figure  1 shows  the  test  article  in  Condition  I.  The  supporting  special 
test  equipment  consisted  of  the  Access  Structure,  Air  Suspension  System, 
Chromate  Water  Transfer  System,  the  Pressure-Purge  and  Vent  System, 

Shaker  System,  Digital  Control  System,  and  Data  Acquisition  system.  The 
Access  Structure  and  Air  Suspension  System  are  shown  in  Figure  1.  A 
closeup  view  of  one  air  bag  pad  is  shown  in  Figure  2.  Figure  3 is  a 
schematic  of  the  air  suspension  system.  Figure  4 is  a schematic  of  the 
Chromate  Water  Transfer  System.  Sodium  niono-chromate  was  used  as  a 
corrosion  inhibitor.  The  concentraLion  was  maintained  between  200  and 
350  parts  per  million.  Figure  5 is  a schematic  of  the  Pressure-Purge  and 
Vent  System.  The  ullage  pressure  was  maintained  at  3.3  psig  + 1.7  psig 
throughout  the  test  except  when  the  1.6  psig  and  8.0  psig  condition  were 
performed.  Figure  6 shows  the  Digital  Control  System  and  Figure  7 is  a 
portion  of  the  Data  Acquisition  System. 

One  major  feature  of  this  test  setup  was  the  ability  to  cant  the  test 
article  as  much  as  thirteen  degrees  from  vertical.  This  discussion  will 
be  limited  to  Condition  I (0°  cant  and  487  inches  of  fluid  with  3.3  psig 
ullage  pressure).  All  of  the  test  conditions  where  modal  data  were 
obtained  are  shown  In  Table  I. 

A half-tank  finite  element  model  of  the  test  article  was  generated  for  the 
analyses  by  the  Martin  Marietta  Corporation.  The  choice  of  the  half-tank 
was  permitted  by  a convenient  plane  of  symmetry  and  was  required  by  the 
time  and  size  constraints  in  the  computer  analysis  that  was  performed. 

A detailed  listing  of  the  modal  grid  and  finite  elements  used  can  be 
found  in  Appendices  A and  B of  the  MMC  report  number  MMC-ET-SE2 I -5  dated 
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October  6,  1978,  The  model  incorporated  the  structural  elements  of  the 
Lox  Tank,  intertank,  load  ring,  air  support  system,  and  the  fluid.  The 
finite  element  computer  programs  used  are  part  of  the  Martin  Marietta 
Aerospace  Library  known  as  FORMA.  A detailed  description  of  the  modeling 
methods  used  and  analysis  procedure  are  found  in  Sections  4.0  and  5*0  nf 
the  MMC  report  number  MMC**ET^SE2 1 *5  dated  October  6,  1978> 

The  Test  Requirements  Specification  (MMC-ET-TM07)  specified  thirteen 
shaker  positions  as  shown  in  Figure  8 > In  past  modal  tests,  it  has  been 
the  policy  of  MSFC  to  select  force  input  points  at  those  locations  where 
mechanically  induced  energy  is  experienced  by  the  structure  to  obtain 
direct  transfer  functions.  An  additional  shaker  position  was  added  at  the 
+ Y Solid  Rocket  Booster  attach  point.  This  shaker  position  was  selected 
to  be  the  principal  data  point  for  the  single  point  random  excitation  at 
the  most  critical  condition  of  testing.  The  reasons  for  this  decision 
were  that  1)  it  would  provide  more  direct  transfer  function  data 
representing  the  SRB  input,  2)  it  would  provide  a more  complete  data 
bank  of  modes  for  future  analysis  that  possibly  could  be  overlooked  during 
multi-point  sine  testing,  3)  it  would  provide  an  excellent  opportunity  to 
compare  the  two  techniques  and  possibly  improve  on  the  current  capability 
of  single-point  random  testing.  Also,  it  was  deemed  necessary  to  incorporate 
the  single-point  random  capability  because  of  the  short  period  of  time 
allowed,  initially,  for  the  entire  test  program.  A total  of  three  ninths 
was  scheduled  for  all  four  conditioris.  It  was  reasoned  that  if  one  long 
delay  was  encountered  after  the  test  Wsis  started,  single  point  random 
testing  might  be  the  only  way  data  could  be  obtained  for  some  of  the  final 
conditions.  Single-point  random  data  was  obtained  at  three  different  input 
points  during  Conditions  II  and  III.  The  driving  points  were  at  the  ogive 
tip  in  the  Z and  Y directions  and  at  the  SRB  c»'05sbeam  attach  point  on  the 
+ Y side  of  the  Lox  Tank. 

This  test  did  have  a few  ‘'firsts'*  for  MSFC.  One,  of  course,  is  the  first 
use  of  single-point  random  excitation  to  obtain  modal  data  on  a large 
scale  structure.  Another  “first"  was  the  canting  of  a large-liquid  filled 
tank  to  4 and  13^  angles  to  obtain  modal  data.  The  cant  angle  did  have 
dramatic  effects  on  bulge  modes  and  bending  modes.  These  effects  are 
covered  in  more  detail  in  the  test  reports  by  MSFC  and  the  MMC. 

A modified  Hewlett  Packard  Model  54513  digital  system  was  used  to  perform 
the  multi-point  sine  test.  This  system  is  shown  in  Figure  6 and  contains 
a mini -computer  and  a fourier  analyzer.  It  has  the  capability  of  controlling 
frequency,  phase  and  amplitude  of  up  to  eight  shakers  simultaneously.  It 
contains  a 32  channel  multiplexer  to  allow  on-line  recording  of  32  force 
and  acceleration  measurements.  At  the  beginning  of  a test  condition,  a 
wideband  sweep  would  be  performed  at  a single  shaker  position.  From  this 
sweep  31  transfer  function  plots  were  made  to  identify  modal  frequencies. 
These  plots  were  of  the  real  and  imaginary  values  of  each  accelerometer . 

This  process  was  repeated  at  several  shaker  locations  to  assure  all  modes 
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were  identified.  Narrowband  sweeps  in  frequency  steps  of  ,01  hertz  were 
performed  where  very  high  modal  density  was  evident.  Based  on  the  sweep 
data  and  plots  of  the  on-line  data,  shaker  positions  were  selected  to  tune 
the  targeted  modes.  To  tune  the  modes,  on-line  measurements  were  monitored. 
Several  parameters  were  used,  but  the  principal  ones  were  the  Co-Quad 
values  and  the  phase  angle  of  these  measurements.  Lissajous  and  the 
driving  point/acceleration  phase  angle  was  used  in  some  cases.  When  the 
mode  was  tuned^all  2Q2  accelerometers , driving  point  forces  and  fluid 
pressures  were  recorded  on  the  Structural  Data  Acquisition  System  shown 
in  Figure  7.  After  the  data  was  recorded,  a soft  shaker  dump  was  used  to 
obtain  modal  decay  data.  The  on-line  measurements  were  recorded  by  the 
HP  545IB  system  and  the  damping  calculations  were  made  by  curve  fitting 
these  decays  automatically. 

Single-point  random  tests  were  performed  at  the  end  of  each  test  condition. 
Excitation  for  the  SPR  testing  was  provided  by  a Hewlett-Packard  5^25 
vibration  control  system.  The  drive  spectrum  was  a shaped  5 to  50  hertz 
bandwidth  ranging  in  composite  force  from  93  to  150  RMS  Force-pounds. 

Data  was  acquired  with  the  same  Hewlett-Packard  5^5 IB  modal  system  used 
to  perform  the  MPS.  Using  the  32  chanfiel  multiplexer,  all  2u2  measurements 
were  recorded  by  seven  sequential  patchboards  thru  a patching  matrix. 
Approximately  thirty  minutes  of  data  were  recorded  of  each  measurement 
and  stored  on  magnetic  tape  for  later  analysis. 

The  data  was  processed  and  analyzed  employing  a least  squares  curve  fitting 
algorithm  to  obtain  the  mode  shapes  and  modal  coefficients.  Tables  II  and 
III  list  the  modal  frequency  and  damping  obtained  from  multi -point  sine 
(MPS)  and  single-point  random  (SPR).  It  should  be  emphasized  that  due 
to  insufficient  time  in  the  schedule,  only  one  excitation  point  for  single- 
point random  was  used  to  obtain  data  for  Condition  I,  The  only  target 
modes  not  well  matched  with  the  mult  I -point  sine  data  are  three  symmetric 
modes.  All  three  unmatched  modes  have  bending  in  the  symmetric  or  Z 
Plane.  Based  on  some  of  the  analysis  of  data  from  Condition  II,  where 
three  excitation  points  were  used,  the  modal  matching  is  even  better, 

A lot  of  shell  modes  are  listed  in  the  SPR  column  that  do  not  appear  in  the 
MPS  column.  Again,  due  to  insufficient  time,  the  numerous  shell  modes  could 
not  be  obtained  using  the  MPS  technique.  This  does,  however,  point  out  the 
great  advantage  of  the  SPR  technique.  The  test  article  is  now  being  installed 
for  the  static  loads  test,  but  we  can  continue  investigating  all  the  inodes 
in  the  tank  by  curve  fitting  the  data  stored  on  tape. 

The  SPR  data  indicates  excellent  correlation  with  the  MPS,  The  damping 
values  agree  very  well  and  tend  to  verify  that  for  this  test  condition 
the  l.ox  Tank  is  a lightly  damped  structure.  It  must  be  realized  that 
many  of  these  modes  will  change  and  have  considerably  higher  damping 
when  the  Lox  Tank  becomes  an  intergral  part  of  the  Shuttle.  There  are 
two  modes  listed  here,  however,  that  will  not  be  affected  significantly 
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in  the  total  assembly.  In  the  opinion  of  the  authors,  the  second  nd 
third  bulge  modes  of  the  aft  dome  (12.76  and  18.95  hertz)  will  still  have 
very  low  damping.  These  two  modes  are  always  the  ones  that  present 
problems  from  the  standpoint  of  'POGO*  and  loads  analysis.  The  second 
dome  bulge  mode  was  the  strong  POGO  initiator  in  the  Saturn  V and  Saturn 
IB  vehicles.  These  modes  should  be  observed  carefully  during  MVGVT 
testing. 

Appendix  A contains  a comparison  of  seven  mode  shapes  corresponding  to 
moda I matches  in  Tables  II  and  III.  A more  complete  set  of  mode  si  jpes 
will  be  published  at  a later  date  and  v/i li  contain  data  from  the  remaining 
test  conditions.  Condition  II  presents  some  very  interesting  data  because 
of  the  effect  of  the  13  degree  cant  angle.  The  math  model  did  a reasonably 
good  job  of  predicting  the  modes  at  the  13  degree  attitude,  but  much  more 
testing  should  be  performed  to  assist  in  improving  the  math  modeling 
techniques  for  in-flight  cant  angles. 
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TEST  CONDITIONS  PERFORMED 
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ET  LOX  MODAL  TEST  DATA 
COMPARISON  OF  FREQUENCY  AND  DAMPING 
MULTI -POINT  SINE  VERSUS  SINGLE  POINT  RANDOM 
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ET  LOX  MODAL  TEST  DATA 
COMPARISON  OF  FREQUENCY  AND  DAMPING 
MULTI -POINT  SINE  VERSUS  SINGLE  POINT  RANDOM 

CONDITION  I - 0°  CANT,  48?  INCHES  FLUID  LEVEL  (LIFTOFF) 
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FIGURE  5 
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SHUTTLE  P/L  BAY  ACOUSTIC  EMVIROHMENT  PREDICTION  ACTIVITY 
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EFFFCT  OF  DELTA  D PAYLOAD  ON  SUBVOLUME 
SPACE-AVERAGED  SOUND  PRESSURE  LEVELS 


One-Third  Octave  Band  Cef.Ter  Frequencies  in  Hz  ' os) 


COA/PARISON  OF  PREDICTED  AND  MEASURED  CHANGES  IN 
LEVELS  (SPACELAB  CONFIGURATION  2) 


SOUND 


Change  In  Sound  Level  When  Payload  Introduced,  dB 


Around  30%  Diometer  Cyllnd 


Predicted 

Measured 

99%  Confidence  Intervals 


Chonge  in  Sound  Level  When  Pc/lood  Introduced,  dB 
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UWllY 


Around  lO^o  Diameter  Cylinder 


Above  Pallet 


Above  95%  Diamete^Cylinder 
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Below  95%  DIometer  Cylindei  and  Pallef 


31.5  63  125  25 

1/4  Scale  Model  One-Third  Octave  Bond  Center  Frequencies  In  Hz  (cp$) 
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8 16  31.5  6: 

Shuffle  One-Third  Octave  Bond  Center  Frequencies  in  Hz  (cps) 

FINAL  COMPARISON  OF  PREDICTED  AND 
MEASURED  CHANGES  IN  SOUND  LEVELSL 
00  DELTA  D PAYLOAD  MODEL 


SHUTTLE  P/L  BAY  ACOUSTIC  ENVIRDNFIENT  PREDICTION  ACTIVITY 


MAJOR  MILESTONES 


• UPDATE  PREDICTION  MATH  MODEL  TO  EXTEND 
COMPUTATION  R.ANGE  DOWN  TO  12  Hz  1/3 

OCTAVE  BAND.  1/79 

• UPD.ATE  HATH  MODEL  TO  INCLUDE  LOW  FREQUENCY 

DYNAMIC  PROPERTIES  OF  SIDEWALL  STRUCTURE.  3/79 

• FINAL  VALIDATION  OF  MATH  MODEL  PRIOR  TO 
SS-1  FLIGHT  USING  .ALL  AVAILABLE  GROUND 

TEST  DATA.  5/79 

( COMPUTER  PROGRAM  RELEASE  - PACES  (VER.  2)  6/79 

• VALIDATION  OF  HATH  MODEL  USING  SHUTTLE 

FLIGHT  DATA  (SS-1  THRU  SS-N)  1/80  - 9/80 
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Ir  y 


SHUTTLE  P/L  BAY  ACOUSTIC  EWIRnrilTENT  PREDICTION  ACTIVITY 


CURRB1T  STATUS  ■ 
(1V78) 


• COHBINATIUI  OF  OV  101 AIC  VH  SCAlf  r<JEL  TEST  RESULTS  HAVE 
VALIIKTED  THE  BASIC  FORMULATIQI  OF  THE  PREDICTiai  miH  MODEL 
FOR  BOTH  THE  BPIY  Ait)  PAYLOAD  OCCUPIED  BAY  CBIDI I lOilS 

• VA  SCALE  f»Da  TEST  HAS  VALIDATED  AC  GIVBI  QUA'TTITATIVE 
EVIDBKE  TO  THE  .CTlQi  HIAT  SOf€  P/L  COiFIGURATIClS  CAN 
PWDUCE  SI&IIFICATT  liICREASES  OVER  THE  EiPIY  BAY  ACOUSTIC 
lEVElS 

■ THE  OV  101  AC  1/A  SCALE  TEST  RESULTS  ARE  CDISlSTEilT  IN 
IlETTIFYIiIG  sort  SI&NIFICAH  MATH  I-m  DEFICIEICIES  1,1 
PORTIQIS  OF  THE  LOW  FREfteiCY  REGIDI  SELaT  201  HZ.  THE 
LAB5E  MAJORITY  OF  THESE  DEFICIEICIES  ARE  IH  THE  UNCON  - 
SERVATIVE  DIRECTIQN. 

• ALL  TESTS  TO  DATE  HAVE  BEEi  CaCUCTED  WITHOUT  P/L  DOOR 
RADIATORS  PRESETT.  TriERlfORE  T.£RE  HAS  BEE!  iO  VALIDATIQT 
OF  THEIR  EFFECT. 
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HRUOA-BENFIELD  METHOD  WERE 


Truncated  component  mode  syntheses  were  performed  employing  (A)  30  DSP/IUS  modes  fixed 
at  the  orbiter  interface  plus  10  constraint  modes  associated  with  orbiter  interface  motion  and  (B)  170 
d.  o.  f.  space  shuttle  mode  truncations.  A simple  170  mode  space  shuttle  truncation  to  30  Hz  serves  as 
model  Bl;  and  a 160  mode  plus  10  orthogonalized  residual  mode  truncation  to  30  Hz  serves  as  model  B2 
for  demonstration  of  the  proposed  method.  Since  the  space  shuttle  model  was  supplied  originally  in  terms 
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"Vibration  Analysis  of  Structures  by  Component  Mode  Substitution"  W.'A.  Benfield  and  R.  F.  Hruda, 
AIAA  Journal,  Vol.  19,  July  1971,  pp  1255-1261. 
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FREOUENCY  COMPARISON 
CASE  B - 7tX  OF  AVAILABLE  MOOES. 

* OIFF  FOR  VARIOUS  MODAL  COUPLING  TECH 
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frequency  comparison 

CASE  P ■ 27X  OF  AVAILABLE  MOOES. 

% OIFF  FOR  VARIOUS  MODAL  COUPLINQ  TECH 
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MODE  SHAPE  COMPARISON 
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VERIFY  THE  SPACECRAFT  STRUCTURE  INTEGRITY 


LOADS  ANALYSIS  BY  SHOCK  SPECTRA 
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0 100  200  300  400  500 

© LIFT-OFF  time  AFTER  LIFT-OFF  (SECOND) 


FLIGHT  MEASURED  VOYAGER  SPACECRAFT 
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FLIGHT  LOAD  HISTORY  OF  MEMBER  71937 
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TIME  (sec) 


ACCELERATION, 


MAX/MIN  FLIGHT  AND  PREDICTED  LOADS 
COMPARISON  FOR  MISSION  MODULE  TRUSS 
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MAXIMUM  MEMBER  LOAD 
FOR  VARIOUS  EVENTS 


-o~ 

■ MEMBER  6(J011 

VOYAGER  A 

-Q— 

• MEMBER  68041 
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■ MEMBER  68011  | 

VOYAGER  P 
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■ MEMBER  68041  1 

LAUNCH  MAX  .q  STG  I IG  STG  I BO  S7G  II  BO  MECO 

DYNAMIC  [VhNTS 


VOYAGER  A INTERFACE 
^ ACCELERATIONS  FOR  TITAN 
STAGE  I BURNOUT  AND 
CENTAUR  MAIN  ENGINE  CUTOFF 


MECO  M 

STAGE  I BURN  OUT 
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TIME  (SECOND) 
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BOTH  SPACECRAFTS 
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TUNNEL  (STT) 


LOADS  METHODOLOGY  FOR  THE  SPACELAB  TRANSFER  TUNNEL 


/. 


The  Spacelab  Transfer  Tunnel  (STT) , located  within  the  orblter  payload 
bay,  provides  access  from  the  orblter  cabin  to  the  Spacelab  (SL)  ex- 
perimental area.  At  the  forward  end,  the  tunnel  adapter  connects  the 
STT  with  the  orblter  cabin.  At  the  aft  end  the  STT  attaches  to  the 
forward  cone  of  the  SL. 

The  configuration  shown  is  known  as  the  short  tunnel.  The  long  tunnel 
configuration,  indicated  by  dashed  lines  for  the  SL  forward  cone,  has 
an  additional  cylindrical  section  added  at  the  front  keeping  the  joggle 
section  next  to  the  SL. 

Design  and  construction  of  the  STT  has  been  contracted  to  the  McDonnel 
Douglas  Technical  Services  Company  (MDTSCO)  by  MSEC  with  MDTSCO  subcon- 
tracting the  flex  sections  to  Goodyear  Aerospace  Corporation  (GAC). 
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LOADS  METHODOLOGY  FOR  THE  SPACELAB  TRANSFER  TUNNEL  (STT) 


The  support  struts  and  keel  fittings  are  designed  by  the  STT  inertia  loads 
during  the  various  phases  of  flight.  The  forward  and  aft  flex  sections  are 
designed  by  loads  due  to  pressurization  and  deflections.  The  flex  sections 
are  Interchangeable  and  are  designed  by  the  envelope  of  loading  conditions. 

Positive  forces  and  moments  are  as  Indicated. 
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LOADS  METHODOLOGY  FOR  THE  SPACELAB  TRANSFER  TUNNEL 


The  basic  STT  design  consist  of  having  the  STT  response  to  Its 
environments  taken  out  through  the  STT  support  struts  and  the 
keel  fittings.  At  the  Orblter/STT  and  Spacelab/STT  Interfaces 
the  only  loads  will  be  those  due  to  pressure  and  the  reslstence 
of  the  flex  sections  to  Induced  deflections. 
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LOADS  METHODOLOGY  FOR  THE  SPACELAB  TRAHSFER  TUNNEL 


Dlsplacenents  and  pressure  requlreaents  for  designing  the  fortrard 
and  aft  flex  sections. 
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f1*f  NO 


AS  EACH  LOAD  PARAMETER  IS  EXPRESSED  AS  A FUNCTION  OF  PRESSURE  AND  FLEX  SECTIC.I  WIDTH 
THE  MOST  EFFICIENT  WIDTH  FOR  MINIMIZING  LOADS  CAN  BE  SELECTED. 
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SPACECRAFT  LOADS  EXPERIBCE  SUfTIARY 
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SPACECRAFT  LOADS  EXPERI3CE  SIM'WY,  CONTIfiUED 
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SPACECRAFT  LOADS  EXPERIENCE  SUtTORY  CONTINUED 
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SPACECRAFT  LOADS  EXPERIENCE  SUffARY,  COfD’INlJED 
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DATE  PROGRAM  OVERVIEW 
by 

W.  Brian  Keegan/William  F.  Bangs 
NASA/Goddard  Space  Flight  Center 


(Figure  1) 

Tne  Dynamic,  Acoustic,  and  Thermal  Environments  (DATE) 
Program  is  an  effort  to  provide  comprehensive  measurements 
of  three  potentially  critical  environments  for  payloads 
contained  in  thf'  STS  cargo  bay,  and  to  coordinate  the 
utilization  of  this  information  such  that  all  STS  payload 
developers  will  benefit  from  it  to  the  fullest  possible 
extent . 

The  DATE  program  was  conceived  because  of  a perceived  void 
in  the  attention  being  given  to  payload  environments  by 
the  STS  developer.  While  this  lack  of  attention  to  payloads 
was  understandable  (given  the  problems  associated  with  the 
development  of  such  a complex  vehicle  as  the  STS),  it  was 
nonetheless  apparent  that  the  payload  community  itself  would 
have  to  initiate  the  effort  to  obtain  the  data  it  felt  it 
needed  if  the  cost  benefits  associated  with  the  use  of  STS 
were  to  be  fully  realized. 

Thus,  the  DATE  program  was  formulated  and  a proposal  was 
made  to  NASA  Headquarters  in  October  1977,  that  the  DATE 
experiment  constitute  a portion  of  the  OEX  Program  whose 
primary  objective  is  "to  augment  the  research  and  technology 
base  for  future  aerospace  vehicle  design  by  utilizing  the 
Space  Shuttle  as  a research  vehicle  to  collect  data  in  all 
related  technology  disciplines,” 


(Figure  2) 

In  liiie  with  this  OEX  Program  objective  therefore,  the 
principle  objective  of  DAVE  is  to  develop  accurate  prediction 
techniques  for  payload  environments  in  each  of  the  three 
critica.l  areas  through  an  iterative  process  of  payload 
response  prediction,  followed  by  the  actual  meas  roment  of  the 
environments,  followed  in  turn  by  refinement  of  tue  prediction 
techniques  until  the  process  of  response  prediction  has  been 
verified  to  be  accurate  through  flight  measurement. 

While  the  STS  is  being  used  to  accomplish  this  objective,  the 
prediction  techniques  so  developed  would  obviously  not  be 
restricted  to  STS  payloads. 


(Figure  3) 


The  approach  to  be  followed  in  meeting  this  objective  is 
to  first  acquire  a set  of  baseline  measurements  that  consider 
the  effects  of  several  potentially  significant  parameters. 

These  include  payl  „d  mass  and  size,  particularly  as  the  size 
affects  the  clearance  with  the  cargo  bay  walls.  Studies 
performed  by  Bolt , Beranek  and  Newman  on  the  STS  cargo  bay 
acoustic  environment  have  predicted  that,  in  certain  cases, 
reduced  clearances  may  have  an  adverse  effect  on  the  acoustic 
noise  environment  encountered  by  payloads. 

Meanwhile,  the  variation  in  dynamic  loads  must  be  assessed 
because  they  will  be  influenced  by  the  payload  location 
within  the  bay  and  by  its  mounting  configuration,  that  is, 
method  of  attachment  to  the  STS  payload  support  structure. 

And  finally,  the  flight-to-f light  variations  brought  about 
by  trajectory  dispersions,  winds,  and  the  like  must  also 
be  measured. 

Simultaneous  with  these  baseline  measurements,  potential 
improved  methods  for  payload  environmental  prediction 
could  be  developed,  the  effectiveness  of  which  could  be 
evaluated  by  comparing  the  measured  responses  with  the  predicted 
ones.  These  analytic  methods  could  thus  be  refined  and 
re-evaluated  until  their  accuracy  has  been  verified,  at 
which  time  these  methods  could  be  utilized  by  ail  payload 
developers. 

The  DATE  program  is  being  managed  by  the  Goddard  Space  Flight 
Center  (GSFC).  However,  because  this  principle  objective  of 
DATE  has  such  broad  technology  implications,  a NASA-wide 
panel  of  experts  has  been  established,  known  as  the  DATE 
Working  Group.  It  provides  guidance  and  direction  for  the 
technology  development  permitted  by  the  DATE  measurements  and 
their  subsequent  application. 


(Figure  4) 


Of  necessity,  this  objectiv^^  of  improved  technology  is  a 
long  term  one.  There  are,  additionally,  some  near  term 
benefits  of  the  data  acquisition  sought  by  DATE  and  these 
are  best  summarized  by  the  secondary  objectives  of  characterizing 
the  STS  payload  dynamic.,  acoustic  and  thermal  environments  and  of 
deve" oping  a set  of  design  and  test  criteria,  directly 
aoplicable  -o  STS  payloads  in  order  to  permit  design  optimization 
so  as  to  belter  utilize  the  full  cwpabilif'es  jf  the  STS. 


(Figure  5) 


Because  of  practical  considerations  necessitated  by  the 
STS  launch  schedule,  the  ’lATE  measurement  program  has  been 
subdivided  into  two  phases.  Phase  I,  while  containing 
fewer  instrumentation  channels  than  is  truly  desirable, 
will  permit  some  data  to  be  obtained  during  the  §TS 
orbital  flight  test  sequences,  namely  on  mission  4,  on  which 
the  payload  will  be  the  Geostationary  Operational  Environmental 
Satellitf'-D  (GOES-D),  a free-flier  that  will  be  boosted  to 
synchronous  altitude  by  the  Solid  Spinning  Upper  Stage  (SSUS) 
and  on  n,'’'.sion  5 in  which  the  payload  will  be  a single 
pallet  that  will  remain  attached  to  the  STS  for  landing  as 
well.  T-,  o Phase  I data  will  indeed  begin  to  characterize 
the  STS  payload  environments  for  the  two  basic  types  of 
payloads,  and  will  serve  ar  a basis  from  which  to  refine  the 
instrumentation  plans  for  specific  follow-on  missions. 

Perm  .ssion  has  been  received  from  the  STS  Project  Office  at 
Johnson  Space  Center  (JSC)  to  utilize  the  development  flight 
instrumentation  system  to  record  this  Phase  I data,  barring 
any  undue  complications  resulting  from  the  earlier  STS 
orbital  flights. 


(Figure  6) 


The  DATE  Phase  II  program  will  consist  of  a far  more 
extensive  set  of  measurements  on  approximately  9 missions 
during  an  18-24  month  period  after  the  STS  becomes  operational. 

For  this  phase,  the  measurement  complement  would  be  expanded 

to  include  force  gages  and  thermal  measurements.  Recording  of  this  data  woul 

utilize  the  Technology  Flight  Instrumentation  System,  which 

will  be  developed  by  the  OEX  Project  Office  for  joint  use 

by  all  OEX  funded  experiments.  The  Phase  II  instrumentation 

plan  shown  here  is  intended  to  conceptualize  the  program. 

Specific  payloads  have  not  yet  been  assigned  to  all  these 
missions  and  the  actual  missions  which  are  instrumented  and 
the  precise  number  of  transducers  used  on  each  will  undoubtedly 
vary  as  the  program  develops. 


(Figure  7) 


While  the  payload  shown  here  is  only  conceptual,  it  can  be 
used  to  illustrate  typical  locations  that  may  be  selected 
for  making  measurements  during  the  DATE  program.  Microphone 
locations  would  be  selected  to  measure  the  spatial  variation 
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througnout  the  cargo  bay,  particularly  as  influenced  by 
the  affects  of  the  payload  configuration.  High  frequency 
accelerometers  would  emphasize  component  responses,  induced 
by  both  payload  bay  acoustics  as  well  as  structure  borne 
random  vibration.  Additionally,  an  attempt  would  be  made  to 
determine  the  magnitude  of  the  random  vibration  directly 
transmitted  to  the  payload  by  the  STS.  Low  frequency 
accelerometers  would  attempt  to  measure  the  forcing  functions 
at  the  payload/STS  interface,  as  well  as  responses  at 
critical  structural  locations  throughout  the  payload. 

The  force  gages  will  be  used  principally  to  measure  the 
forces  transmitted  across  the  payload/STS  interface. 

Finally,  the  thermal  measurements  would  be  used  to  measure 
temperatures  at  critical  locations  on  the  structure  and 
to  measure  thermal  fluxes  incident  at  various  locations  on 
the  payload  as  a result  of  the  sun,  as  well  as  from  other 
payloads  and  the  STS  itself. 


(Figure  8) 


Problems  currently  being  encountered  by  DATE  are  outlined 
here.  The  problem  of  payload  manifesting,  that  is,  which 
payload  will  be  launched  on  which  mission,  impacts  our 
program  because  for  any  particular  mission,  the  STS  data 
recording  system  must  be  properly  configured  for  the 
instrumentation  complement  that  is  contained  on  a specific 
payload.  Changing  the  payload  to  a different  mission 
creates  obvious  problems  of  coordination.  As  an  example  of 
the  problems,  the  payload  planned  for  mission  number  4 
has  changed  three  times  in  the  last  two  months.  While 
obviously  not  an  insurmountable  problem,  it  does  create 
headaches. 

The  schedule  is  rapidly  becoming  a problem.  Despite  the 
fact  that  the  STS  launch  schedule  has  slipped  somewhat,  our 
ability  to  start  the  necessary  preliminary  activities  has 
also  been  delayed  because  of  delays  in  funding  authority. 
This  is  creating  a schedule  compression,  which,  while  not 
yet  a serious  problem,  will  rapidly  become  one  if  funding 
commitments  are  not  soon  forthcoming. 
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(Figure  9) 


The  DATE  funding  requirements  are  outlined  here,  and  include 
funds  for  the  acquisition,  preparation  and  calibration  of 
the  flight  measurement  system  including  the  transducers, 
cabling,  and  signal  conditioning  equipment.  Additionally, 
funds  are  included  for  reduction,  analysis,  and  reporting 
of  all  data  obtained  by  the  DATS  program. 

Expenses  associated  with  the  integration  of  DATE  hardware 
onto  the  payloads,  a relatively  modes’;  expense,  would  be 
borne  by  the  payloads  themselves.  The  expenses  of  integrating 
DATE  hardware  into  the  STS  itself  would  be  borne  by  the 
STS  program  office. 

It  must  be  emphasized  that  the  funding  outlined  here  does 
not  include  the  development  of  a data  recording  system  and, 
in  fact,  assumes  essentially  free  usage  of  an  existing  system. 
Additionally,  the  costs  associated  'vith  the  developmeit  of  the 
associated  new  technology  are  not  included  and  would,  therefore, 
have  to  be  funded  through  the  NASA  Research  and  Technology 
Operating  Plan  (RTOP)  process. 

Some  continuing  measurement  prograi.:  beyond  Phase  II  would 
probably  be  desired  at  a significantly  reduced  scope  in  order 
to  evaluate  vehicle  modifications  ji.'id  the  like.  Funding  for 
such  an  effort  is  not  included  in  the  currently  requested 
resources , however . 

(Figure  10) 


DATE  is  certainly  not  the  only  program  in  existance  that  has 
as  one  of  its  objectives  the  measurement  of  STS  payload 
environments.  In  closing,  then,  a brief  review  of  the  status 
of  such  programs  was  felt  to  be  in  order. 

First,  DATE  is  officially  unfunded  for  FY  79  and  beyond.  While 
some  indications  of  potential  FY  79  funding  have  been 
received,  no  dollars  have  thus  far  been  received. 

The  Payload  Wideband  Data  System  (PWDS)  is  an  effort  by  the 
Shuttle  Payload  Integration  and  Development  Project  Office 
at  NASA/ JSC,  to  provide  payload  associated  environmental 
measurements  after  the  STS  becomes  operational.  To  the 
best  of  our  knowledge,  this  effort  remains  unfunded  at  this 
time. 
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Similarly,  Martin-Marietta  has  proposed  development  of 
the  Environmental  Response  Instrumentation  System  (ERIS) 
to  the  Air  Force  for  the  measurement  of  STS  payload 
environments,  but  again,  to  the  best  of  our  knowledge  the 
effort  remains  unfunded. 

Two  programs  that  will  make  STS  payload  environmental 
measurements,  at  least  on  a limited  basis,  are  LDEF/SBEM 
and  VFI. 

The  Shuttle  Bay  Environment  Monitor  (SBEM)  will  make 
reasonably  extensive  acoustic  and  vibration  measurements 
during  the  flight  of  the  NASA-Langley  Long  Duration  Exposure 
Facility  (LDEF).  This  will  be  limited,  however,  to  a single 
mission  on  what  is  probably  a not-too-typical  STS  payload, 
thereby  limiting  the  application  of  the  acquired  data. 

Meanwhile,  the  Verification  Flight  Instrumentation  (VFI)  is 
being  developed  by  NASA/Marshall  Space  Flight  Center  for  use 
on  the  first  two  STS  spacelab  missions.  While  here  also  the 
dynamic  and  acoustic  instrumentation  is  extensive,  it  must  be 
noted  that  its  use  will  be  limited  to  two  missions  and  much 
of  the  data  acquired  will  apply  to  rather  specialized 
Spacelab  hardware. 

In  conclusion,  therefore,  it  must  be  emphasized  that  there 
are  no  funded  programs  which  have  as  their  objectives  the 
general  characterizations  of  environments  for  the  various 
classes  of  typical  STS  payloads.  It  is  felt  that  this  situation 
should  be  remedied  quickly,  if  the  somewhat  haphazard  and 
after-the-fact  methods  that  sometimes  accompanied  the 
characterization  of  payload  environments  on  conventional 
launch  vehicles  is  to  be  avoided. 
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Figure  2 


PRINCIPLE  OBJECTIVE 

To  Develop  Accurate  Prediction  Techniques 
For  Payload  Environments  Through  The 
Iterative  Process  of  Prediction^  Measurement^ 


And  Refinement 


Figure  3 
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APPROACH 


0 ACQUIRE  BASELINE  MEASUREMENTS  ON  STS  PAYLOADS 

considering: 

PAYLOAD  MASS  AND  CONFIGURATION  EFFECTS 
BAY  LOCATION  AND  ATTACHMENT  METHOD  EFFECTS 
FLIGHT-TO-FLIGHT  VARIATIONS 

0 DEVELOP  IMPROVED  ANALYTIC  METHOD  FOR  PAYLOAD 
ENVIRONMENTAL  PREDICTION 

0 TEST  IMPROVED  METHODS  AGAINST  MEASURED  DATA 
0 REFINE  AND  ITERATE  METHODS  AS  APPROPRIATE 
0 DISSEMINATE  PREDICTION  METHODS  AND  DATA  FOR  APPLICATION 
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0 TO  CHARACTERIZE  THE  STS  PAYLOAD  DYNAMIC/  ACOUSTIC 
AND  THERMAL  ENVIRONMENTS 

0 TO  DEVELOP  INTERIM  DESIGN  AND  TEST  CRITERIA  FOR 
STS  PAYLOADS 


504 


Figure  5 


PHASE  I INSTRUMENTATION  PLAN 
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Figure  9 
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Ai;  IMPEDANCE  TECHNIQUE 
FOR 

DETERMINING  LOW  FREQUENCY 
PAYLf  AD  ENVIRONMENTS 


Kenneth  R.  Payne 

Martin  Marietta  Corporation 
Denver  Division 


SUMMARY 


Analysis  of  various  payload  configurations  is  a very  lengthy  and 
expensive  task*  Individual  extensive  raodels  of  the  payload  and  boos- 
ter must  be  mathematically  coupled  for  the  final  system  models  and 
then  elaborate  time  domain  response  analyses  conducted*  The  task  of 
the  analysis  integration  for  all  the  organizations  involved  as  well 
as  the  manpower  involved  in  model  coupling  and  loads  computations 
drive  these  costs  skyward.  With  the  expected  payloads  exhibiting  num- 
erous variations  in  configurations  and  experiments,  the  detailed 
approach  of  the  past  will  not  be  a viable  cost-effective  technique. 

A preferable  technique  would  eliminate  the  necessity  of  creating 
detailed  coupled  models  as  well  as  eliminating  the  need  of  an  integra- 
tion task.  If  possible,  the  technique  would  allow  the  payload  organi- 
zation of  designers  and  dynamicists  to  generate  with  model  information 
from  the  booster  organization,  their  own  payload  response  and  loads 
predictions* 

This  study  reported  in  this  paper  was  conducted  to  deteiinine  the 
feasibility  of  a new  impedance  technique  for  determining  payload  low 
frequency  environments.  By  accounting  for  the  dynamic  coupling  of 
the  payload  and  booster  in  the  equation  of  motion  in  the  frequency  do-  \ 
main,  the  analytical  effort  is  diminished  by  eliminating  the  final 
eigensolutions  as  well  as  reducing  the  equations  to  simple  complex 
transfer  function  multiplications.  In  addition,  the  model  require- 
ments of  the  booster  consist  of  free-free  unloaded  interface  irodal 
characteristics.  Therefore,  the  task  of  integrating  the  loads  analy- 
sis can  he  acc  ompl  i shed  h\  (^hlnininp  a t of  ''siancliircl*'  booster 
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model  data  and  the  payload  organizations  computing  chr^ir  ovn  loads 
analysis  cycles. 

The  information  presented  in  this  paper  includes  results  of  the 
use  of  the  impedance  technique  on  Titan  flight  data  as  laell  as  pre- 
dictions of  the  low  frequency  environments  for  a proposed  Shuttle 
payload.  The  requirements  for  implementing  the  in^edance  techniques 
and  it*s  feasibility  are  discussed. 
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THIS  GRAPH  SHOWS  A SCHEMATIC  OF  THE  M.ATH  MODEL  OF  PAYLOAD  PL-13A.  NODES 
781  THROUGH  786  ARE  MASS  POINTS  ON  THE  SPACECRAR.  THE  DAMPER  IS  ATTACHED 
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SPACE  SHUTTLE  PAYLOAD  LOAD  ALLEVIATION 
USING  BILINEAR  LIQUID  SPRINGS 


R.  G.  Huntins^on/  R.  E.  Martin**  and  W.  M.  Dreyerf 
General  Dynamics  Convair  Division 


l1  Dynamics  Convair  Division 

San  Diego,  California  ^ ^ ^ 


ABSTRACT 

A method  has  been  developed  for  attenuating  payload 
response  to  Space  Shuttle  Orbiter  landing  loads,  ^he  attenua- 
tion system  consists  of  preloaded  bilinear  liquid  springs  acting 
between  the  payload  and  Orbiter.  The  spnngs  provide  high 
stiffness  during  liftofT  when  Shuttle  input  loads  are  of  low  fre- 
quency and  moderate  magnitude.  During  landing,  when  Or- 
biter-induced  loads  are  high*  the  spring  preload  is  exceeded, 
lowering  the  system  frequency  and  attenuating  payload 
response  to  the  higher-frequency  input.  A nonlinear  analysis 
procedure  is  employed  and  results  are  correlated  with  test 
data  obtained  from  a one-third  scale  dynamic  model. 

INTRODUCTION 

Two  important  goals  of  the  Space  Shuttle  program  are  to: 

1 . Provide  the  simplest  possible  interface  with  a multiplicity 
of  payloads. 

2.  Accept  existing  space  payloads  with  little  or  no  redesign. 

From  the  standpoint  of  structural  loads,  these  goals  have  been 
found  to  conflict,  particularly  for  the  larger  payloads  to  be  ear- 
ned beyond  low-earth  orbit  by  the  Interim  Upper  Stage  (lUS). 
It  IS  desired  that  the  I US  be  designed  to  be  compact  ^allowing 
maximum  payload  sue)  and  that  all  payloads  attach  only  to  a 
standard  interface  on  the  lUS  (not  directly  to  the  Orbiter). 

D/namic  studies  have  shown  that  cntical  design  loads  for 
mo8»^  spacecraft  will  result  from  either  the  Space  Shuttle 
lifto  f transient  or  from  landing  of  the  Orbiter  after  an 
aborted  mission.  Further,  normal  ipitch -plane)  load  factors  on 
the  payload  substantially  exceed  the  2 to  5g  typical  of  current 
launch  vehicles  and.  in  some  cases,  the  allowable  attachment 
loads  on  the  Orbiter  are  exceeded. 

This  paper  describes  a bilinear  liquid  spring  for  attaching 
the  lUS  to  the  Orbiter  that  provides  the  proper  dynamic 
characteristics  to  attenuate  both  liftoff  and  landing  loads  to 
levels  consistent  with  current  expendable  launch  vehicles.  The 
nonlinear  dynamic  analy.sis  approach  and  its  validation 
through  dynamic  model  tests  are  also  presented, 

DYNAMIC  ENVIRONMENT 

Figure  1 illustratej  a hypothetical  spacecraft  10  feet  m 
diameter,  40  feet  long,  and  weighing  5 000  lb  that  is  typical  of 
one  class  of  payloads  requinng  an  lUb  mr  boost  to  higher  or- 
bit. The  potential  for  high  normal  accelerations  of  the 
spacecraft  and  high  local  loads  at  the  lUS  attachment  'loints 
IS  obvious  from  the  geometry. 

Both  critical  design  conditions  for  payload  major  struc- 
ture and  the  attachments  are  transient  loading  conditions 
Liftotf  ioaus  result  from  worst -case  thrust  buildup  rates  and 
thrust  differentials  among  the  three  Orbiter  liquid  rockets 


• Chief  of  Structural  Dynamics.  Member  A lA  A 

**  Manager  of  Structures  Technology,  Associate  Fellow  AlAA 
t Senior  pyoamics  Engineer 


OnilTERAnACHMCIIT 
ACCEL  TIME  HISTORIES 

Figure  I.  Trend  model. 


and  the  two  large  solid-propellant  booster  rockets  as  the  com- 
plete Shuttle  vehicle  lifts  off  the  launch  pad.  The  second  criti- 
cal load  condition  results  if  the  mission  is  aborted  and  the  Or- 
biter  lands  with  a complete  payload  at  the  design  sink  rate  of 
10  fps.  The  analytically  predicted  shock  spectra  at  the  lUS-to- 
Orbiter  attach  points  for  these  conditions  are  shown  m Figure 
2.  A shock  spectrum  is  a plot  of  the  peak  response  of  a single- 
degree-of -freedom  system  to  the  load  transient  as  a function 
of  the  system's  natural  frequency.  .Thus,  payload  mounting 
natural  frequencies  m the  4 to  5 Hz  region  will  cause  high 
dynamic  responses  for  the  liftoff  conditi.  ''^ile  the  17  Hz 
region  :s  critic?*.!  for  abort  landing. 

Trend  studies  were  performed  using  l..  .iiple  mode!  of 
Figure  1 to  determine  potential  load  levels  as  a function  of  the 
mounting  stiffness.  Due  to  the  large  overhang,  most  of  the 
response  c'«»mes  from  the  fundamental  pitching  mode.  Figure  3 
shows  the  peak  normal  acceleration  at  the  spacecraft  center 
of  gravity  as  a function  of  the  fundamental  pitching  frequen- 
cy. 

It  IS  seen  that  for  liftoff  a frequency  greater  than  7 Hz  is 
required,  but  this  frequency  range  produces  unacceptably 
high  accelerations  at  landing.  Conversely,  for  landing,  a fre- 
quency hi  w .5  Hz  IS  desired  but  gives  unacceptable  liftoff  ac- 
celerations A mounting  frequency  well  below  5 Hz  is  ruled  out 
due  to  large  deflections  at  the  tip  of  the  spacecraft  and  the  po- 
tential for  coupling  with  Shuttle  control  system  or  propulsion 
.\vstem  modes.  Frequencies  above  the  17  Hz  region  are  not 
practical  due  to  the  excess  structural  weight  needed  to  achieve 
that  degree  of  stiffne.ss.  Thus,  the  .Shuttle  dynamic  environ- 
ment requires  different  mounting  stiffnesses  .n  liftoff  and 
abort  landing  to  hold  payload  accelerations  to  the  desired  low 
levels 
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Ftfiurp  3.  Spacecraft  c^peak  acceleration  due  to  liftoff  and  land- 
ing cernus  fundamental  pitch  frequency. 


LOAD  ALLEVIATION  CONCEPT 

Several  concepts  were  evaluated  for  providink?  the  dual  stiff- 
ness characteristics  required  to  minimize  payload  response. 
The  selected  design  consists  of  a stiff  lUS  structure  mounted 
to  the  Orbiter  by  bilinear  liquid  springs.  These  springs  provide 
high  stiffness  under  the  moderate  liftoff  loading  condition  and 
lov*  stiffness  during  severe  landings. 

High  structural  stiffness  is  achieved  through  the  support 
cradle  system  shown  in  Figure  4.  The  cradle  acts  as  a 
“strongback"  for  the  lUS  and  carries  loads  from  the  IDS  to  the 
Orbiter  bay  cargo  attachment  points  Figure  5 shows  the  six 
cradle*to-Orbiier  support  locations  and  the  reaction  force 
directions. 


LOWER 

CRADLE  SPArECRAFT 

SEGMENT  MOUNTING 

FORWARD  INTERFACE 

FRAME 

Figure  4.  Bilinear  liquid  spring  cradle-to-Orbiter  attachment 
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Since  landing  response  is  primarily  in  the  Orbiter  vertical 
(pitch)  plane,  liquid  springs  were  considered  only  for  the  four 
cradle-tO’Orbiter  Z attachments.  Trend  studies  showed, 
however,  that  the  aft  spnngs  contributed  little  to  dymamic 
load  reduction.  Therefore,  liquid  spnngs  were  incorporated 
only  at  the  two  forward  Z attachments 

Liquid  spnngs  were  selected  for  this  application  because 
they  are  light,  compact,  reliable,  and  have  been  space- 
qualified.  A lujuid  spnng  stores  energy'  by  fluid  compression 
and  Its  piston  orifices  provide  clamping  Liquid  spring  double- 
action IS  achieved  b>  mechanical  caging,  as  shown  in  Figure  6. 
With  this  arrangement,  the  internal  force  due  to  differential 
piston  area  keeps  the  spring  nulled  until  preload  is  exceeded 
Preloading  is  achieved  by  silicone  fluid  pressurization  at 
assembly.  Until  the  preload  is  exceeded,  no  fluid  compression 
takes  place  and  only  the  structural  stiffness  of  the 
lUS/cradle 'Orbiter  system  is  active.  VVhen  Mie  liquid  spring 
preload  is  exceeded,  fluid  is  compressed  and  :he  system  stiff- 
ness drops  to  essentially  that  of  the  liquid  spnng  Reference  1 
presents  a detailed  explanation  of  the  liquid  .spring  .support 
sy.stem. 

ANALYTICAL  APPROACH 

Dynamic  analysis  of  the  load  alleviation  system  uses  the 
method  of  normal  modes  calculated  with  an  assumed 
linearized  spring  and  incorporates  the  ipnng  nonlinearity  as 
part  of  the  generalized  force  in  the  equations  of  motion  Figure 
7 show.s  the  bilinear  force  deflection  curve  characterizing  the 
liquid  spnng.  As  pointed  out  in  Reference  2.  it  us  extremely  im- 
portant to  account  for  this  bilinear  effect  in  the  dynamic 
analysis  t rather  than  use  an  average  stiffness  value i to  deter 
mine  the  n*spon.se  accurately 
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Figure  6.  Liquid  spring  double-action  capaOility. 


Figure  7 IVS  liquid  spring  characteristics. 


For  a linear  system,  the  equations  of  motion  in  matnx 
notation  are  given  by 

|m|  |x|  + lc(  |x]  + Ik|  ;x]  - lF(t)|  (P 

where  !ml,  (cl.  and  !kl  are  the  mass,  damping,  and  stiffness 
matrices,  respectively;  while  | F(t)j  is  the  time -dependent  ap- 
plied force  vector.  The  variables  x,  x.  and  x are  displacement, 
velocity,  and  acceleration,  respectively 

These  equations  are  solved  using  normal  iwith  respect  to 
both  the  mass  and  stiffness  matrices)  modes  as  generalized 
coordinates  by  making  the  substiiutior, 

|xi  - I 4>1  !q)  i2) 

wheie  I 1 IS  the  matrix  of  modal  vectors  and  |q)  is  the  vector 
of  time -dependent  generalized  c<x>rd;nate  :.  This  yields  the 
equation  of  motion 

(q|  +l2U>l,q|  ^-(w^llq)  - (Ml- 1 ( Fa)  | (3) 

where  (Ml  is  the  generalized  mas‘.  matrix,  (o  is  modal  circular 
frequency,  f is  modal  dam^ping  <c'C(.|-)  and  the  matrix 
superscripts  -1  and  T reprc.sent  inverse  and  transpose,  respec- 
tively 

Eq  3 IS  solved  numerically  by  standard  integration  tech- 
niques for  {q}.  and  its  derivatives  as  a funct:v>n  of  time.  Dis- 
placements. velocities,  and  accelerations  are  obtained  by 
substitution  into  Eq  2 


For  a nonlinear  Fy stem,  the  nonlineanties  are  treated  as 
external  forces  and  introduced  into  the  right  side  of  Eq  3 so 
that  the  nght  side  becomes: 

- (M1-'I*I>|T|F(u*Pi  t4) 

where  F is  the  nonlinear  fc*’ce.  The  modal  parameters  (wj, 
|M)*i  and  (<l)  I are  based  on  the  assumed  linear  spring  r»-.pre- 
sented  by  the  solid  lino  in  Figure  7. 

Nonlinear  forces  are  desr;ibed  in  terms  of  forces  and 
damping  coefficients  whose  magnitude?  are  functions  of  the 
relative  displacement  oetween  two  points  on  the  structure. 
Specifically,  at  any  two  points  on  the  structure,  i and  j,  the 
total  nonlinear  forces  are  given  by: 

F,  « P (Xj  - Xj)  + C (x;  • xj) » (Xj  - Xj)  (5) 

Fj  * -Fi  (6) 

where  P ixj  • xj)  is  the  force  as  a function  of  relative  displace- 
ment between  i and  j and  C (xi  • xj)  is  the  damping  coefficient 
as  a function  of  relative  displacement  between  i and  j Ob- 
viously, since  the  nonlinear  forces  are  functions  of  displace- 
ment, a sufficient  number  of  modes  must  be  included  in  the 
analysis  to  desenbe  the  displacements  adequately. 

In  this  case,  x\  and  xj  represent  the  displacements  of 
nodes  in  the  finite -element  model  on  eitlie:  side  of  the  liquid 
spnng  and  are  elements  of  ( x } . The  noniintar  spring  force,  P. 
is  obtained  by  table  lookup  of  the  force  difference  between  the 
solid  and  dashed  lines  of  Figure  7 for  the  given  relative  deflec- 
tion, xj  - Xj.  The  liquid  spnng  has  a nearly  constant  damping 
coefficient  so  that  C (Xj  - x^)  «=  C and  the  addition  of  the  liquid 
spnng  damping  to  the  generalized  force  is  just  C-'Xj  - Xj). 
Modal  structural  damping  is  estimated  from  expcnence  and 
included  in  the  left  side  of  Eq  3.  Eq  3 with  the  right-hand  side 
modified  by  Eq  4 to  include  nonlinear  forces,  is  solved  incre- 
mentally in  time  by  the  same  numencal  irtegration  techni- 
ques as  in  the  linear  case.  At  each  tir.ie  increment,  the  non- 
linear forces  are  computed  as  a function  of  the  appropriate 
relati'’e  deflections  and  added  to  the  linear  (time-dependent) 
forces. 

The  usual  checks  are  made  :o  determine  that  a email 
enough  integration  t.me  increment  is  ised  to  get  acceptable 
accuracy  In  addition,  our  expenence  has  shown  that  better 
accuracy  is  obtained  if  the  arbitrary  lineeiized  spnng  cons- 
tant IS  closer  to  the  lower  of  the  bilinenr  .spring  constants  If  it 
IS  too  near  the  higher  spring  constant,  loo  many  modes  are  re- 
quired to  get  accurate  deflections 

Both  liftoff  and  landing  analyses  considered  only  the 
3>'mmetnc  i pitch  plane)  response.  Liftoff  forcing  functions 
consisted  of  Orbiter  engine  and  solid  rocket  motor  thrust, 
stand  reactions,  and  applied  wind  loads.  Landing  forcing  func- 
tions included  the  longitudinal  and  vertical  components  of 
mam  landing  loads  at  the  fore  and  aft  gear  attachment  points 
for  a synimetnc  landing  condition. 

The  spacecraft  analyzed  was  assumed  to  be  rigid.  40  feet 
long,  and  to  weigh  5,000  lb.  It  wa.s  cantilevered  from  the  for- 
ward end  of  the  lUS  The  lUS/cradle  system  was  16  feet  long 
and  V eighed  36,500  lb. 

Modes  for  the  spacrcraft/lUS/cradle  and  Shuttle  were  ob- 
tained separately  and  then  combined  via  mode  svTithesis.  For 
the  Shuttle,  the  lowest  48  fret-free  symmetric  modes  were 
employed  without  payload  and  uiuonstrained  at  the  cradle  at- 
tachment points  The  spacecraft- lUo/cradle  modes  consisted 
of  the  lowest  13  symmetnc  ui^vies  with  the  system  assumed 
fixed  at  the  cradle/Orbiter  atta^'Kinent  points  Each  liquid 
spring  was  linearized  to  a value  of  10'  Ib/in  for  mode  com- 
putation. as  shown  in  Figure  7 Total  system  modes  were  ob- 
tained bv  the  component  mode  synthesis  method  described  in 
References  3 and  -y 
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In  the  response  analysis,  the  f»rst  50  system  normal 
modes  were  included.  This  covered  the  frequency  range  up  to 
24  Hz  for  liftoff  and  38  Hz  for  landing.  Liquid  spring  bilinear 
stiffness  charactenstics  are  shown  m Figure  7.  The  damping 
coefficient  for  each  liquid  spring  was  assumed  to  be  a constant 
880  Ib-sec/in. 

ANALYSIS  RESULTS 

Typical  time  histones  of  spacecraft  center  of  gravity  (eg)  tran- 
sient accelerations  dunng  Orbiter  landing  are  shown  m 
Figures  8 and  9.  As  seen  in  Fip»r«  8,  maximum  acceleration 
with  no  load  alleviation  is  5.6g  0:*-^  jccurs  0.25  second  after 
touchdown.  Peak  negative  acceleration  is  -2.7g  at  t =*  0 74 
second.  Figure  9 shows  the  corresponding  respc*nse  based  on 
analysis  using  the  liquid  '.prmg  model.  It  can  be  seen  that  the 
nature  of  the  response  is  quite  different.  Maximum  positive 
acceleration  still  occurs  at  t » 0.25  second,  but  is  rec'uced  in 
magnitude  to  3.4g.  Peak  negative  acceleration  is  mcreised  to 
-3.5g  and  occurs  earlier  at  t » 0.40  second. 


Figure  3.  Landing  response  of  spacecraft  eg  without  load 
alleviation 


Figure  9 Landing  response  of  spacecraft  eg  (vtth  Imd  u//eria- 
tion. 


Figure  10  presents  corresponding  time  histories  of  rela- 
tive dellectjons  between  the  lUS  cradle  and  the  Orbiter  for- 
ward Z-attachment  for  the  linear  and  nonlinear  systems.  For 
the  linear  system  (no  load  alleviation),  the  peak  relative 
deflection  is  about  0.072  inch  (compression),  producing  a 
force  on  each  Orbiter  forward  attachment  of  72  000  lb.  With 
the  liquid  spnng  connecting  *he  cradle  to  the  Orbiter  at  the 
forward  supports,  the  peak  relative  deflection  increases  to  - 
0.42  inch,  but  the  maximum  reaction  force  at  each  Orbiter  for- 
ward Z-support  IS  reduced  to  37.100  lb.  The  liquid  spring 
preload  was  set  at  21.750  lb.  corresponding  to  0.022  inch  rela- 
l4'"e  deflection.  This  corresponds  to  approximately  2g  load  fac- 
tor Jt  the  forward  Z-?.ttachments.  The  nonlinear  analysis 
shows  ih’^t  a substantial  portion  of  the  total  response  time  is 
spent  on  the  roft  portion  of  the  stiffness  curve  during  landing. 
Liftoff  response  calculations  for  this  configuration  predict 
that  peak  reactions  never  exceed  the  liquid  spring  preload. 


Figure  10  Liquid  .spring  deflection  during  landing 


Load  and  acceleration  reductions  due  to  the  bilinear 
spnng  are  summarized  m Table  1 for  the  baseline  system.  Sig- 
nificant reductions  in  acceleration  and  attachment  load  can 
be  seen  when  the  liquid  spring  is  incorporated  in  the  landing 
analysis. 


Table  1.  Liquid  spring  effect  on  haselire  system  response. 


Rcitults  of  llvnamic  ArtBlv<iis 

Without 

With 

Itf  m 

Condition 

Liquid  Springs 

Liqi  id  St  rings 

SpttrecrwftcB 
Pruk  Pitch  izl 

Liftoff 

J.5g 

*.og 

Acceleration 

Landing 

5.6g 

3.^g 

Forward  7 
Kcurtion 

Landing 

72.WH)  Ih 

.17.1001b 

Peak  Klai«*ic  7.- 
OeniTtion  t 

Liftoff 

1 14 

1.14 

Sparrrraft  Tip  On, I 

Lauding 

1.44 

Initial  concern  about  increa.sed  ftpacecraft  defieciion  due 
to  liquid  spring  flexibility  was  found  to  be  unwarranted.  As 
shown  in  Table  1.  .spacecraft  maximum  tip  deflection  is  only 
1 44  inches.  The  relatively  low  defiection  raises  the  question  of 
reducing  tlic  liquid  spring  stiffness  for  greater  attenuation  of 
landing  loads  .A  landing  .inalysis  assuming  a 10,000  Iban. 
spring  rate  instead  of  .39. .500  Ibon  . keeping  the  preload  cons- 
tant. increased  the  spacecraft  tip  deflection  to  2.1  inches  while 
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reducing  spacecraft  eg  peak  acceleration  to  2.5g.  Liftoff 
response  is  unaffected  by  reduction  m spring  constant  since 
preload  is  not  exceeded  during  liftoff. 

While  the  5.000-lb.  rigid  spacecraft  was  the  basis  for  the 
trend  studies  and  for  sizing  the  liquid  spring  system, 
spacecraft  flexibility  was  also  considered.  An  existing  Air 
Force  spacecraft  was  analyzed  by  including  the  spacecraft 
cantilever  modes  in  the  mode  synthesis  and  then  performing 
the  nonlinear  dynamic  response  analysis  in  the  same  manner 
as  for  the  rigid  spacecraft.  Results  verified  the  effectiveness  of 
the  liquid  spring  system  in  reducing  landing  response. 

The  effect  of  spacecraft  flexibility  was  also  evaluated 
through  trend  studies  in  which  the  5,000-lb.,  40-foot-long 
spacecraft  was  assumed  to  be  a uniform  cantilever  beam. 
Figure  11  shows  the  effect  of  spacecraft  fundamental  can- 
tilever bending  frequency  on  spacecraft  eg  peak  acceleration 
due  to  landing  when  the  liquid  spnng  system  is  employed.  It 
can  be  seen  that  peak  landing  response  is  relatively  insensi- 
tive to  spacecraft  frequency. 

NONUNEAR  ANALYSIS  VERIFICATION 

To  obtain  better  definition  of  the  analytical  uncertainties  as- 
sociated with  the  Shuttle-induced  dynamic  disturbances  and 
to  venfy  the  nonlinear  analysis  procedure,  a model  test  pro- 
gram was  conducted.  A one-third  scale  beam-type  model 
(Figure  12)  representing  the  spacecraft/IUS/cradle  system 
was  built  and  tested.  Model-to-full  scale  ratios  are  listed 
below: 

Length  1/3 

Frequency  1.0 

Mass  and  stiffness  0.056 


Figure  11.  Effect  of  spacecraft  fundamental  cantilever  bending 
frequency  on  eg  peak  acceleration. 


The  0,056  stiffness  ratio  was  based  on  an  available  liquid 
spring.  lUS/cradle  support  stiffnesses  were  modeled  with  leaf 
springs.  Two  spacecraft  were  modeled,  one  stiff  and  one  flexi- 
ble. having  interchangeable  tip  masses.  This  gave  a spacecraft 
cantilever  frequency  range  of  4 to  35  Hz. 

The  model  was  supported  on  two  Ling  Model  B335 
electromagnetic  shakers  to  ciraulate  Orbiter  acceleration  in- 
puts at  the  forward  and  aft  Orbiter-to-cradle  Z-attachments. 
Three  interchangeable  forward  support  assemblies  were 
tested: 
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1 A ri^d  support  link  to  act  as  a baseline  for  attenuation 
device  comparison. 

2.  A liquid  spnng  assembly  with  an  average  spring  rate  of 
4,500  Ib/in.  having  an  adjustable  preloau  range  from  300 
to  3,000  lb. 

3.  A dual  parallel  liquid  spring  assembly  with  an  average 
spring  rate  of  1.200  Ib/in.  and  a nonadjustable  2 400  lb 
preload.  Figure  13  shows  the  measured  force/stroke 
characteristics  for  this  liquid  spnng  assembly. 


Figure  13,  Forcelstroke  characteristics  for  dual  liquid  spnng  as- 
sembly. 

Figure  14  ^hows  the  test  setup  for  the  dual  liquid  pring 
assembly  and  the  stiff  ‘spacecraft  beam  with  tip  mass 


Test  results  correlated  well  with  analycis.  Peak  responses 
were  generally  predicted  within  20%  and  almost  always  were 
conservatively  higher  than  the  test  values.  Figure  16.  a typical 
comparison  of  predicted  and  measured  vertical  accelerations, 
shows  good  correlation  between  analysis  and  test.  Our  ex- 
penence  is  that,  particularly  in  nonlinear  systems,  it  is 
difficult  to  mathematically  model  all  damping  mechanisms; 
this  results  m the  piedicied  response  being  higher  than 
measured.  Of  course,  at  specific  points  of  lower  amplitude 
(usually  near  a node  of  an  important  mode)  the  percentage  er- 
ror of  the  calculation  may  be  higher,  but  this  is  usually  of  little 
design  consequence. 


CONCLUSIONS 

Bilinear  liquid  springs  are  an  effective  means  of  reducing 
spacecraft  loads  due  to  Orbiter  landing  without  "r'reasir.g 
liftoff  response.  This  attenuation  system  also  substantially 
reduces  the  vertical  reaction  forces  acting  at  the  Orbiter /cargo 
support  points  dunng  landing. 

The  analytical  method  for  handling  spring  nonlinearity 
has  been  validated  by  dynamic  model  tests  Based  upon  cor- 
relation with  te.st  results,  analyticallv  predicted  accelerations 
tend  to  be  conservative. 

Spacecraft  peak  acceleration  response  to  Orbiter  landing 
is  relatively  insensitive  to  spacecraft  cantilever  bending  fre- 
quency when  the  liquid  spnng  system  is  incorporated. 

For  a given  liquid  spring  preload,  spacecraft  peak  ac- 
celeration response  to  Orbiter  landing  can  be  reduced  by 
reducing  the  spring  constant  at  the  expense  of  increased 
spacecraft  dynamic  deflection  Spacecraft  deflections, 
however,  are  surprisingly  small  for  practical  values  of  spring 
constant. 
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OBJECTIVE:  DEVELOP  METHODOLOGY  TO  REDUCE  PAYLOAD  CONTRACTORS  NEED  TO 
INTERFACE  WITH  LARGE  DETAILED  SHUTTLE  STRUCTURAL  ANALYSES. 
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PRCDICTIONS/CONFIDEHCE  SHOUU)  IMPROVE  WITH  EXPERIENCE 
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PAYLOAD  EXAMPLE  USING  "EXACT"  SHOCK  SPECTRA 
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FREQUENCY  IN  HERTZ 


PAYLOAD  EXAMPLE  USING  "EX.ACT"  SHOCK  SPECTRA  (CONTINUED) 
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X-DIRECTION  ACCELERATION  PREDICTIONS  ARE  BEING  FURTHER 
EVALUATED  TO  STUDY  "STEADY-STATE"  ACCELERATION  CONTRIBUTIONS 
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ENVELOPE  OF  PAYLOAD  AnACH  POINT  MAXIMUM 
Z-DIRECTION  ACCELERATION  RESPONSES  FOR 
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ESTIMATION  OF  PAYLOAD  LOADS  USING 
RIGID  BODY  INTERFACE  ACCELERATIONS 
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SAMPLE  PROBLEM 
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FOR  STAGE  0 IGNITION 
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THE  CLAW  MODEL  SHOWED  VARIABLE  RADIAL  STIFFNESS 
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STIFFNESS  CURVES  WERE  DERIVED  USING  THE  MODEL 
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FORCE-DEFLECTION  CURVES  ARE  BASED  ON  TRUNNION  CENTERED  AT 
SPHERICAL  BEARING  AND  O.OOA  INCH  DIAMETRAL  CLEARANCE. 
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THE  APPROXIMATE  FREQUENCY  OF  VARIOUS  PAYLOADS  WOULD  BE 
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WITH  25  PERCENT  HASS  AT  EACH  OF  FOUR  INTERFACES 
A 56  RESPONSE  LOAD 

AN  AVERAGE  STIFFNESS  FROM  THE  PRECEDING  CURVES 


^-1  o 

o d 

cd  o 

O fJ 


P-  c 

cd  bop. 

*0  s 

cd 

<D  p 6 

•r|  *H 

•H  3 

to 

CO 

P P CO 

•H  CO 

0 

0 P.C0 

P *H 

p 

CO  cd 

cd 

•d 

P P 

P 

0)  0>  CO 

0 X 

p X-H 

60 

P PX 

t3 

P 

3 H 

CD  • 

•H 

0 *d 

CO  P 

P 

^ P 

cd  P 

04 

cd  • 

X P 

CO 

CO  »d 

X 

•H  CO  cd 

P 0 

p 

- 0 

j-i  id 
Cd  l/)  o • 
X :3  j.:  4-» 

u o -H  a 

•H  rH  0) 
W)  > i 
♦H  c <D 
X }h  OMh 
P. 

•H 

(d  cd  o'd 
o X 

4->  to 
O -H 
IS  4^ 

O V)  O CS 
X <u  o 
to  > 0«H 
J I V) 

O D 3 Cd 
Si  fJ  Cd  ^ 
cd  0:3 
CJ  o t>0 
1>  O X ‘H 
U-H  4^ 
C *j  O ts 

<0  u o o 
« <u  c o 

O tH  rt 
«fl»W  C 00 

<1>  lU  O c 

W‘H 
• <U  ■P 
T)  <U  p a 
« O 3 
O p (U  o 
«:<  o c E 

>.*P  .H 

« G t) 
0.0  h (0 
. X di  O 

»+^  4a  4J 
O O >H 
P 'O  Cd 
(O  O 04 
C!  ^ 
Ot)  rHX 

•H  c o>  o 
fj  cd  4J  «d 
cd  cd  o 
6 tn  Si 
•htJ  3 0) 
X cd  u V) 

000:3 

$H  »H  cd  cd 
P4X  o 
o^cd  o o 

< P<4JX 


T)  OOtH  ncJ 
Cd  to  cd 
O CD  O 

fH  (O  (/)  rH 
>s*H  C >s 

Id  o cd 

04  O 04  04 
(A  to 

0)  e:  CD  fH 
X O Si  rH 
•P  p4  Cd 

to  £ g 

4i  CD  3 CO 
O M e 
•H  p 

P *0  X o 

Ct  cd  Cd  t+4 

(DOS 

U r-C  X 
P X44  r-| 

0 Cd  O «H 
P4  p4  cd 
P »H 
10  0)  3 O 
CM  X*H  CD 
POP. 
66  CViin 
P p d> 

•H  Cd  0) 

B X X CO 
3 P P CD 
CO  *H 

to  'd'O  o 
cd  CD  C p 
p cd  o 
P cd  3 

O ti*  o 

•H  p CD 
'd  p vD  p 

O to  N «P 

CO  CD 

cd  3 P 
X CO  o>  o 
•H  CD  X 
<D  3:  63 
p 4->  P *H 
Cd  H-l  ID  X 
X 

to  3 

CD  • P *H 
> CO  c3 
P - CD  CO 
3 3 P 

O X *H  fH 
nix  3 
CD  CP  CO 

X to  0) 

H ^ *H  p 


COHSTRAINT  LOADS  ON  THE  PAYLOAD  CONSIST  OF 


CO 

c 


1—4 


I 

S 

a 

1 

i 

►—4 

o 

s 

a 

a 

C-> 

o 

a 

a 

Mi-4 

zr 

•—4 

s 

o 

s 

s 

S 

LU 

i 

s 

u. 

II 

< 


693 


' RADIAL 

18150  INCH  POUNDS  MAXIMUM 


CO 

9 

r^  d 

Cd  o d 

u 4:j  o 

Q> 

•H  4H  a 

W) 

d Pi 

^ <0  <D 

O O <0  o 

Cd  ^ ^ 

OO 

OH  *H  CO 

d P.EO 

iH  Ci  C/) 

d <0  *H  iH 

X C! 

d PJ  CM 

cc$  ••t:^  C)  3 

Pi  <C  rH 

o.  p;  45  o 

4H  h4  rH 

O cd  H P4 

(X  *C4H 

O *H 

0 0,00 

^ a W)  o 

fd 

4J  C*H  • o 

O C4H  9 

O X bOrH 

rd  O Cd 

C cd  CM 

tH  P>  O 

o +j  iH 

O O rH 

a V. 

d d 

T3  O Cd  V) 

d o rH  t:: 

<u  a *H  o *H 

O Cd  CO 

O Cd  c -O 

•ri  to  ^ *H 

!3  d <D 

p»  d <0 

'td  < O 4->  d 

'J  H f=  O 

.-3  ^ Pi  rH 

d 4-  d9 

•H  fj  o cd 

doe 

a.  pu  > 

tM  cd  *H  P> 

(D  O Om 

X *H  • 

p.  d.cj  :d  e 

Cd  *0  Cd  C to 

03  <D  4J  t;^  d 

cd  fc:  ‘H  9 

O £= 

to  O rH  d 

c/)  9 

•H  rH  O d 

•O  M 4i  o X 

fd  < o 

t/i 

CJ  <U  »H  cd 

d d 4J  i~j  a 

TJ 

O Xir-i  C B 

O O Pi 

C[5 

rH  0>  d 

•H  .H  9 Of  ^ 

O 

'tJ  iH  :3  'O 

d 4->  d u 

»-3 

c;  td  i-4  >H  d 

d O cd  cd  d 

O O Cd  p>  Cd 

d 'H  *H 

•H  rH  Pi 

Pi  Pi  O Pi 

c 

p>  cd  iH  nd 

4->  tM  Pi  O O 

•H 

o rH  Oitd  cd 

O CJU  O 

rj 

•H  cd  O O 

9 O N O 

U 

Pi  «H  *^0  ‘H  fH 

cd»d  00 

Mh  -TJ  ;d  Pi 

o 4J  d • r-4 

v> 

<d  O O iH 

r-4  OX 

P>  PirH^  Cd 

XMh  o 4J  X 

O 

<d  di*H 

cd  O |5  "H  fH 

.U 

^ +J  O (/)  *9 

ix  4H  o o 

P>  UTJ  cd 

X 0*H  P> 

TJ 

CD  tC  < p< 

Cd  4->  ^ Pi  Cd 

cj 

t/)  Pi  c/)  h4 

•H  P4  fi 

o 

•H  *H  Di  CD 

d c)  o d*H 

iH 

•9  cd 

O r d O X 

>s 

P»  4H 

P H U O 

Cl 

Pi  cd  O Q> 

t0  4J  Cd  U Pi 

Cd  Pi  X«4H 

d d > o p4 

4:2  m Cd  P>  o 

mH  O Of 

u o 

4J  U XtD  <d 

V)  P>  4J 

u u dad 

(/)  4.ind  cd  d • 

cd  o cd  4->  to 

-r4  ^4  cd  (1)  C/> 

«rl 

,d  rJ  o 9 u 9 

4-4  O O to 

4J  »y)  fH  O Pi  cd 

uxx  o u 

o u o o 

o 4H  d d 

c;  Pi  9 d d-*“< 

e d-H  o 

O 0 9 

O 9 cd  4H  g 

M cd  u d o p> 

U d u o 

*H  nd  ‘H  rH  *H 

01  u 

<0 'd  a 

d P Cd 

o cd  d 9 o •»”* 

O 9 »-•  r4  d 

to  *H  cd  rH 

•H  9 9 d 

c*  1 O 

p O cd  £; 

<0  d d rH  d d 

O iH  • p *H 

v>  o O cl  O 

*f  1 to  X 

d)  *H  mH  4H  U 

Pi  rH  pt  4J  fd 

4H  9 9 

tpi  cd  'H  u tj 

r:  tj  0 9 0 

•H  P O 

o r.d‘H  E;  cd  V) 

0 9 Cd  Pi  « 

^ Pi  C O cd 

.d  Cd  o *H  ^ 

pi  4H*44 

H J^l  rO  9 P 

694 


IF  FRICTION  IS  NEtiLECTEB 


f / 


4*. 


o 

CO 

CO 

<c 


CO 

o 

o 

CO 


I 

a 


£ 


695 


If. 


/V 

f ■/ 


c 

o 

*r^ 

fJ 

u 

•H 

V4 

li. 

w 

c 

•H 

U 

<U 

f-H 

fcC 

o 

■TJ 

oj 

O 

r: 

o 

•fH 

c 

c 

:rJ 

H 


bO 

<u 

4J  p: 

4J 

C5  'H  <D 

d 

<D  nd  <;  ♦-> 

S »-5 

•H  4-> 

oo^  H 

4J  O 

U »H  -H 

*H 

cd 

3 

4->  0 > 'H 

to 

4J  5 

to  «H 

Ct5  H 4J 

'd 

O 

d to 

4J  4J  ^ 

d -H 

rt  • Cj  +J 

o x: 

‘H 

O 0)  to  TJ 

Oa  4J 

a-»^  j c: 

O +J 

U cd  d 

o d 

<D  <U  O 

O JO 

^ r-~t  rH  to 

to 

•H  W5  »0 

O /-'^ 

<y  o d 

iH  < 

d 4J  rj 

c;i  c-3  O 

T3  < 

to  s P. 

C ^ 

pH  *1-1  »H 

Oj 

O +J  O 

X 

^ P r-^  O 

4J  rH 

O p O 

•r^  4:4 

4J  T-«  M 

f;  £; 

P^  44  TJ  Csl 

•H  0 

d U p »“H 

r4  to 

.Xi  CTJ 

to 

U Ph  to 

(0  < 

4J  ‘H 

*d 

to  »H 

C rH 

•rt  +j  g 

D (D 

j;  rt  -ri  H-i 

0 <U 

■P  ^ r-«  CsJ 

p a. 

rt  p-i 

0 0) 

0 boo  d 

0 > 

d4-l 

\0  -H 

'd  *H  d 

to  4-> 

U)  P-i  to  0 

0 

43  0>  0 

< 

•H  *i:j  *d  ^d 

0) 

d d 

^ d 

0 to  4^  0 

d d 

to  t::  *H  iH 
0 0 d 

rH  bO 

TJ  0 bOrH 

0 d 

d d 

0>  *H 

0)  X6  -H 

•i<  »d 

M d 

•H 

d ,0  Tj  d 

P-i  > • 

d P^ 

0)  0 >H 

to  d • 

4h»  Ph  0 

*d  «j  u to 

•H 

d 0 cS^H'd 

rO  d 

0 rH  )S  f4  d 

^ W P4 

rH  XO*H  d 

0 ‘f-l 

d.d  rH  0 

to 

d Oito  Oi 

0)  d*«H 

0 tU 

.X!  Ojd 

•H  0 rd  0 

4.)  *H  4-A 

d 0 p<  H 0 

4-> 

d p d 0 

add 

d ‘H  0 

0 P4*H 

P-«  pO  <0  *00 

0 

P4  d <U  iH 

to  04*d 

0 0 u 

*d  P^  <.j 

♦d  *H  d to 

d 0 

d n>  4->  44  *H 

0 U to 

Opd  u P^ 

iH  <U 

rH  4.A  (D  4)  Tj 

rH  H 

X P4  fA  d 

0 r-4  rcj 

rt  o«H  d 0 

,d  d 'd 

Oi  4->  *d  *H  i-t 

H rt  d 

696 


o 


o 

•H 

4-J 

U 

•H 

U 

bt 

•H 

•H 

cr 

C 

o 

u 

tA 

*03 

O 

o 

•H 

a 

CJ 

3 

Jh 

H 


♦J 

^ cd 

^ a 

•H  :3  .H 

c:  o ^ 
oim  G 
P4  o 
moo 
ox 

U +J  rH 

o cd 
<d  ^ o 

cd  -H 

to  +J 

p:  10  -H 

4J  Jh 

TJ  C O 
Cd  o 

O G fJ 
rH  O to 

e o 

cd  "O 
•H  C <U 

*\3  cd  X 
cd  <-) 
5^^  to 
’Td 


4-i  Cd  o 

U G 
O rH 

'Td 

u o 

"OOP 
*H  cd 
X P CO 
P U 'H 
•H  »H  P 

> P to  • 

«p  0)  P 
p:  >-H 

to  c: 

o o *H  p: 
X > o 
to  r~f  o X, 
0X0 

a > 

•H  PJ  "O  'O 
cd  *H  rH  <D 

to  o to 
Cd  bO  O O 
p:  X p. 

to-r-l  <0  s 
•H  Td  *H 
Cd  "O 

p o cd  (0 

P^rH  0»0 
O rH  cd 
O 0)  >,  O 
X rj  rH 


Cd  O O 
O • X ‘H 
rH  bCf-*  P 
P O 
Cd*H  -H 
01*0  • P 

Cd  totp 
4J  O P 
PrH 

•H  -H  p 
O P O Cd 

PX, 

•H  rH 


0)  p p Cd 
> O O-H 
•H  *H  -H  *t3 

m ^ c!  cj 
M-l  C5  »H 

<u  a 

r=:  X f < *w 
H 43  4->  O 


698 


Payload  Trunnion  Contact  Stresses 


:s 


Cd 

p 

to 

•H 

p p 

o 

iH 

0 

p 

P o 

o 

u 

to 

0 Xi 

o 

to 

p 

o 

O tp 

0 

cd 

0 o 

CVJ 

o 

p 

0 

P o 

p 

P o 

p 

0 

to 

p 

to  lO 

p 

to 

w 
<U  4h 
U O 


U P. 
•H 

m . 


V) 

(D 


<D 
U 
to  U 
O (D 


O 

u 

<d 


0 P 

^ o 

Cl}  *H 

to  o 
O »H 
to  Jh  _ 

to  M-l  »t3 

(D 

4->  O 
to  O 

to  Jh 
4-» 

O p 
c\i  <U  O 
^ 'H  4-i 

c:  u 

O H T} 
u rt 

O 


o 

JQ 
ctJ  U 


P <D 
CtJ  to 


0) 
tu  A-: 
<D  w H 
u 
0} 


u 

• • 

P to 

X 

to  0 *H 

0 

0 

'll}  u to 

0 to 

p p p. 

O P 

p 

p tp 

p 0 

•H 

O P o 

tp  p 

P.  0 o 

u 

0 

P o 

p 

p 

O r;  o 

P *H 

p 

O -H  o 

p 

to 

O 

u ^ 

to 

rH  < to 

•H  P 

o; 

CM  p . 

p 

rH 

•H  > 

p. 

CU,  p 

p 

4-t  t^rH 

bO  0 P 

O O ‘iH  O 

♦H  to 

u 

P rP 

to  P 

p 

od  0 

ctj 

O 


to 

c}  p 


0 mh 
^ o 
0} 

to 

P (O  0 
0 bO 
0 pna 
0 
to 

to  u 
<4H  O 0 


cj 
•H 

0} 

p 

X 
cd  P 


o o 

p 

Mh  <d 


P p 
pup 
o 


0 4-> 

to  P O 


0 cd 
P 


cd  p 
P td 


P 'P 

<0  0 O P 


to  P 
0 o 


U 0 
P 


P 0 U 
P^'r^  Cd  *H 
O 0 


o 

o 

Cu  0 

•H 

u 

to 

u 

p< 

x: 

P 

p 

to 

p 

P 

p 

p 

0 

0 

P 

• 

•H 

p 

rP 

•H 

bO 

p 

P 

• 

p 

+ » 

p 

P 

p 

O 

•H 

o 

•H 

o 

•H 

4-J 

to 

iH 

rH 

T} 

O 

P 

P. 

0 

■H 

P 

Cd 

bO 

> 

,.o 

O TJ 

O 

p 

o 

•H 

r-l 

P 

XO 

p 

X 

P 

P 

P 

o 

0 

• H 

0 

O 

n 

P 

P 

o 

> 

CO 

r-H 

t; 

to 

n 

'O 

o 

701 


GERUS 


1— 

• 2H 

UJ 

Q. 

O 

UJ 

a 

S 

UJ 

H- 

or: 

CJ 

c 

<t 

a 

-j 

o 

QC  _l 
LU  UJ 
h-  Q 
^ O 

s: 


CO. 

LU 

«: 

CD  ns 

z o 

►-M  .• 

oo.  oc 


S !Si 

> >i  'O' 

•— 

CkS'  UJ. 
. LU  -J; 

CO 

5 g 

3 5 


« 


702 


oor^cMroiAooir<orot-i 


CO 

> 

> 

« 

H- 

CD 

o 

!5 

oc 

s 

O 

-J 

u. 

UJ 

U.' 

UJ 

UJ 

UJ 

&; 

fe 

HEAO 

1 

ATS 

OAO 

I 

h- 

•—1 

O 

»M-i 

o. 

h- 

E 

O 

Hi 

u. 

703 


TOTAL  SWPLE 


UJ 

CJ 

OL 

Lju 

x; 


Cxl 

o 

>- 

UJ 

a:: 

CO 


<£ 

CD 


CO  CO 

§ = 

^ CD 

^ a 

CD  Z 

I—  H- 

GO 


Z CO 
1^  LU 
CO  CO 
>: 


a 

o 


CO 

<a:  I— 
zr  z 
<C  LU 

H-  CO 

s 1 

CO  h- 


LU 

LU 

a 

LU 

1—4 

H-H 

o 

oe: 

z 

CO 

CO 

=3 

1— 

CD 

CM 

CD 

rH 

a 

a: 

u. 

o 

CD 

a 

CO  t- 
LU  Z 
CO  IXl  g 


<c 


-J  CO  2 

a. 


1— 

LU 

-J 

a 

CO 

<C 

<c 

C-> 

GO 

y 

p 

LU 

•— I 

oc 

QcC 

(Z 

I— 

3- 

CD 

LU 

O 

CO 

►—4 

o 

CD 

y 

z 

LU 

CO 

a 

LU 

h- 

CO 

o 

ZD 

ZD 

O- 

CO 

PQ 

<c 

00 

f— 

CD 

LU 

OC 

LU 

O 

CD 

H- 

H— 

1 

<C 

a 

z 

ZD 


a 


o 

o 


a: 

LU 


o 

c 


oc: 

o 


s 


Ou 

CD 


OC 

Q. 


■ ^ 


704 


705 


SIGNIFICANT  REDUCTION  IN  LAUNCH  TRANSIENT  LOADS  WERE  DEMONSTRATED  IN  GROUND 


> 

o 

o 

h- 

H- 

z 

**^ 

CD 

LU 

»— • 
-J 

a 

tu 

O 

C/5 

z 

-J 

o 

1— 

a 

u. 

CO 

CO 

s 

CO 

LU 

C-> 

CO 

LL. 

O 

LU 

1— 

=> 

Z 

Z 

CO 

CO 

e: 

CD 

QC 

LJU 

CO 

a 

3 

CO 

CO 

<a: 

O 3 

UJ 
<C  GO 
h-  PQ 
<C  ' =3 
O CO 


CD  <C 

►—4 

^ z 
u-  o 


706 


n 


4>'- 


:i 


. 1 


C->  00 

O h- 

CO  ^ 

ca  ^ 
z: 

<c  I— 

LU  CD 


CD 
<C  <C 


o o 

« CO 


□c 

CO 


CD 


tn 

UJ 

<-p 


CD 

CO 

CO 


CD 


CD 

a 


t 

<c 

LU 


s 


CL.  < 

3 (5 

LU 

> LU 
LU  CD 
d)  Z 


O CO 


z 

CKr 

CU 

pa 

> 

oc 

LU 

LU 

o 

h- 

IZ 

CO 

a 

o 

c 

z 

LU 

CO 

LU 

pa 

H- 

h— 

K- 

o 

h- 

CO 

h- 

<c 

o 

CO 

LU 

LU 

<c 

<c 

LU 

LL. 

1— 

pa 

o 

■5^ 

pa 

CD 

O 

cC 

o 

3 

O 

•—4 

•V; 

a 

CD 

pa 

CD 

CO 

2^ 

LU 

CU 

o 

X 

< 

O 

CD 

CD 

CO  . 

LU 

-U 

CD 

*T" 

h- 

<C 

ZD 

IMi-4 

CD 

1— 

<C 

< 

CD 

ca 

QC 

h- 

pa 

<c 

CD 

►—4 

z 

LU 

a: 

Z! 

Q 

CU 

LU 

z 

Q 

<C 

UJ 

pa 

LU 

uu 

> 

LU 

<C 

CDCI 

> 

oc 

a 

a 

LU 

LU 

CD 

O 

h-i 

H- 

LU 

> 

QC 

3 

lU 

CO 

<y> 

QC 

ZD 

3 

CO 

< 

PQ 

LU 

CO 

CD 

z 

CD 

Q 

O 

CO 

3 

CU 

<C 

z; 

pa 

CU 

CD 

O 

h— 

<L 

LU 

o 

DC 

_l 

CO 

CO 

CD 

OP 

LU 

LU 

z 

LU 

h^ 

h- 

LU 

CD 

hr 

LU 

o 

CO 

Z 

1—4 

h“ 

C 

h- 

Ou 

CD 

Ou 

CO 

h- 

LU 

CO 

o 

•— 1 

O 

CO 

CO 

ST* 

LU 

CU 

LU 

O 

z 

O 

CC 

C/0 

LU 

• 

O 

pa 

o 

CSl 

LU 

a 

O 

►-H 

o 

LU 

LU 

QC 

LU 

CD 

o:: 

a. 

-D* 

O 

s 

o 

z 

LU 

LU 

3 

LU 

•—4 

►—4 

CD 

H- 

CD 

ZI 

h- 

OP 

a 

CD 

Q. 

3 

o 

CU 

CD 

*‘Z 

q:: 

LU 

<C 

LU 

• 

pa 

LU 

CJ 

CO 

Z 

CU 

h- 

CD 

DP 

OP 

CO 

c 

§ 

CD 

CD 

CD 

h" 

o 

LU 

CD 

ca 

1— 

2c 

z 

<C 

CO 

OP 

CO 

a 

o 

»— H 

<C 

OP 

OP 

•— 4 

< 

CO 

-U 

CU 

>- 

h— 

h- 

O 

LU 

o 

CD 

o 

>- 

CO 

O 

CD 

z 

OP 

_i 

<c 

-D 

-J 

z 

LU 

CO 

O 

a. 

<C 

<t: 

o 

DP 

LU 

LU 

LU 

*?g 

1— 

CD 

CD 

CD 

<C 

pa 

ct: 

LU 

LU 

1—4 

»— I 

LU 

h- 

LU 

o 

H- 

Z 

LU 

o 

h- 

h- 

CkP 

<c 

—1 

LU 

E 

»— 4 

cc 

_l 

pa 

I 

o 

3 

p 

CD 

CD 

ill 

1—4 

<C 

.U 

h— 

z 

O 

1 

Ou 

> 

LU 

CD 

h— 

h— 

1— 

CO 

LU 

CU 

CO 

O 

CO 

O 

CD 

S 

.3 

< 

o 

C=) 

CO 

- -4 

oc 

O 

>- 

LU 

A 

• 

0 

9 

• 

707 


C 21 
Q. 

a o 


I—  « 

m 2: 

« o 

•— « »>~4 

ac 

LU 

§ t 

or: 

UJ  o 

CD  CO 

2 ^ 

-J  o 


o 

CD 


Z 

>-  =3 

C3  a 


CO 

o^: 

LU 

h- 

CO 

o 

UJ 

o 

CJ 

CO 

!sy* 

»>"-i 

CD 

CO 

<c 

dki: 

I— 

<c 

CO 

a 

CO 

1— 

LU 

rc 

CJ 

CD 

Ckl 

1—4 

O 

-J 

LL. 

UL. 

LU 

Z 

CJ 

o 

<C 

a 

CO 

LU 

on 

03 

<c 

H- 

PQ 

CO 

>- 

a 

_l 

UJ 

<ac 

o 

>- 

CO 

LU 

a 

o. 

o 

_1 

LU 

LU 

CO 

O 

P 

K- 

o 

►-4 

Qe: 

y 

LU 

CO 

CO 

CO 

o 

o 

Q- 

LU 

CO 

o 

>- 

Cu 

GO 

LU 

<c 

q:: 

=D 

708 


i 

I 


CD 

CO 

CD 

Li. 

o 

a 

— J 

►--4 

od 

UJ 

oc 

zc 

> 

709 


PONSE  OF  PAYLOAD  AT  LIFT-OFF  WAS  ABOUT  3-SIGMA  ALTHOUGH 
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I BUFFETING  RESPONSE  LEVELS  FAR  IN  EXCESS  OF  ALL  PREVIOUS  TITAN/CENTAUR  EXPERIENCE 
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SPACELAB  structur;:  assessment 
BY  COMPARATIVE  ANALYSIS  - DR.  ROBERT  L.  MANN, 
MCDONNELL  DOUGLAS  TECHNICAL  SERVICES  COMPANY 


This  paper  describes  a technique  for  determining  whether  the  Spacelab 
structure  is  structurally  capable  of  accommodating  specific  payload 
configurations.  The  process  involves  use  of  finite  element  models 
and  post  processor  programs  to  compare  mission  peculiar  loads  with 
maximum  as-run  test  loads.  Loads  analysis,  model  verification, 
capabil ity  n;atrix  construction,  post  processor  program  and  output  format 
are  described. 
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OBJECTIVE  OF  LOADS  ANALYSIS 
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GENERALIZED  SHOCK  SPECTRA  METHOD 


» 

f 

f 

■ ! 
I 

t 


B c« 

c/j 

£ Z 

Q 2 


UJ 

CO 

o 

X 

-J 

z 

H 

Q 

X CO 

X 

< 

X 

X 

z 

s 

CJ  X 

H* 

Q 

CO 

H 

CO 

UJ 

CJ 

Z X 

o 

X 

CO 

z 

X 

z 

3 X 

X 

z 

X 

»-i« 

UJ 

o 

3 

< CD 
-J  X 

h- 

X 

H 

CO 

X 

X 

< 

X 

o 

X 

U 

>! 

H 

X 

X 

a 

O) 

X 

X 

UJ 

o 

F 

< 

X 

H* 

X 

u. 

< 

o 

o 

K" 

X _J 

CO 

1- 

u. 

z 

J— 

3 

J 

»-  < 

z 

X 

< 

< 

«J 

g 

X 

o 

CJ 

> 

UJ 

< 

Q 

X 

< 

z 

CQ 

h- 

CO 

ZD 

Q 

O CO 

K 

Q 

Q 

o 

X* 

O 

> 

z 

C/D 

Q 

< 

X 

X 

z 

H 

X X 

o 

X 

UJ 

o 

o 

Z 

w 

Z 

CO  X 

CJ 

1- 

J 

H 

o 

CO 

z 

m 

o 

o 

cc 

< 

H 

O X 

< 

CO 

X 

J 

K 

CO 

Q 

X X 

X 

CJ 

X 

X 

o 

-I 

< 

X 

X 

CO  H 

•->« 

X 

UJ 

uu 

X 

o 

a 

CO 

. 

X 

o 

H* 

< 

> 

C-> 

UJ 

H 

o 

O 

X 

X X 

X 

CO 

Q 

CO 

-J 

z 

X 

X 

Q 

o 

X 

X 

X 

o 

o 

UJ 

1- 

X 

O 

-J  X 

CO 

X 

_J 

u 

X 

u 

U. 

X* 

O 

X 

< 

X 

h- 

z 

H 

z 

< 

z 

-1 

o 

a z 

M 

H- 

ZD 

UJ 

< 

< 

X 

3: 

CJ 

UJ 

X 

o < 

o 

CJ 

5 

X 

CJ 

O 

n 

J 

Z X 

c 

< 

-J 

CO 

UJ 

LU 

a 

X 

o 

1— 

X 

X 

< 

< 

CJ 

a 

Ul 

CO 

X 

CO 

< 

CO 

UJ 

X 

< 

< 

z 

> 

1— 1 

X 

CJ  X 

X 

X 

1— 

H 

U- 

X 

UJ 

X 

< X 

-J 

u 

H 

o 

u. 

X 

CO 

X 

z 

> 

X X 

< 

UJ 

< 

UJ 

o 

CJ 

O 

Q 

z 

Z 

a. 

X 

h- 

m 

X 

z 

X 

X < 

o 

o 

o 

CO 

u 

z 

a. 

CO 

3 

CJ 

O X 

z 

UJ 

1— * 

H 

< 

a 

X 

K 

> 

Q 

u 

CO 

-J 

z 

3 

X 

u 

h- 

UJ 

o 

< 

UJ 

3 

< 

Q < 

X 

z 

X 

U 

o 

<L 

X 

v-J 

UJ 

X 

-J 

X X 

H 

3 

X 

3 

X 

CO 

CO 

CO 

z 

X 

H-  h- 

X 

X 

CO 

o 

h- 

CD 

X 

3 CJ 

o 

a. 

UJ 

LL 

Z 

u 

X X 

• 

CD 

X 

UJ 

X 

O 

% 

H 

X 

< 

z X 

X 

z 

X 

CO 

X 

H 

< 

o 

CJ 

X 

C CO 

o 

t— 1 

1— 

UJ 

UJ 

X 

X 

CJ 

< 

o 

-J 

u. 

CD 

H 

UJ 

X 

o 

X 

X 

X 

-J 

o 

z 

-J 

u 

X 

o 

CO  u 

X 

o 

< 

U. 

m 

< 

X 

UJ 

< 

X 

« o 

X 

X 

X 

< 

UJ 

UJ 

X 

•-M 

o 

iU 

X 

X 

H 

X 

cc 

CJ 

Q 

u 

u. 

CJ 

-J 

X 

< CO 

z 

X 

> 

u 

< 

O 

< 

o 

CJ 

CJ 

X 

X 

o 

X 

UJ 

u. 

UJ 

UJ 

H 

X 

h-  «J 

H 

< 

u 

a: 

X 

X 

UJ 

X 

o 

a < 

X 

z 

X 

< 

UJ 

X 

h- 

H 

o 

CD 

X 

X 

X Q 

X 

Q 

< 

o 

CL 

h- 

o 

< 

z 

> 

X o 

h- 

z 

X 

<o 

Z 

< 

CD 

X 

X 

Q 

CO  z 

< 

(O 

CJ 

•-< 

UJ 

Z 

o 

X 

Q 

X 

X 

X 

< 

-J 

Ul 

Z 

CJ 

M 

o 

X 

X 

< 

UJ 

CD 

o 

UJ 

H 

o 

z 

**•4 

o 

X 

X 

CO 

U1 

X 

z 

ID 

CJ 

z 

X 

3 

-J 

O X 

z 

3 

CO 

QC 

H 

Q 

z 

CO 

< 

< 

X ^ 

o 

h- 

z 

X 

H 

O 

UJ 

X 

u 

X 

CO  CJ 

Q 

o 

o 

•im 

< 

u. 

X 

<! 

X 

CO 

CO 

X 

X 

•— « 

X 

3 

X 

1- 

o 

UJ 

h- 

UJ 

X 

X 

z 

X X 

X 

X 

CO 

z 

X 

X 

X 

o 

X 

X 

X X 

X 

H 

X 

X 

Ul 

UJ 

CO 

1- 

X 

1- 

u 

1- 

CD 

H > 

X 

CO 

X 

o 

X 

CO 

H- 

3: 

z 

<o 

1 

1 

1 

1 

1 

o 

l-M 

Q- 

CO 

CO 

X 

lU 

o 

X 

a 

766 


I 


GENERALIZED  SHOCK  SPECTRA  METHOD 


LU 

QC 

O 

LU 

<C 

* 

*J- 

LU 

oc 

UJ 

h— 

1 

UJ 

<C 

CJ 

LU 

zn 

O 

> 

» 

=D 

LU 

<c 

h- 

CD 

2 

LU 

UJ 

O 

H- 

25 

CU 

2! 

C/5 

CkC 

CU 

LU 

ZZD 

ZSl 

1— 

oc 

CU 

2 

M 

CO 

o 

LU 

< 

CL. 

LU 

CO 

O 

oo 

_l 

a 

H— 

<c 

LU 

2^ 

o 

QC 

<c 

O 

cu 

LU 

CO 

C>0 

> 

y 

cu 

LU 

G 

LU 

o 

•—4 

CU 

LU 

QC 

UJ 

CO 

C/) 

LU 

o 

— J 

<c 

Id 

CO 

Q- 

<c 

LU 

s 

LU 

s: 

a 

CO 

> 

LU 

OC 

o 

o 

1— 

CU 

o 

y 

cu 

UJ 

a 

-4J 

LU 

LU 

CO 

CU 

2 

o 

I— 

O 

ID 

<c 

UJ 

» V* 

o 

<?: 

21 

LU 

CO 

<c 

cu 

-J 

CU 

> 

QC 

CO 

o 

CO 

1— 

LU 

C>0 

f— 

CL. 

UJ 

•U 

QC 

o 

CO 

CO 

LU 

C/5 

UJ 

Q. 

_l 

<c 

oc 

o 

C/5 

<c 

h- 

cu 

UJ 

Q. 

n 

CO 

UJ 

CO 

o 

o 

o 

UJ 

cu 

LU 

4 

<c 

QC 

o 

U- 

QC 

c/5 

UJ 

<c 

H- 

a 

LU 

o 

C/5 

K-M 

2E! 

•W'. 


"S 

6 


767 


LAUNCH  VEHICLE  FORCING  FUNCTION  DECOMPOSED  INTO  IMPULSE 
DELTA  FUNCTIONS  FOR  EACH  LAUNCH  VEHICLE  MODE 
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GENERALIZED  SHOCK  SPECTRA 
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PANEL  SESSION  - QUESTIONS  AND  RESPONSES 


♦ 


Card  to  Gatto 


Do  you  feel  there  is  an  active  program  to  reduce  cost  of  integration  on 
the  Shuttle? 

Gatto: 


This  is  a broad  question.  In  the  area  that  I have  cognizance  of,  yes, 
there  is.  In  the  real  world,  it  is  not  structured  for  high  priorities  and 
funds.  To  make  tools  more  efficient  requires  service  of  our  better  people 
in  the  program.  It  is  not  easy  to  take  these  people  off  of  critical  project 
efforts  and  assign  them  to  methodology  improvement  programs.  Yahata 
is  very  active  in  this.  This  problem  has  been  brought  to  the  attention  of 
management.  Management  has  been  very  receptive.  We  are  trying  to  do 
our  part;  we  are  getting  some  efficiency  in  calculation  techniques. 

Stone: 


Rockwell  has  devoted  discretionary  resources  to  this  topic.  Use  of 
interactive  graphic  computers  was  added  to  the  design  effort  and  has  been 
successful.  We  are  trying  to  cut  down  on  turnaround  time  for  doing 
analysis.  Interfaces  have  been  set  up  for  direct  communication  between 
people  involved. 

Hill  to  Card: 


The  papers  on  reducing  loads  focused  primarily  on  landing  cantilever  loads 
in  the  Orbiter.  Might  you  be  able  to  design  an  optimal  mechanical  device  to 
give  an  index  to  performance?  Same  kind  of  device  that  has  been  done  for 
flutter  based  on  measure  of  index  for  performance. 

Card: 


In  terms  of  Langley's  work,  these  techniques  have  been  applied  in 
a practical  sense  but  there  could  be  something  new  still  out. 

Ryan: 

Thir  optimized  index  of  performance  has  been  done  in  control  work,, 
and  we  are  getting  beneficial  results  and  insights.  Using  the  work  done  on 
active  flutter  suppression  and  control  as  a basis,  we  could  and  should 
develop  performance  indexes  for  loads  on  payloads  and  go  from  there. 
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Wade  to  Herzbergt 


Is  there  a need  for  payloads  to  underscand  criteria  that  Shuttle  has  used 
in  order  to  understand  how  much  conservatism  i^.  used  in  loads  and  forcing 
functions?  There  was  one  workshop  with  payload  users  which  described 
criteria  used  on  Shuttle.  Would  you  suggest  we  continue  this  as  a normal 
service? 

Herzberg; 

Not  sure  which  workshop  you  are  talking  about.  I'm  thinking  more 
of  a straightforward  interchange  between  people  who  are  critically 
responsible  for  the  payload.  Do  not  think  a workshop  is  the  right  forum; 
just  having  access  to  technical  information  is  all  that  is  required. 

Wade; 

k 

T 

The  Structural/Mechanical  Working  Group  perhaps  would  be  a good 
forum  for  anyone  who  has  a payload  to  fly  on  Shuttle.  These  are  technical 
Interchange  meetings.  I suggest  it  be  done  there. 

Herzberg; 

In  the  Marshall  involvement  on  Shuttle  at  Lockheed,  emphasis  has  been 
placed  on  uncertainty  factors  on  Shutth.  and  payload  model.  Recent 
modal  test  results  should  help  reduce  these.  Can  we  expect  to  receive 
this  type  results  of  ground  vibration  test? 

Ryan; 

If  people  want  it,  NASA  and  MSFC  and  other  Centers  would  be  happy 
to  respond  to  the  payload  community.  Particularly,  let  Don  Wade  or  me 
know.  I'm  sure  we  can  do  this.  Information  on  dynamic  test  and  follow-on 
test  could  be  provided  or  discussed  in  another  meeting. 

Wade; 


For  Air  Force  payloads,  post  test  data  correlation  has  been  sent 
when  completed  to  the  PIC  contractor.  Aerospace,  and  SAMSO. 
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Simondi: 


You  reach  CDR  on  payloads  two  years  prior  to  launch.  Changes  after 
CDR  are  expensive  and  not  desired.  Effect  on  cargo  element  and  cargo 
dynamic  analysis  are  at  the  18  months  point.  Why  is  this  acceptable? 

Wade: 


It  is  not  until  18  months  before  the  flight  manifest  is  defined; 
consequently,  it  is  not  practical  to  run  a verification  analysis  until  the 
manifest  of  the  configuration  is  defined.  This  is  why  we  shc,v  18  months. 
The  final  verification  load  analysis  is  planned  to  be  finished  6 months 
prior  to  flight.  In  some  cases,  much  of  the  verification  of  models  for 
payloads  are  not  accomplished  until  then;  therefore,  judgement  has  been 
standard.  Special  requests  for  payload  verification  at  some  other  time 
are  run.  For  example,  we  are  willing  to  do  this  on  TDRS.  We  are 
running  somewhat  early  to  support  ground  test  program  to  be  run  at 
9 months  prior  to  flight.  We  try  to  be  flexible  within  reason,  but  we 
have  to  establish  the  manifest  first. 


Simondi: 


POOR  QUALITY 


In  final  design  analysis  just  prior  to  CDR,  if  you  don’t  have  knowledge 
of  cargo  at  this  point,  doesn’t  this  give  problems? 


Wade: 


Problem  is  one  primarily  of  interaction  of  payload  with  whatever 
else  is  riding  --not  a big  driver.  The  more  sensitive  thing  appears  to 
be  the  position  in  the  payload  bay,  which  depends  quite  a bit  on  dynamics 
the  payload.  Some  payloads  are  insensitive  to  position;  others  are 
sensitive.  You  put  more  than  one  model  in  the  payload  bay  in  different 
positions  so  the  position  can  be  checked  on  loads.  Our  standard  service 
is  to  perform  final  verification  analysis  only.  Design  analysis  of  the 
payload  description  is  not  what  you  want  run  when  you  want  it  run; 
therefore,  you  should  x*un  payloads  in  several  places  of  the  bay  with 
spectrum  of  payloads  from  very  low  to  high  frequencies. 


Pcngelley; 


In  the  past  ten  years,  there  has  been  almost  infinite  increase  in  our 
ability  to  handle  linearized  eigenvalue  solutions  for  real,  undamped 
problems.  During  this  time,  hardly  anyone  has  learned  anything  on 
nonlinear  structures.  We  can't  get  rid  of  clearances  and  joints  in 
payload  design  and  fabrication.  We  are  still  seeing  nonlinearities, 
parametric  excitation,  damping,  structural  damping,  viscous  damping. 
We  need  work  on  damping  predictions  and  nonlinearities. 


Wada; 

I have  put  together  a paper  for  the  university  on  how  to  handle 
damping.  Everybody  uses  simplified  approach--!  percent  for  modal 
coupling.  Ninety-nine  percent  use  this  approach.  Nonlinearities  cannot 
be  defined.  Different  type  of  damping  for  each  different  joint.  Almost 
nobody  can  define  input  parameters  for  damping, 

Pengelley: 

Is  it  a lack  of  knowledge  of  what  to  put  in? 

Ryant 

Everyone  has  taken  the  easy  way  out  with  1 percent.  We  need  simple 
models  generated  on  physical  insight  instead  of  large  finite  element 
systems.  We  can't  depend  entirely  on  computers,  but  must  depend  on 
physics. 

Herzberg; 

The  problem  of  damping  goes  back  as  far  as  1870;  no  progress  has 
been  made.  We  will  never  make  any  progress  unless  there  is  a large 
volume  of  experisne'ntal  data.  This  is  one  area  where  the  universities 
must  help  us.  ^ 

Craig; 


True,  the  universities  would  like  to  make  contribution. 
Hill: 


Mathematical  tools  to  model,  in  principle,  have  been  around.  To 
handle  large  eigenvalue  problt?ms,  these  tools  are  not  practical.  Computers 
can  handle  problems  without  damping,  cai\  handle  mass.  Trying  to  develop 
best  linear  model  for  damping. 


' T 
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Slmondl  and  Keegan; 


Noticed  Shuttle  is  paying  a lot  of  attention  to  liftoff  and  landing 
loads.  Talked  about  buildup  rate  of  liquids,  etc.  There  is  an  absence 
of  discussions  on  transonic  buffeting  and  how  it  might  affect  spacecraft 
loads.  Af. sr  Voyager,  we  wonder  what  thoughts  are  on  pursuing  this? 

Wade: 


Ran  combined  flutter  buffet  model  of  configuration  for  both  ascent 
and  descent.  Looked  at  buffet  levels.  Probably  biggest  concern  are 
loads  on  vertical  tall  as  opposed  to  payloads.  Load  checks  are  made  at 
transonic  Mach  number.  Transonic  loads  were  low;  therefore,  we 
thought  there  is  enough  n.argin.  Were  able  to  tolerate  buffet  levels 
we  expected.  No  significant  overall  vehicle  response  to  buffet  from 
tests. 


Gatto: 


Op 


This  is  similar  to  what  we  have  done.  During  descent,  coming 
back,  operate  with  rudder  deployed  open.  Buffeting  is  something  we 
will  have  to  wait  on  flight  experience  to  see  if  tasting  has  been  adequate. 
Liftoff  and  landing  are  where  payload  dynamics  are  of  interest  to 
Shuttle  vehicle  design.  These  are  only  two  out  of  many  flight  events 
that  are  important. 


Wade; 


Everyone  should  use  other  conditions  to  be  sure  they  e not 
overlooked  something,  particularly  high  q loading  conditions.  Loads 
man  should  look  at  overall  load  cases. 


From  the  floor: 


For  payload  performance  capability  increases,  will  you  be  going  to 
weight  reduction  or  thrust  ?ugmentation?  Do  you  foresee  affect  on 
loads  analysis  or  acoustic  environment? 


Wade; 

Don't  envision  there  will  be  very  much  effect.  Will  be  some  changes 
in  bridge  fairings  on  Orbiter  to  take  weight  out  of  bridges.  Bridge  design 
has  been  based  on  Orbiter  capability,  not  on  payload  requirements.  This 
will  allow  some  weight  reduction  in  some  of  the  bridges.  Different  primary 
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solid  rocket  boosters  are  being  considered,  where  the  thrust  for  rocket 
is  equal  to  their  weight.  Ignition  acoustics  and  overpressure  of  solids 
environment  could  be  worse  at  liftoff.  We  will  need  to  look  at  max  q 
alleviation  again.  We  may  have  to  reshape  trajectory  based  on  additional 
thrust  and  weight.  Should  not  affect  payload  world  very  much. 

Simondi  and  Keegan; 

CSFC  presentation  on  DATE  appeared  to  support  dynamics.  If  this  is 
the  case,  considering  state  of  funding,  would  the  panel  make  comments 
to  support  or  oppose  operational  measurement  system? 

Card; 

Already  cast  vote. 

Ryan; 

We  all  want  flight  data.  First  Shuttle  flights  have  instrumentation 
geared  to  remove  uncertainty  in  vehicle  parameters  as  we  move  down- 
stream. Primarily  for  getting  Shuttle  ready  for  operational  flight  first. 

We  have  requested  special  instrumentation  on  first  Spacelab  flights  for 
dynamic  data. 

Wade; 

r 

We  do  not  have  enough  flight  instruments.  Not  successful  in  getting 
everything  desired  in  payload  world.  During  OFT  flights,  biggest  void 
was  in  getting  enough  accelerometers  in  the  payload  bay.  System  proposed 
on  LDEF  has  a problem  with  Z to  50  Hz  accelerometers  which  cannot 
measure  transient  loads.  Were  able  to  get  some  0 to  50  Hz  in  payload 
bay  area.  Encourage  more  thought  be  given  to  instrumentation  on  DATE 
especially  in  payload  area.  Also,  get  0 to  50  Hz  loads  transducers  on 
LDEF  payload. 

From  the  floor: 


■y/hat  is  the  status  of  payloads  to  flight  on  first  Shuttle  flights?  Are  they 
already  designed? 


Wade; 


• First  launch  will  carry  DFI  to  check  out  Shuttle  system  itself. 
Second  flight  is  TRS  to  rescue  Skylab,  being  developed  at  MSFC.  CDR 
on  TRS  is  this  week.  Third  flight  will  contain  OSTA-1.  Experiments 


are  now  being  integrated  into  OFT  payload  through  CDR.  Payloads 
are  now  in  development.  Payload  integration  plans  are  being  written. 

From  the  floor; 

Are  there  any  guidelines  to  reduce  loads? 

Wade: 


Frequencies  that  stay  away  from  2.  5 Hz  system. 
Is  there  a publication  of  these  guidelines  forthcoming? 
Wade: 


Two  attempts  were  made  to  do  this,  MCR  1612  and  a study  with 
the  Air  Force  on  load  alleviation.  We  put  out  some  guidelines  which 
Boeing  is  using. 
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